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QUASI-ONE-DIMENSIONAL COMPRESSIBLE FLOW ACROSS 
FACE SEALS AND NARROW SLOTS 
II - COMPUTER PROGRAM 
by John Zuk and Patricia J. Smith 
Lewis Research Center 

SUMMARY 

The computer program presented in this report calculates the properties of com- 
pressible fluid flow with friction across shaft face seals and throi^h narrow slots. The 
mathematical model on which this program is based is developed in a companion report. 
A summary of this analysis is included. 

The program requires the seal geometry, the reservoir conditions, and the gas 
properties as input. It will then calculate mass and volume leakage across the seal; 
mean friction factor; force; center of pressure; and distributions of pressure, temper- 
ature, density, friction parameter, velocity, and Mach number across the seal for both 
laminar and turbulent flow regimes and for choked and nonchoked flow. 

The program is written in FORTRAN IV. Input, internal calculations, and output 
are in either the International System of Units (SI) or the U. S. customary system. Au- 
tomatic plotting of several of the output parameters is an available option. 

This report includes a description of the program, details of the mmierical tech- 
niques required to solve the set of equations in the mathematical model, flow charts of 
the main prc^ram and the essential subprograms, a complete description of the input 
variables, and sample input and output. 


INTRODUCTION 

The computer prc^ram presented in this report, QUASC, calculates the properties 
of compressible fluid flow with friction across shaft face seals and through narrow slots. 
It is based on a mathematical model which is developed in a companion report (ref. 1). 

A summary of this analysis and the governing equations of the model are given in the 
next section. 


The analysis includes fluid inertia, viscous friction, and entrance losses. The 
model is valid for both completely subsonic flow and choked flow. It is also valid for 
both laminar and turbulent flow regimes. 

The computer program has proved useful in seal design, where much of the phys- 
ical information of interest is difficult to determine experimentally. A mathematical 
model tested on a computer can generate this needed information. The model on which 
QUASC is based has been shown to be a good model for conventional face seals which 
are pressure (force) balanced and for seals with self-actii^ lift pads. These two types 
of seals are shown in figures 1 and 2. 

Some of the physical parameters of interest in designing a seal are the leakage, the 
pressure distribution across the seal, and the opening (separating) force. These and 
other parameters are determined by the prc^ram QUASC for specified seal geometries, 
reservoir pressures and temperatures, and gas properties. The program also requires 
two additional parameters. One of these is the variation of a mean Fanning friction 
factor with Reynolds number. The other is the entrance loss coefficient. 

Results of the computer analysis agree quite well with experimental results. The 
comparison of these results is discussed in detail in reference 1. 

This report is intended to serve three purposes: (1) to give a summary of the quasi- 
one-dimensional analysis of compressible flow with friction; (2) to give a detailed de- 
scription of the computer program QUASC which performs this analysis; and (3) to serve 
as a user's guide for QUASC. 

SUMMARY OF QUASI-ONE-DIMENSIONAL ANALYSIS 

The mathematical model on which the program QUASC is based is for face seals 
which consist of two parallel, coaxial, circular rings separated by a very narrow gap 
(fig. 3). A pressure difference exists between the rings' inner and outer radii. The 
cavities on either side (i. d. and o. d. ) of the sealing dam are assumed to be constant 
pressure reservoirs. 

The flow is assumed to be entirely radial. It has been shown that the effects of 
rotation of the two seal faces are negligible for most applications (ref. 2). Since the 
radial distance across the seal is very small compared to the mean radius of the rings 
(i. d. /o. d. 0. 98), the area expansion is negligible. The analysis is also valid for flow 
throi^h a narrow slot. 

In reference 1, the analysis is sej^rated into two parts. The first part is the flow 
in the passage itself, where the flow is assumed to behave as a constant-area adiabatic 
flow with friction. The second part is the flow in the entrance r^ion. These two parts 
are discussed separately herein also. 
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Constant-Area Adiabatic Flow with Friction 


It is assumed that the flow in the seal leakage flow region behaves as a constant- 
area adiabatic flow with friction. A quasi-one-dimensional approximation is made 
wherein it is assumed that the flow properties can be described in terms of their cross- 
sectional averages. 

The followii^ assumptions have been made in the analysis: 

(1) The area expansion due to radius increase is neglected. (In most main-shaft 
face seals the i. d. /o. d. is about 0. 98. ) Hence, the analysis is also valid for narrow 
slots. 

(2) The flow is adiabatic. 

(3) No shaft work is done on or by the system. 

(4) No potential energy difference is present such as caused by elevation differences. 

(5) The fluid behaves as a perfect gas. 

(6) The sealing surfaces are parallel. 

With these assumptions, the flow is commonly known as Fanno line flow. 

In reference 1, under the model assumptions, the governing conservation equations 
reduce to the following (all symbols are defined in appendix A): 

Conservation of mass: 


5e + l*!!=o 

P 2 ^2 

Conservation of momentum: 

-A dP - dA^ = M du 

Conservation of energy: 


Perfect gas relation: 

P p T 


( 1 ) 


( 2 ) 


(3) 


( 4 ) 
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With the introduction of the mean Fanning friction factor f 


Y_ -dP/dx 
2pu^/D 


and the hydraulic diameter D 



the conservation equations can be combined into one equation 

^ yM^ ^ 2 ~ 

Integrating equation (7) from the point of interest to the point of choking gives 






(1 - M^)dM^ 
jl + 1 (y - 1 )m2 


( 5 ) 


( 6 ) 


(V) 


( 8 ) 


(Since eq. (7) is a first-order differential equation, one boundary condition must be 
specified. The only known condition is that Mach number is unity at the critical choking 
length L (fig. 4). ) Carrying out the integration gives 


^{h-x) = 
D 


1 1 - 
^ m2 


+ 2Li_l 

2y 


In 


- (y + 1 )m 2 

J 

1 + i(y - 1 )m2' 


= B(M) 


(9) 


The mean Fannii^ friction factor does not vary with flow length and is known to vary 
with Reynolds number according to the relation (ref. 1) 



Re 


n 


( 10 ) 
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The hydraulic diameter for both coaxial rings and parallel planes (narrow slots) is 
twice the film thickness, 


D = 2h 


( 11 ) 


Further manipulation of the conservation equations and the perfect-gas relation re- 
sult in the following equations: 


P 


- [l + (y - l)M^dM^ 


M' 


1 + ^ (y - l)M^j 


( 12 ) 


- (y - l)dM^ 
dT 2 


T 1 + 1 (y - 1)M^ 
2 


(13) 


which, when integrated from the point of interest to the point of choking, give 


^ = _L 

P* M 


- (y + 1) 
2 


1 + -^ (y - 1)M^ 


1/2 


(14) 


- (y + 1) 

2 


1 + - (y - 1)M^ 
2 


(15) 


Other equations that are needed are those for Re 3 molds number and Mach number: 

PjUjD 


Re = 


Ml 


(16) 


c yy^T 


(17) 
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Entrance Flow 


As stated in reference 1, the entrance region flow can be assumed to be isentropic 
and adiabatic. Since real entrance flows are not isentropic because of viscous friction, 
turnii^ losses, and so forth, the entrance loss is accounted for by introducii^ an em- 
pirically determined entrance velocity loss coefficient into the isentropic flow equa- 
tions. The resulting entrance temperature equation is 


T 


1 “ 


T 


0 


(r - m\ 
1 + 

2CI 


(18) 


while the entrance pressure equation is 


^1 = 


2iy/(r-i) 

(y - l)Mj^ 


1 + 


2C 


L J 


(19) 


When = 1> ih© How is isentropic. 


Solution of Equations 

The equations developed in the preceding sections, plus the equation of state written 
as 


R = ( 20 ) 

P 

must now be solved. The independent variable is chosen to be x, the distance across 
the face of the sealing dam (fig. 3). The known parameters are 

Pq sealed gas pressure (upstream reservoir pressure) 

P„ ambient pressure (downstream reservoir pressure) 

0 

Tq sealed gas temperature (upstream reservoir temperature, also the stagnation tem- 
perature) 



Cj^ entrance velocity loss coefficient 
AR distance across face of sealing dam 
h film thickness 

k, n constants in friction factor - Rejmolds number relation 

These two constants apply to laminar flow (Re :< Re^) and turbulent flow (Re > Re^. 
For flow in the transition region (Re^ < Re < Re^.), variation of the friction factor with 
Reynolds number is determined by the method described in appendix B. 

The set of equations in the mathematical model cannot be solved explicitly. Con- 
sequently, an implicit solution, that is, an iterative procedure, must be used. 

There are three types of flow considered: (1) critical flow, when the exit Mach 
number is 1 and the exit pressure equals the ambient pressure; (2) supercritical flow, 
when the exit Mach number is 1 and the exit pressure is greater than the ambient pres- 
sure; and (3) subcritical flow, when the exit Mach number is less than 1 and the exit 
pressure is equal to the ambient pressure (flow is entirely subsonic). Since the bound- 
ary conditions on the mathematical model are slightly different for each type of flow, 
the solution of the equations is slightly different for each type. The simplest set of 
equations is for critical flow. Hence, this flow case is considered first. 

Critical flow . - For critical flow, the boundary conditions are as follows: 

Exit pressure 1?2 = P3 

Choking pressure P* = Pg 

Exit Mach number Mg = 1 

Flow length L = AR 

For this case it is necessary to calculate choking film thickness h*, entrance Mach 
number Mj^, and mean friction factor f. 

Start the implicit solution of the equations by guessing a value for Mj. Substitute 
for Mj in equations (14) and (19) giving 




1 

- (y + 1) 
2 

p* 

Mj ^3 

Ml 

1 +-(y - l)Mf 

L 2 \ 


1/2 


( 21 ) 
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and 




1 


^1 + 1 (y - 1 )(Mi/Cl)^] 


2]y/(y-i) 


( 22 ) 


Dividing equation (21) by equation (22) gives 


P1/P3 

Pl/Po 


p / [l+l(y - 1 )(M^/Cl)^ 

/y + 1 1 L 2 ^ ^ -I 

P„ \ 2 M. r 1 21^/2 

3 1 |l+l(y-l)Mjj 


(23) 


Equation (23) can be solved for Mj by the method of iteration, a standard numerical 
analysis technique (see appendix C), It should be noted that equation (23) does not have 
a solution for certain combinations of Pq, Pg, and y. In that case, the flow is as- 
sumed to be subcritical. See appendix B for further discussion of this point. 

Once Mj^ has been found, it can be substituted into equation (9) giving 

(24) 

Since x = 0 at the entrance, the flow lei^th L is e^al to AR. In equation (24), both 
f and h* are unknown. However, it is known that f = k/Re (eq. (10)). Substituting in 
equation (23) for f, Reynolds number (eq. (16)), and hydraulic diameter (eq. (11)) gives 

(25) 

^ h* \pjU^2h*y 

Entrance density pj is calculated from the perfect-gas relation (eq. (20)) 





where Tj and Pj are calculated from the entrance equations (18) and (19). Leakage 
flow velocity Uj is calculated from the definition of Mach number (eq. (17)). Equa- 
tion (25) can now be solved for h*. Iteration is required to determine whether the flow 
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is laminar, transition, or turbulent. A block diagram of this solution is shown in fig- 
ure 5. 

When Mj and h* have been found, it is then possible to calculate the distribution, 
across the face of the seal, of pressure, temperature, Mach number, velocity, density, 
and friction parameter B. First, calculate T* from equation (15) for M = Mj. 

Choose a value of the independent variable and solve equations (1), (9), (14), (17), 
and (20) simultaneously. Since the solution cannot be obtained explicitly, an iterative 
solution must again be used. First, substitute for x in equation (9), giving 

= B(M^) (26) 

Solve equation (26) by the method of iteration for M^. Once M^ has been determined, 
find from equation (15). Substitute and in the definition of Mach number 
(eq. (17)) to give u^. Substitute u^, Uj, and pj into the conservation-of-mass equa- 
tion 

^I'^l = ^x"x (27) 

to give Determine from the perfect-gas relation (eq. (20)). 

Supercritical flow . - For supercritical flow, the film thickness h is greater than 
h*. The boundary conditions for this type of flow are as follows: 

Exit Mach number Mg = 1 

L = AR 

Again, the model equations cannot be solved explicitly. Start the implicit solution by 
guessing a value for Mj. Calculate and Pj from the entrance equations (18) 
and (19). Calculate pj from the perfect-gas relation (eq. (20)) and Uj from the def- 
inition of the Mach number (eq. (17)). Calculate the Re 3 molds number from equation (16). 
Use this Reynolds number to determine the friction factor 7. Substitute for h and F 
at X = 0 in equation (9), which gives 


4f AR 
D 



(23) 


Substituting the value of Bj into equation (9) and solvii^ for M by the method of itera- 
tion gives a new value for Mj. If the new and old values of Mj are equal, the set of 
equations is solved. H not, use the new value of Mj as the guess for Mj and repeat 
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the procedure until the two values for are equal. A block diagram for this scheme 
is shown in figure 6. The distribution of pressure, temperature, Mach number, veloc- 
ity, density, and friction parameter B are found in the same way as for critical flow. 

Subcritical (subsonic) flow . - For subcritical flow, the film thickness h is less 
than h*. For this type of flow, the boundary conditions are as follows: 


L = Li2 + AR 


^2 = ^3 


where Lg is the distance from the end of the seal to the imaginary point of choking. 
The additional unknown parameter Lg, which depends on Mg, is zero in the two other 
types of flow. The flow lei^hs involved in this type of flow are shown in figure 4. 

As before, the set of equations in the model for subcritical flow cannot be solved 
explicitly. The implicit solution is similar to that for supercritical flow. Start by 
guessing M^; then calculate T^, Pj, p^, Uj, Re, and f. Now, however, calculate 
= B(Mp from equation (9). Calculate L from equation (9) at x = 0 


BjD 


(29) 


From the second boundary condition, calculate Bg, 

Lg = L - AR 


Calculate Bg from equation (9) 



(30) 


(31) 


Substitute Bg into equation (9) 

Bg = B(Mg) (32) 

and solve equation (32) for Mg by the method of iteration. Calculate P for Mj^ from 
equation (14) and Pg from equation (14) for P and Mg. If Pg and Pg are equal, 
that is, if the second boundary condition is satisfied, the set of equations is solved. If 
they are not equal, a new Mj must be guessed and the procedure must be repeated. 

The successive values of M^ and Pg are saved, and the correct value of Mj is de- 
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termined graphically for P 2 = Pg. On the computer, a curve-fitting scheme is used to 
determine the correct Mj from the stored values of Mj and Pg. Figure 7 shows a 
detailed block diagram of this calculation scheme. The distributions, across the seal 
face, of pressure, temperature, Mach number, velocity, density, and friction param- 
eter are determined the same way as for critical flow. 

COMPUTER PROGRAM 

The prc^ram QUASC which performs the analysis of quasi-one-dimensional com- 
pressible flow with friction across shaft face seals and narrow slots is written in FOR- 
TRAN IV for the computer at the Lewis Research Center. This computer is an IBM 
709411/7044 direct couple computer using IBSYS version 13. 

QUASC is applicable for the type of seals shown in figures 1 and 2 and for narrow 
slots. The program must be supplied with the geometry of the seal, gas properties, 
reservoir conditions, constants for determining the variation of mean friction factor 
with Reynolds number, and certain logical variables which control output and plotting. 
Input and output can be in either the International System of Units (SI) or the U. S. cus- 
tomary system of units. 

Automatic plotting of some of the output variables is an available option. The ver- 
sion of QUASC in use at the Lewis Research Center uses the DD80 microfilm plotter 
(ref. 3). Users of other computing systems would have to adapt that part of the program 
to fit the equipment available. 

In general, QUASC performs the following operations in analyzing the flow across a 
seal: It reads the input data and checks that these data are consistent. For instance, 
there are three input variables which can determine the flow length. Since only two are 
necessary, the third must be made consistent with the other two. When the input data 
have been read, QUASC determines the minimum film thickness for choked flow h*. 
Then, for each given film thickness h, QUASC determines which type of flow is asso- 
ciated with this h, that is, which solution of the equations is required as described in 
the section Solution of Equations (subsonic, critical, or supercritical). 

When the type of flow has been determined, QUASC calculates the following param- 
eters: entrance Mach number; mean friction factor; Reynolds ntimber; distributions, 
across the seal face, _of pressure, temperature, density, velocity, Mach number, and 
friction parameter 4f(L - x)/D; mass and volume flow rates; Knudsen number; sealing 
dam force; center of pressure; and where appropriate, rotational Reynolds number, 
variables associated with power dissipation, and axial film stiffness. 

When these data have been calculated for all the given film thicknesses, QUASC 
writes the output data with appropriate labels and headings. 
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To increase program efficiency and to facilitate program writing, QUASC includes a 
number of subprograms. Figure 8 shows the hierarchy of the subprogram calls. Vari- 
ables are transmitted between the main program and the subprograms through labeled 
COMMON storage. 

A more detailed description of QUASC and descriptions of the subprograms are 
found in appendix D. 

The formulas for many of the parameters calculated by QUASC are listed in table I. 
Not listed are the formulas for h*, Mj, f, and Re, which are discussed in detail in the 
section Solution of Equations. It should be noted that when the flow is in the transitional 
or turbulent flow regime, the power loss due to rotation is not calculated. Also, there 
is no Reynolds number criterion for determining turbulence due to rotation. A complete 
list of the variable names used by the program is found in the program listing in appen- 
dix E. More details on the numerical integration by Simpson’s rule, which is used to 
find the sealing dam force and center of pressure, and on the numerical differentiation,' 
used to find film stiffness, are found in appendix D. 


REMARKS ON ADDITIONAL PARAMETERS CALCULATED BY COMPUTER PROGRAM 

As previously stated the program will calculate rotational Reynolds number, vari- 
ables associated with power dissipation, and axial film stiffness. These additional pa- 
rameters can be useful in selecting the seal design and seal operatii^ design gap. 

The rotational Reynolds number is determined by usii^ the gap as the characteristic 
length, since gap film seals generally operate with merged rotational boundary layers 

The Knudsen number can be found from (ref. 4) 


_ Molecular mean free path 2. 96 M^^^ 
Mean film thickness 


Pp 


Under conditions of very small film thicknesses, the IQiudsen number may be greater 
than 0. 01, and this continuum analysis would no longer be valid. (In that case, a slip 
flow analysis must be used. ) 

The total power is found by considering the viscous shear due to rotation only 
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Power = RJ2 


X Shear force = ^ ^ 
h 



h 


where Rfi is the mean rotational velocity. 

A rough estimate of the gas temperature rise through the film can be found by as- 
suming that all the energy dissipated by viscous shear is added to the gas film: 

^ _ /iR^n^A 

film,av 

This calculated film temperature rise will be higher than actually occurs since the pre- 
dominant mode of heat transfer, conduction by the walls, is neglected. 


INPUT DATA 

Input to QUASC is by punched cards. The NAMELIOT feature of FORTRAN IV is 
used to read data pertaining to the seal, the gas properties, and the reservoir condi- 
tions. The feature allows the individual variables to be named and eliminates compli- 
cated card formats. Also, it permits changing just a few variables on each READ state- 
ment, which keeps data cards at a minimum in parametric studies of candidate seal 
designs. 

The first card required by QUASC is a title card. The title identifies the data and 
uses columns 1 to 72 of one card. It is read by format (12A6). 

The next cards required by QUASC contain the parameters in NAMELICT/SDATA/. 
Variables in /SDATA/ are listed in table II. When the data from NAMELIST/SDATA/ 
have been read, QUASC checks the value of SKPH. If SKPH is . FALSE. , cards with 
film thickness data are needed. The first card contains the variable NJ, which is the 
number of film thicknesses. It must be in the first three columns of the card in for- 
mat (13). Then there must be NJ film thicknesses punched six to a card in format 
(6F12.6). 

The next cards required contain the variables from NAMELIRT/PDATA/ . The 
variables in /PDATA/ are listed in table HI. 


OUTPUT 

Computer output consists of the input data and calculated parameters. If input is in 
SI units, output is also in SI units. If input is in U. S. units, output is also in U. S. units. 
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The printed output parameters are identified and the units of each are also printed. A 
sample of the output data appears in appendix F. 

Page one of the output contains the title which identifies the data; the input data as 
they are read from cards; the checked input data; the calculated parameters seal face 
area and gas constant; a list of what optional parameters will be calculated, namely 
power, normalized force and center of pressure, and distributions across the seal face; 
and a list of parameters that will be plotted. The legend at the bottom of the page iden- 
tifies the flow regime associated with a particular film thickness. The legend is as 
follows: 

indicates the critical film thickness 
+ indicates the flow is in the laminar -turbulent transitional regime 

/ indicates the flow is in the turbulent regime 

Page two contains lists of parameters that vary only with film thickness. These pa- 
rameters are mass flow rate, standard volume flow rate, Knudsen number, mean free 
path of the gas molecules, axial film stiffness, sealing dam force, center of pressure, 
normalized sealii^ dam force, normalized center of pressure, friction factor, pressure 
flow Reynolds number, rotational flow Reynolds number, choking pressure, choking dis- 
tance, power, heat generation due to viscous shearing, apparent temperature rise due 
to power dissipation, and torque. Choking pressure and choking length are not printed 
for supercritical flow since they do not apply to that calculation. Power and parameters 
based on power dissipation are not printed for the transitional or turbulent flow regions 
since they are not calculated. 

Page three and the following pages contain lists of parameters that vary with radial 
distance, namely, pressure, Mach number, friction parameter 4fL/D, density, veloc- 
ity, and temperature. The film thickness to which each set of lists pertains is printed 
above the lists. 


SAMPLE PROBLEM 


An example of the use of the computer program is given here with the following con- 
ditions: Air at 4. 48x10^ N/m^ (65. 0 psia) is to be sealed from ambient air at 1. 03x10^ 

2 

N/m (15. 0 psia) by a seal operating in the externally pressurized mode. The mean 
temperature of the pressurized air is 310. 9 K (100° F). The sealing dam outside diam- 
eter is 0. 1684 meter (6. 630 in. ), and the inside diameter is 0. 1658 meter (6. 530 in. ). 
The design surface speed is 61 meters per second (200 ft/sec). 

It is desired to find a design mean film thickness which is large enough so that 
power dissipation and viscous heating temperature rise are sufficiently low, yet small 
enough so that the mass leakage is tolerable. From our experience, the best method 
is to try mean film thickness inputs of 2. 54 to 25. 4 micrometers (0. 1 to 1. 0 mils) in 
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increments of 2. 54 micrometers (0. 1 mil). However, to give a sample output for tran- 
sition flow and turbulent flow, the range of film thicknesses has been increased to 45. 72 
micrometers (1. 8 mils). The output desired includes the mass and volume flow rates, 
sealing dam force due to the pressure drop, center of pressure, power loss, and ap- 
proximate temperature rise due to shearing. A check on the validity of the analysis is 
made by examining the Knudsen number and the rotational flow Reynolds numbers. For 
this study, isentropic entrance conditions are assumed. Thus, the prc^ram input will 
include 


Mean rotational velocity, V, m/ sec (ft/ sec) 

Molecular weight of gas, m, kg/kg-mole (Ibm/Ib-mole) . . . 
Reservoir pressure (highest pressure), Pq, N/m (psia) . . 
Ambient pressure (lowest pressure), Pg, N/m^ (psia) . . . . 
Specific heat at constant pressure, Cp, J/kg-K (Btu/lbm- R) 
Film thickness (increase in increments of 2. 54 pm 

(0. 1 mil)), h, pm (mil) 

Reservoir temperature, Tq, K (°F) 

Inner radius, Rj, m (in. ) 

Flow length, m (in. ) 

Specific heat ratio, 

h 



kt 

"t 

Loss coefficient, Cj^ 


61 (200) 

.... 28.9660 (28.9660) 

. . . . 4.48x10"^ (65.0) 

1.03x10^(15.0) 

10^ (0. 24) 

2. 54 to 45. 72 (0. 1 to 1. 8) 

310 (100) 

0.0829 (3.265) 

.... 1. 29x10"^ (0. 050) 

1. 4 

24 

1.00 

0.079 

0.25 

1.00 


The data sheet for this sample problem is given in table IV. The sample problem is 
shown in both SI and U. S. units. The choking film thickness is 13. 104 micrometers 
(0. 516 mil); film thicknesses of 27. 94, 30. 48, and 33. 02 micrometers (1. 1, 1. 2, and 
1. 3 mils) are in the transition flow regime; and film thicknesses of 35. 56 micrometers 
(1. 4 mils) and greater are in the turbulent flow regime. The distributions of pressure, 
Mach number, friction parameter, density, radial velocity, and temperature are also 
shown. Plots (f^. 9) are shown for film thickness ratios and also for film thicknesses 
of 2.54, 13.1, and 45. 72 micrometers (0. 100, 0.516, and 1. 80 mils). This problem 
ran in approximately 0. 5 minute on the Lewis computer. 


CONCLUDING REMARKS 


A summary of the quasi-one-dimensional analysis of compressible flow with friction 
has been presented. Also, a detailed description of the computer program QUASC 
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which performs this analysis is given. This computer program has proved to be very 
useful in seal design. An example mainshaft seal problem required about 0. 5 minute to 
analyze 18 different film thickness conditions. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 29, 1972, 

132-15. 
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APPENDIX A 


SYMBOLS 


A area 

B friction parameter, 4f(L - x)/D 

entrance velocity loss coefficient 

c speed of sound 

c specific heat at constant pressure 

P 

D hydraulic diameter (D = 2h for coaxial rii^s and narrow slots) 

F sealii^ dam force 

F dimensionless sealing dam force 

f mean Fanning friction factor 

h film thickness (gap) 

k numerical constant in friction factor - Reynolds number relation 

Kn Elnudsen number 

L flow length from entrance to point of chokii^ 

M Mach number, u/'^y^T 

M mass flow 

m molecular weight of gas 

n exponent in friction factor - Reynolds number relation 

P pressure 

wetted flow perimeter 
Q volume leakage flow rate 

R mean radius 

AR sealing dam radial width (physical flow length), R 2 - Rj 
dt gas constant (universal gas constant/molecular weight) 

Re Reynolds number 

Re(p) Re 3 molds number of leakage flow (defined in table I) 

Re(r) Reynolds number of rotational flow (defined in table I) 



r radius 

T temperature 

u leakage flow mean velocity (x-direction) 

V mean rotational velocity 

W flow width 

X radial coordinate direction 

x^ center of pressure 

Xp dimensionless center of pressure 

y specific heat ratio 

X mean free path of gas molecules 

fji absolute (dynamic) viscosity 

p density 

T shear stress 

Subscripts: 

av average 

I laminar flow 

max maximum 

rot rotational 

t turbulent flow 

w wetted surface area 

x location along flow leakage length from entrance 

0 sealed (reservoir) conditions 

1 entrance conditions 

2 exit conditions 

3 ambient sump conditions 
Superscript: 

* critical flow condition 
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APPENDIX B 


DERIVATIONS AND REMARKS 
Transition Flow Friction Factor Derivation 


For laminar flow 


f = k^(Re)"^ (Bl) 

For turbulent flow 

f = k^(Re) ^ (B2) 

Equations (Bl) and (B2) are straight lines when plotted on logarithmic graph paper. 

Taking logs (to any arbitrary base) we have log f = log k^ - n^ log Re and 
log f = Ic^ kj - n^. log Re. Define x = log Re and y = log f . Therefore, for laminar 
flow, y = log k^ ■ X and for turbulent flow, y = log k^. - n^x. 

Let Xj be the value of x such that flow is laminar for x < Xj. Let Xg be the 
value of X such that flow is turbulent for x > Xg. Note that in general Xj # Xg. To 
smooth the discontinuity between Xj^ and X 2 , define a new function y 

3 2 

y = ax +bx +cx + d 


subject to the constraints 


ylx2 = ^ - V2 


dx „ dx 
^1 


*2 


= 0 


These constraints result in four linear equations in a, b, c, and d which can be solved 
simultaneously. These equations are 


^ 2 

axj + bxj + cxj + d = log k^ - 
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3 2 

ax 2 + bXg + CX 2 + d = log kj, - 1 x^X 2 
3ax^ + 2bXj + c = 0 
3ax2 + 2 bx 2 + c = 0 

The determinant of these equations is nonzero; and the solution for a, b, c, and d 
gives y as 

2[iog(k./kp +(n^Xj^ ~ x^(xj^ + xg) + 3 xjX 2X - 1/2 x^CSxj - X2)l 

y = log kj - n^.X 2 

(X2 - xj)3 


Comment on Solution of Mach Number Equation for Critical Flow 

It was mentioned in the main text that equation (23) does not always have a solution. 
To illustrate, if equation (23) is rewritten as 


g(Mj ) = 




-1 


y/(y-l) 


2 M. 


,-|l/2 




(B3) 


and plotted for PQ/P 3 = 3.0, y = 1. 4, and Cj^ = 0. 8 , it is obvious that g(Mj) does not 
cross the Mj-axis and, hence, does not have a solution in the range of acceptable Mj’s 
(see fig. 10 ). 
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APPENDIX C 


DETAILS OF NUMERICAL TECHNIQUES USED IN PROGRAM 
Method of Iteration 


This method is a standard numerical analysis technique that is used for finding the 
roots of an equation. It consists of writing an equation 


g(x) = 0 


(Cl) 


in the form 


X = G(x) 


(C2) 


Any X which satisfies equation (C2) is a root of equation (Cl). 

Equation (9) in the main text is an equation whose roots must be found numerically, 
so it can be used to illustrate the method of iteration. Equation (9) can be written as 


g(M)=lJLz_M! + 21±iin 






-{y + 1)M^ 
2 


1 + 1 (y - 1)M^ 


B = 0 


(C3) 


where B is a constant 


g _ 4f(L - x) 

2h 

Rewritii^ equation (C3) so it is in the form of equation (C2) gives 



(C4) 


(C5) 
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The first step in finding an M which will satisfy equation (C5) is to guess a value for 
M. Then it must be substituted in G(M), which will give a new value for M. If the two 
values of M agree, the equation is satisfied. If they do not, the new value of M is 
substituted in G(h^ to give yet another value for M. This procedure is repeated until 
M = G(M). Table V gives the steps in the solution of equation (C3) by the method of 
iteration for y = 1. 4, B = 1, and an initial guess for M of 0. 1. By examining a plot of 
g(M) against M, as shown in figure 11, it is readily observed that g(M) crosses the 
M-axis at M = 0. 5087. Hence, M = 0. 5087 is a root of g(M). 

The program QUASC uses the method of iteration to solve equations such as equa- 
tion (C3). 


Numerical Differentiation for Film Stiffness 

The differentiation required for the calculation of axial film stiffness, -dF/dh, is 
performed numerically by Lagrange's method (ref. 5). The general formula for La- 
grange differentiation is 


n 


L'(x) = ^ L|(x)y. 
i=0 


where 


and 



P.(x) = (X - x.)-l 


j=0 


For the special case of the derivative at one of the tabulated points Xj^, 
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where D.j^ = (xj^ - xpp.(xp, i k. 

Reference 5 describes a computing scheme which is followed in QUASC. A matrix 
is defined whose elements are a-- = x. - x.. The product of the off-diagonal elements 

IJ I J +1-| 

of each row is multiplied by Xj^ - Xj^, except for the k row. This defines D.j^. The 
sum 


n 



i=0 

i^tk 


is formed and multiplied by the product of the negative of the off-diagonal elements in 
the column. This gives the first term in the formula. The second term is formed 
by summing the reciprocals of the off-diagonal elements in the k*"^ column and multi- 
plying the sum by yj^. 


Numerical Integration by Simpson's Rule 

Simpson’s rule is a common method of numerical integration. The formula is 
written 



y dx = y(xQ> + 4y^Xg + + y(xj)j 


(C7) 


where Ax = x^ - Xq. The total interval of integration can be broken up into smaller sub- 
intervals. This is written as 




y dx 


23 



The total Interval must be broken up into enough subintervals to ensure that the value of 

f^i 

J y dx is accurate. The subroutine which does the numerical integration by Simp- 
le 

son's rule makes smaller subintervals where the value of the integrand is changing 


rapidly. 
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APPENDIX D 


DETAILED DESCRIPTION OF PROGRAM 

This section contains details regarding checking the consistency of the input data, 
traps for invalid data, and the subprograms. The variables mentioned in this appen- 
dix are defined and listed either in the section INPUT DATA or in the program 
listings. 


Main Program QUASC 

The main program QUASC controls the complete quasi-one-dimensional analysis. 

It defines the labeled COMMON storage for transmission of data among the subprograms. 
It defines some constants and labels for the output. It reads the input data, checks that 
these data are consistent, and rejects cases for which the input data are invalid. It 
calls subprograms to determine h*; M^; f; Re; and the distributions across the seal 
face of pressure, temperature, density, velocity, Mach number, and friction parameter 
4f(L - x)/D. It calculates the other parameters associated with the flow, such as mass 
flow rate, and writes the output. It then calls a subroutine to do the automatic plotting. 
QUASC then transfers to read new input data. 

The labeled COMMON storage defined by QUASC contains constants needed by the 
subprograms and the output data. The variables in the COMMON blocks are listed in 
tables VI, Vn, and Vni. The variable names and array names are not listed because 
the names are not the same in all the programs. 

There are two kinds of parameters in the COMMON block /ARRAYS/ (table VII): 
one is the distributions across the face of the seal, the other includes some of the pa- 
rameters that vary with film thickness. The arrays are dimensioned for 20 film thick- 
nesses and 11 points across the seal face. The two exceptions are the distance array 
and the pressure array P. These are dimensioned for 101 points across the seal face. 
The finer mesh is necessary for the numerical integration of the pressure function in 
calculating force and center of pressure (see table Vn). 

The parameters in common block /CONVRT/ are the constants needed by the pro- 
gram for the internal calculations. The array CONV has dimension (12, 2). The first 
column contains constants for calculations in SI units. The second column contains con- 
stants for calculations in U. S. customary units. The variable lUNITS is 1 for SI units 
and 2 for U. S. units. 

QUASC first reads the input data and checks that these data are consistent. One 
parameter it checks for consistency is the flow length, RDIF. There are three input 
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variables that can be used to determine RDIF. These variables are RlIN, R2IN, and 
RDIFIN. Since only two are necessary, the third must be made consistent with the other 
two. The check is made as follows: If RlIN * 0 and R2IN * 0, QUASC sets Rl = RllN 
and R2 = R2IN and then calculates RDIF = R2 - Rl. If RlIN ^ 0, R2IN = 0, and 
RDIFIN * 0, QUASC sets Rl = RlIN and RDIF = RDIFIN and then calculates 
R2 = Rl + RDIF. If RlIN = 0, R2IN * 0, and RDIFIN * 0, QUASC sets R2 = R2IN and 
RDIF = RDIFIN and then calculates Rl = R2 - RDIF. Any other combination of RlIN, 
R2IN, and RDIFIN is considered an error since there is not enough information avail- 
able to determine a nonzero flow length. If this error condition exists, QUASC writes 
a message and transfers to read new data from NAMELIST/SDATA/ . 

Another parameter that is checked for consistency is the pressure ratio. Three 
input variables are available for determining the pressure ratio, but only two are neces- 
sary. The check is made the same way as the check for flow length. If POIN * 0 and 
P3IN # 0, QUASC sets PO = POIN and P3 = P3IN and then calculates PRAT = P0/P3. If 
POIN * 0, P3IN = 0, and PRIN # 0, QUASC sets PO = POIN and PRAT = PRIN and then 
calculates P3 = PO/PRAT. If POIN = 0, P3IN 0, and PRIN * 0, QUASC sets 
P3 = P3IN and PRAT = PRIN and then calculates PO = PRAT*P3. Any other combination 
of POIN, P3IN, and PRIN is considered an error. In that case, QUASC writes an error 
message and transfers to read new data according to the input code. 

Another parameter that must be checked for consistency is the rotational velocity. 

If SPEED # 0 and CAPV = 0, CAPV is calculated from SPEED. If CAPV * 0 and 
SPEED = 0, SPEED is calculated from CAPV. If both SPEED and CAPV = 0, the sys- 
tem is considered static and the logical variable PWRSKP is set to . TRUE . to omit cal- 
culations involving power. 

The final check for consistency is made on the lexical variables. If any of the logi- 
cal variables are set to . TRUE. , which means a calculation is to be omitted, the corre- 
sponding element in the PLTSKP array must also be set to . TRUE. 

QUASC also tests for pressure ratios which do not yield a solution to equation (23). 
In that case, the flow is assumed to be subcritical and h* is set to 10^®. 

For valid pressure ratios, QUASC determines the critical film thickness h by 
calling subroutine START. For other film thicknesses, QUASC compares h with h . 

If h > h*, subroutine CHOKE is called to determine M^, f, and Re for supercritical 
flow. If h < h*, subroutine NCHOKE is called to determine Mj, f, and Re for sub- 
critical flow. And subroutine DIST is called to calculate the distributions across the 
face of the seal. QUASC then calculates the other flow parameters. Function subpro- 
gram SIMPS2 is used for the numerical integrations in the calculation of force and cen- 
ter of pressure. When all the data for all the film thicknesses have been calculated, 
QUASC calls subroutine STFNSS to determine the axial film stiffness. 
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When all the calculations are complete, QUASC writes the input data, the "checked" 
input data, and the output data with appropriate headings and labels. Then QUASC calls 
subroutine GRAFIC to do the automatic plotting. 

The final command in the prc^ram is a transfer to read new input data. 


Subroutine START 

Subroutine START determines the critical film thickness h* and Mj, i, and Re(p) 
for h*. It uses three arithmetic function statements which define Pj, Tj, and 
4f(L - x)/D (eqs. (19), (18), and (9)). The calculation scheme is described in the main 
text section Solution of Equations. There are some intermediate variables defined in 
the iteration procedure which are not mentioned in the solution of the equations. Their 
purpose is mainly to delineate the individual steps in the iterative solution. 


Subroutine CHOKE 

Subroutine CHOKE determines Mj, f, and Re(p) for supercritical flow. It uses 
three arithmetic function statements which define P^, T^, and - x)/t> (eqs. (19), 
(18), and (9)). The computing scheme is described in the section Solution of Equations. 
The iteration stops when L = AR. 

In the iteration to determine Mj, CHOKE saves three values of each L and Mj. 

It then does a linear curve fit (least squares) on the three points (L, M^) and picks from 
the line a new value of Mj that corresponds to L = AR. This is used to generate 
a new L, which is checked against AR. If they are not equal, the point such that 
I L - AR I is largest is discarded and the new point is substituted for it. Then a new 
linear curve fit is made and the process is repeated \mtil | L - AR| < 10”"^ 

If the curve fit produces a new Mj which is less than zero, Mj is chosen to be 
one-half the smaUest of the three existing M^^'s. 

As its final step, CHOKE calculates T*, which will be needed in calculating the 
temperature distribution across the face of the seal. 


Subroutine NCHOKE 

Subroutine NCHOKE determines Mj, f, and Re for subcritical flow. It uses five 
arithmetic function statements to define Pj, Tj, 4f(L - x)/D, P/P*, and Bg (eqs. (19), 
(18), (9), (14), and (32)). The computation scheme is described in the section Solution 
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of Equations. The fferation stops when Pg = P 3 . 

In its iteration, NCHOKE saves three values of each Mj and Pg. It does a linear 
(least squares) curve fit on these three points (Pg, Mj) and picks from the line a new 
value of Mj that corresponds to Pg = P 3 - Then a new Pg is generated and checked 
against P 3 . If they are not equal, the point such that | Pg - Pg | is largest is discarded 
and the new point is substituted for it. Then a new Pg is calculated and the process 

is repeated until 1^2 ” ^ 3 ! ^ ^tol- c 

A cut-off criterion of P^.^, = 5P,xlO ® is necessary because of roundoff errors on 

the computer. For the special cases of h* = 10^^ and for h very close to h (approx- 
imately 1 10'^ h* |), a three-point curve fit is not necessary. Two points were found to 
be sufficient for the straight-line curve fit since the Mach numbers involved are very 
low. 

The values of Mj used are screened systematically by the subroutine. If 
generates an L < AR, is lowered by AMj until an Mj^ is found which will produce 
an L > AR. 

The final calculation in NCHOKE determines T*, which is needed in calculating the 
temperature distribution across the face of the seal. 


Subroutine DIST 

Subroutine DIST determines the distributions of pressure, temperature, density, 
Mach number, velocity, and friction parameter 4f(L - x)/D across the face of the seal. 
The subroutine uses two arithmetic function statements to define Tj and Bg (eqs. (18) 
and (32)). The computation scheme is described in the section Solution of Equations. 

Pressure is calculated at 101 points since a fine mesh is needed for the numerical 
integration in the force and center-of-pressure equations. The other parameters are 
calculated at 11 points across the face of the seal. 


Subroutine STFNSS 

Subroutine STFNSS performs the numerical differentiation of force with respect to 
film thickness, which gives the axial film stiffness. 

There are four regions in which the differentiation is performed. These four re- 
gions are for subcritical laminar flow, subcritical turbulent flow, supercritical laminar 
flow, and supercritical turbulent flow. The four regions are necessary because the 
force varies differently with film thickness for the different types of flow. The axial 
film stiffness in the transition region is not computed because the exact values of the 
sealing dam force are unknown. 
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The subroutine first arranges the data in order of ascending x. It chooses the five 
data points x., y. (i = k - 2, k - 1, k, k + 1, k + 2) for use in equation (C6). If there 
are less than five data points in the set, all are used in the differentiation. 

The computation scheme for STFNSS is described in detail in appendix B. 


Subroutine SIMPS2 

Function subprogram SIMPS2 uses Simpson's rule to evaluate the integrals in the 
force and center -of-pressure equations. A statement such as 

F = SIMPS2(J, XO, XF, G, K) 

gives F as the definite integral 


F = 



The integrand is evaluated at interior points by the fvinction subprogram G names in the 
calling vector. The parameter J is a constant needed by G. In the case of QUASC, J is 
the index of the film thickness being examined. 

The range of the integration is not divided uniformly. More subdivisions are made 
in the region where the integrand is changing rapidly. If SIMPS2 cannot evaluate the in- 
tegral to an accuracy of 10 ^ times the value of the int^rand with 200 subdivisions, the 
parameter K is raised by 2 to indicate the integral is inaccurate. 


Subroutines FFUNC and XCFUNC 

Function subprc^ram FFUNC evaluates the integrand in the force equation. Function 
subprogram XCFUNC evaluates the integrand in the center-of -pressure equation. Both 
FFUNC and XCFUNC call subroutine PRESS to evaluate the pressure for the j film 
thickness at any distance x across the seal face. 


Subroutine PRESS 

Function subprc^ram PRESS uses a three-point Lagrange interpolation to determine 
the pressure at any distance x across the face of the seal. One-himdred and one values 
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of X and P are stored for this interpolation. Subprogram PRESS first searches the x 
array to find the three stored values of x closest to the given x and then uses these 
three points in the interpolation. 


Subroutine GRAFIC 

Subroutine GRAFIC makes the followii^ plots: 

(1) Power against film thickness 

(2) Center of pressure against film thickness 

(3) Force against film thickness 

(4) Pressure against x 

(5) Temperature against x 

(6) Density against x 

(7) Mach number against j£ 

(8) Friction parameter 4f(L - x)/D against x 

All the plots are dimensionless. The film thickness is normalized by the maximum h. 
Radial coordinate direction x is normalized by AR. Power is normalized by the maxi- 
mum power. Center of pressure is normalized by AR. Force is normalized by 
(Pjj - Pj)ARffR. Pressure is normalized by Pq. Temperature is normalized by Tq. 
Density is normalized by Pq = Pq/RTq. Mach number and friction parameter are al- 
ready dimensionless. 

The subroutine presented here is just a skeleton. It has no commands that will pro- 
duce plots. Each user can add the plotting commands appropriate to the equipment 
available. The subroutine does, however, Include commands for the normalization of 
the plotting parameters and labels for the plots as they appear in this report. 


Subroutine SORT 

Subroutine SORT rearranges the ordered pairs of numbers (X(l), Y(l)), 

(X(2), Y(2)), . . . , (X(N), Y(N)), which are stored in arrays X and Y in order of ascend- 
ing X. The (X, Y) pairing is preserved. 
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APPENDIX E 


PROGRAM LISTING AND FLOW CHARTS 


PROGRAM FOR OUASl-ONE DIMENSIONAL COMPRESSIBLE FLOW WITH FRICTION 
AND PARALLEL FILMS 


INPUT VARIABLES 


TITLE - ALPHA-NUMERIC IDENTIFICATION OF THE DATA 

RUN - INNER RADIUS OF SEAL 

R2IN - OUTER RADIUS OF SEAL 
ROIFIN - FLOW LENGTH 
WIDTH - FLOW WIDTH 
MOLWT - MOLECULAR WEIGHT 
CP - SPECIFIC HEAT 
MJ - RESERVOIR VISCOSITY 
GAMMA - RATIO OF SPECIFIC HEATS 
SPEED - ROTATIONAL SPEED 

CAPV - SURFACE VELOCITY 

XLAM - EXPONENT IN RE - F RELATION FOR LAMINAR FLOW 

CONLAM - CONSTANT IN RE - F RELATION FOR LAMINAR FLOW 

XTURB - EXPONENT IN RE - F RELATION FOR TURBULENT FLOW 

CONTRB - CONSTANT IN RE - F RELATION FOR TURBULENT FLOW 

RELAM - UPPER LIMIT OF REYNOLDS NUMBER FOR LAMINAR FLOW 
RETURB - LOWER LIMIT OF REYNOLDS NUMBER FOR TURBULENT FLOW 
PWRSKP - LOGICAL VARIABLE, IF TRUEt SKIP POWER CALCULATIONS 
NRMSKP - LOGICAL VARIABLE. IF TRUE, SKIP NORMALIZATION OF FORCE 
AND CENTER OF PRESSURE 

PRSSKP - LOGICAL VARIABLE. IF TRUE, SKIP PRINTOUT OF DISTRIBUTIONS 
PLTSKP - LOGICAL ARRAY, IF TRUE, SKIP THE APPROPRIATE PLOT 
PLTSKP(l) - POWER VS FILM THICKNESS 
PLTSKP(Z) - CENTER OF PRESSURE VS FILM THICKNESS 
PLTSKPO) - FORCE VS FILM THICKNESS 
PLTSKPI4) - PRESSURE VS DISTANCE ACROSS SEAL 
PLTSKP(5) - TEMPERATURE VS DISTANCE ACROSS SEAL 
PLTSKP(6) - DENSITY VS DISTANCE ACROSS SEAL 
PLTSKP(T) - MACH NUMBER VS DISTANCE ACROSS SEAL 
PLTSKPIB) - 4FL/D VS DISTANCE ACROSS SEAL 

SKPH - LOGICAL VARIABLE, IF TRUE, SKIP READING OF FILM THICKNESS 
DATA 

NOSI - LOGICAL VARIABLE, IF TRUE, INPUT AND OUTPUT ARE IN U.S, 
CUSTOMARY UNITS. IF FALSE, INPUT AND OUTPUT ARE IN 
INTERNATIONAL UNITS. 

NJ - NUMBER OF FILM THICKNESS IN ADDITION TO H* 

H - ARRAY OF FILM THICKNESS 

POIN - SUPPLY PRESSURE 

P3IN - AMBIENT PRESSURE 

PRIN - SUPPLY TO AMBIENT PRESSURE RATIO 
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JOIN - TOTAL TEMPERATURE 
LOSS - VELOCITY LOSS COEFFICIENT 
INCODE - INPUT CODE 

INCODE = 1 MEANS TRANSFER TO READ NEW TITLE CARD 
INCODE * 2 MEANS TRANSFER TO READ NEW FILM THICKNESS 
DATA 

INCODE = 3 MEANS TRANSFER TO READ NEW RESERVOIR DATA 
OUTPUT VARIABLES 


Rl - INNER RADIUS OF SEAL 
R2 - OUTER RADIUS OF SEAL 
RDIF - FLOW LENGTHf R2-R1 
P3 - SUPPLY PRESSURE 
P3 - AMBIENT PRESSURE 
PRAT - PRESSURE RATIO, PO/P3 
TO - TOTAL TEMPERATURE 
AREA - SURFACE AREA OF SEAL 
RGAS - GAS CONSTANT 

MOOT - MASS FLOW RATE 
Q - VOLUME FLOW RATE 
KN - KNUDSEN NUMBER 

LAMDA - MEAN FREE PATH OF GAS MOLECULES 
STIFF - AXIAL FILM STIFFNESS 
FORCE - SEALING DAM FORCE 
XC - CENTER OF PRESSURE 
XCBAR - NORMALIZED CENTER OF PRESSURE 
FBAT - NORMALIZED FORCE 
FRICT - MEAN FRICTION FACTOR 
REP - PRESSURE REYNOLDS NUMBER 
RE^ - ROTATIONAL REYNOLDS NUMBER 
PSTA^ - PRESSURE AT POINT OF CHOKING 
LSTAR - DISTANCE FROM ENTRANCE TO POINT OF CHOKING 
POWER - POWER DUE TO VISCOUS SHEARING 
HSHEAR - SHEAR HEAT 

DELT - APPARENT TEMPERATURE RISE DUE TO POWER DISSIPATION 
TORQUE - TORQUE DUE TO POWER 

X - DISTANCE ACROSS SEAL 
P - PRESSURE DISTRIBUTION ACROSS SEAL 
MACH - MACH NUMBER DISTRIBUTION ACROSS SEAL 

FF - FRICTION PARAMETER 4FL/D DISTRIBUTION ACROSS SEAL 
RHD - DENSITY DISTRIBUTION ACROSS SEAL 
V - VELOCITY DISTRIBUTION ACROSS SEAL 
T - TEMPERATURE DISTRIBUTION ACROSS SEAL 
CALC - ARRAY OF LABELS 

PROGRAM VARIABLES 


C - PERMANENTLY STORED ARRAY OF LABELS 
PI -- 3.1415927 

RUNIV - UNIVERSAL GAS CONSTANT 


TEMP - TEMPORARY STORAGE FOR A CONSTANT 

_ MACH NUMBER FOR WHICH P3/P0 IS A MAXIMUM 
CC - TEMPORARY STORAGE IN COMPUTING CRITICAL 


PRESSURE RATIO 
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PRCRIT - PRESSURE RATIO THAT CORRESPONDS TO CRITH * MAXINUN 
PRESSURE RATIO FOR WHICH EQ. 23 HAS A SOLUTION 

M - CHARACTERISTIC OF LOG(IO) OF P3 
PTOL - TOLERANCE FOR ROUND-OFF ERROR IN ITERATIVE SOLUTION OF 
EQUATIONS FOR SUB-CRITICAL FLOW 
DELX - STEP SIZE ACROSS FACE OF SEAL 

MUM - VISCOSITY AT MID-SEAL 

K - INDICATOR THAT NUMERICAL INTEGRATION IN FORCE CALCULATION 
IS INACCURATE 

K< - INDICATOR THAT NUMERICAL INTEGRATION IN CENTER OF 
PRESSURE CALCULATION IS INACCURATE 
J - FILM THICKNESS INDEX 

H(J=1» IS THE CRITICAL FILM THICKNESS H* 

REAL MOLWT,MDOT,KN,LAMDA, HUM, MACHrLSTAR, LOSS 
LOGICAL PWRSKP,NRMSKP,PRSSKP,PLTSKP(8),SKPH,N0SI 
REAL MU 

COHMON/ARRAYS/X( ion ,P(20, 101 ) , MACHI 20, ll ) , VI20, 11 ) ,T{ 20, 11 ) , 

1 RHOI 20,11 ) ,FF( 20,1 i) ,H(20) ,XCBAR(20) ,LSTAR(20) ,PSTAR(20) , 

2 REP(20»,FRICT( 20) , POWER (20), FORCE! 201 .STIFF (20) ,FBAR(20) 
COMMON/CONSTS/GAMMA,RDIF,XLAM,CONLAM,XTURB,CONTRB,TO,PO,P3,PTOL, 

1 RGAS, LOSS, RELAM, RETURB, MU, TSTAR 

COMMON/CONVRT/CONV(12,2) .lUNITS 

NAMELIST/SDATA/Rl IN, R2IN,RD IF IN, WIDTH,MOLWT,CP,MU, GAMMA, SPEED, 

1 CAPV,XLAM,CONL AM, XTURB.CONTRB, RELAM, RETURB, PWRSKP.NRMSKP, 

2 PRSSKP,PLTSKP,SKPH,NOSI 
NAMELIST/PDATA/POIN,P3IN,PRIN,TOIN,LOSS,INCODE 

DIMENSION TITLE(12),C1(4),C2(4).C3(4),C4(4),C5(4),C6(4) ,C7(4) 
DIMENSION CALC(7,4),RUNIV(2),MOOT(20).0(20I,KN(20),LAMDA(20), 

1 RER(20),DELT(20),HSHEAR(20) ,T0RQUE(20) ,XC(20) 

DATA (Cim,I = l,4)/6H ,6H P0W.6HER ,6H / 

DATA (C2(I),I=«1,4)/6HDIMENS,6HI0NLES,6HS QUAN,6HTITIES/ 

DATA (C3(I),I=1,4)/6H PRESS, 6HURE 01 ,6HSTRIBU,6HTI0NS / 

DATA IC4(I), 1*1, 4)/6HP0WER ,6H ,6H XC(,6HBAR) / 

DATA (C5(I),I*1,4)/6HF0RCE ,6H ,6H PRES.6HSURE / 

DATA (C6m,I = l,4)/6HTEMPER,6HATURE ,6H 0EN,6HSITY / 

DATA (C7(I),I*1,4)/6HMACH N,6HUMBER ,6H 4,6HFL/0 / 

DATA BLANK/6H / 

DATA PI,RUNIV(1),RUNIV(2)/ 3.1415927, 8.31436E3, 1545.4 / 

EXTERNAL FFUNC,XCFUNC 
MU = D. 

READ INPUT DATA - SEAL IDENTIFICATION 

SEAL DIMENSIONS, PLOTTING COOES, PHYSICAL 
PROPERTIES OF GAS 
SEVERAL GAPS BETWEEN PLATES 

100 READ (5,1) TITLE 
READ (5,SDATA) 
lUNITS*! 

IF(NOSI) IUNITS=2 

CHECK CONSISTANCY OF SEAL DIMENSIONS AND LOGICAL VARIABLES 

IFIRIIN.EQ.O.) GO TO 101 
IF(R2IN.EQ.O.) GO TO 102 
R1 = RUN 
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R2 * R2IN 
RDIF* R2-R1 
GO TO 104 

101 IF(R2IN*EQ.O*) GO TO 103 

IF (RDIFIN.EQ.O. ) GO TO 103 
RDIF= ROIFIN 
R2 ^ R2IN 
R1 = R2-RDIF 
GO TO 104 

102 IF (RDIFIN.EQ.O. » GO TO 103 
ROIF= ROIFIN 

R1 = RUN 
R2 = Rl+RDIF 
GO TO 104 

103 WRITE (6y32) R1 I N» R2 IN, RDI FIN 
IF CS<PH) GO TO 106 

READ (5,2) NJ 
NJ = NJ>1 

IF (NJ.GT.l) READ (5,3) (H ( J) , J:=2, NJ ) 

106 READ ( 5,P0ATA) 

GO TO 700 

104 IF(NRMSKP) PLTSKP(3)= .TRUE. 

IF(.NOT.PRSSKP) GO TO 110 

DO 105 I»4,8 

105 PLTSKP(I)* .TRUE. 

110 IF (SKPH) GO TO 120 
READ (5,2) NJ 
NJ = 

IF (NJ.GT.l) READ (5,3) ( H ( J ) , J>2, N J ) 

120 READ (5,PDATA) 

H(l)* l.E-3 

IF (NJ.LT.2) GO TO 130 
DO 125 J»2,NJ 

125 H(l) * AMINKHdltHU)) 

DEBUG H(l) 

WRITE INPUT DATA 

130 WRITE (6,10) TITLE 

IF (NOS I ) WRITE (6, 11) R1IN,P0IN,M0LWT,R2IN,P3IN,CP,C0NLAH,RDIFIN, 

1 PRIN,GAHHA,XLAM, WIDTH, TOIN, MU, RELAH, LOSS, SPEED, CONTRByCAPV, 

2 XTURB,RETURB 

IF(.NOT.NOS1 ) WRITE(6,13) R1 IN,POIN, M0LWT,R2IN,P3INtCP,C0NLAM, 

1 RDI F IN, PRIN, GAMMA, XL AM, WIDTH,TOIN, HU, RELAH, LOSS, SPEED, 

2 CONTRB ,C AP V , XTURB , RETURB 

CHECK CONSISTANCY OF PRESSURE DATA 

200 IF (P0IN.NE.0«) GO TO 240 

IF (P3IN.NE.0.) GO TO 210 

WRITE (6,30) 

GO TO (100,110,120), INCOOE 
210 IF (PRIN.NE.O^) GO TO 230 
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c 


220 

WRl 

[TE 

(6,31) P0IN,P3 

IN 


GO 

TO 

(100,110,120), 

INCODE 

230 

PO 

s 

PRIN*P3IN 



P3 

= 

P3IN 



PRAT » PRIN 
GO TO 260 

240 IF {P3IN,NE.0.I GO TO 250 

IF (PRIN.EO.O.) GO TO 220 

P3 = POIN/PRIN 
PO = POIN 
PRAT * PRIN 
GO TO 260 

250 PO = POIN 
P3 = P3IN 
PRAT » PO/P3 

CALCULATE CRITICAL PRESSURE RATIO 

H* CANNOT BE DETERMINED FOR P3/P0( GIVEN) GREATER THAN P3/P0 
(CRITICAL) - I.E., EQ. 23 DOES NOT HAVE A SOLUTION 
ASSUME NON-CHOKED FLOW FOR THIS CONDITION 

260 IF (GAMMA. EQ. 1. ) GO TO 300 

TEMP » (GAMMA+1. )/2. /(GAMMA-1.) 

CRITM = SQRT(LOSS**2-TEMP*SQRT(TEMP**2+LOSS**2*(LOSS**2-1.))) 

CC = CRITM*SQRT( l.+. 5*( GAMMA-1. )*CRITM**2)/( L0SS»*2+.5*(GAMMA-1. ) 
1 *CRITM**2)**(GAMMA/(GAMMA-1.)) 

PRCRIT = CC*SQRT(2. /(GAMMA+1. ))*L0SS»*(2.*GAMMA/(GAMMA-1.)) 
IF(P3/P0.LT.PRCRIT) GO TO 300 
WRITE (6,56) P0.P3, GAMMA, LOSS 

H(l) » 1.E30 

WRITE PROGRAM CONSTANTS AND RESERVOIR CONDITIONS 

300 RGAS » RUNIVI IUNITS)/MOLWT 
TO = TOIN+CONVd, lUNITS) 

FOR AIR, CALCULATE RESERVOIR VISCOSITY FROM SUTHERLANDS FORMULA 

IF(GAMMA.EQ.1.4) MU-CONV ( 2 , lUN ITS )*T0**1. 5/ ( TO+CONV ( 3 , lUNI TS ) ) 

TO =• T0-C0NV(1,IUNITS) 

N AL0G10(P3) 

PTOL » 5.*10.«*(N-6) 

IF (SPEED. EQ.O.) GO TO 330 

CAPV» PI*SPEED*(R1+R2 )/C0NV(4, lUNITS) 

GO TO 340 

330 IF (CAPV.EQ.O. ) GO TO 335 

SPEED* C0NV(4,IUNITS)*CAPV/PI/(R2+R1) 

GO TO 340 

335 PWRSKP - .TRUE. 

CP * 3 • 

340 IF (PWRSKP) PLTSKP(1)*.TRUE. 

AREA * PI*(R2**2-R1**2) 

IF (WIDTH. EQ. 0. ) WIDTH=PI*(R1+R2) 

DELX s ROIF/100. 

DO 350 1-1,101 
350 X(I) - FLOAT! 1-1 )*DELX 

IF (PWRSKP) PLTSKP(1)=.TRUE. 
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DO i60 1=1,7 

DO 363 J=l,4 

360 CALCCI , J)= BLANK 

DO 383 1=1,7 

GO TO (361,362,363,364,365,366,367), I 

361 IF(PW^SKP) GO TO 380 

DO 1361 J=i,4 

1361 CALCn , J)= Cl( J) 

GO TO 380 

362 IF(NR'ISKP) GO TO 380 

DO 1362 J=l,4 

1362 CALC(I,J)= C2(J) 

GO TO 380 

363 IF(PRSSKP) GO TO 380 

DO 1363 J=i,4 

1363 CALC(I,J)= C3(J) 

GO TO 380 

364 IFCPLTSKP(l) ) GO TO 1364 
CALC(I,1)= C4(l) 

CALC(I,2)= C4(2) 

1364 IF(PLTSKP(2) ) GO TO 380 
CALCd ,3)= C4(3) 

CALC(I,4)= C4(4) 

GO TO 380 

365 IF(PLTSKP(3) ) GO TO 1365 
CALCd, 1)= C5(l) 

CALCd, 2)= C5(2) 

1365 IF(PLTSKP(4) ) GO TO 380 
CALCd, 3)= C5(3) 

CALCd, 4)= C5(4) 

GO TO 380 

366 IF(PLTSKP(5) ) GO TO 1366 
CALCd, 1>= C6(l) 

CALCd ,2)= C6(2) 

1366 IF(PLTSKP(6n GO TO 380 
CALCd, 3)= C6(3) 

CALCd, 4)= C6C4) 

GO TO 380 

367 IF(PLTSKP(7) ) GO TO 1367 
CALCd ,1)= C7(l) 

CALCCI,2)= C7(2) 

1367 IFCPLTSKP(8)) GO TO 380 
CALCd ,3)= C7(3) 

CALCd, 4)= C7(4) 

380 CONTINUE 
C 

IF(NOSI) WRITE(6,12) R1 , PO, MOLWT, ( CALC ( 1 , I) , 1=1 ,4) , R2,P3,CP, 

1 (CALC(2t I) ,I=1,4),RDIF,PRAT,GAMMA, {CALCC3, I) ,1=1,4) , WIDTH, 

2 TO, MU, (CALC (4, I ), I =1,4) , AREA, RGAS, SPEED, (CALC(5,I ) ,1=1,4), 

3 L0SS,CAPV,(CALC(6,I),I=1,4),(CALC(7,I),I=1,4) ^ 

IF( .N3T.N0SI ) WRITE( 6,14) R 1 , PO, MOLWT, ( CALC ( 1 , I ) , I =1 ,4 ) , R2 , P3 ,CP , 

1 (CALC(2,I) ,1=1,4 ),RDIF, PRAT, GAMMA, (CALC (3, I), I =1,4), WIDTH, 

2 T0,MU,(CALC(4,I),I=1,4),AREA,RGAS,SPEED, (CALC( 5 , I ) . I =1 , 4 ) , 

3 LOSS.CAPV, (CALC(6, n,I=l,4),(CALC(7,I),I=l,4) 

C 

TO = TO + C0NV(1,IUNITS) 

C 

C H(l) IS H*, THE CRITICAL FILM THICKNESS. CALCULATE IT FIRST. 
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C IF H(J) IS GREATfR THAN H*. ASSUME CHOKED FLOW. IF H(J) IS LESS 

C than H*, ASSUME NON-CHOKFD FLOW, 

C 

400 DO 473 J=1,NJ 

IF (J.NE.l) GO TO 410 

IF (HI 1).GE.1,E30) GO TO 470 

CALL START 

IF(H(l) ,LT, 1.E30) GO TO 415 
GO TO 470 

410 IF (HI J).GT.H(l) ) GO TO 420 
412 CALL NCHOKE(J) 

GO TO 415 
420 CALL CHOKEIJ) 

415 CALL DIST(J) 

C 

440 MDOT(J) = WIDTH*H(J) *RHO( J, 1 ) * V ( J , U *CONV ( 5 , 1 UNI TS ) 

Q(J>= M00T(J)*C0NV(6, lUNITS) 

KNIJ) = 2.96*MACH(J, 11)/REP( J) 

LAMDAI J) = KN( J)*H(J > 

MUM = MU 

IF ( GAMMA. EQ. 1.4) MUM=C0NV ( 2 , I UN I TS ) * ( T ( J , 6 ) +CONV ( 1 , I UNI TS ) )»*1,5 
1 /(T(J,6)+C0NV( 1,IUNITS)+C0NV(3,IUNITS) ) 

RER(J) = RHOI J,6)*CAPV*H( J)/MUH*2./C0NV(7,IUNITS) 

IF (PWRSKP) GO TO 460 

C 

C CALCULATE POWER, SHEAR HEAT, TORQUE, AND TEMPERATURE RISE DUE TO 

C POWER DISSIPATION 

C 

IF( REP (J).LE, RELAM) POWER ( J ) =MUM*AREA*CAPV**2/H ( J )*CONV ( 8 , 1 UNITS ) 
DELTU )= C0NV(9, lUNI T S )*POWER( J )/ ABS (HOOT ( J ) I /CP 
HSHEARI J)= C0NV(9,IUNITS)*P0WER( J) 

TORQUE! J)= POWER! J)*CONV( 10, lUNITS )/SPEEO 
C 

C CALCULATE FORCE AND CENTER OF PRESSURE 

C 

460 K=0 

KK = D 

FQRCE(J) = SIMPS2I J,0.,RDIF,FFUNC,K)*WIDTH 
XC(J) = SIMPS2(J,0.,RDIF,XCFUNC,KK)/F0RCE( J)*WIDTH 
IFIK.NE.O.OR.KK.NE.O) WRITE(6,63) J,K,KK 
470 CONTINUE 

CALL STFNSS(NJ) 

C 

C WRITE OUTPUT DATA 

C 

SCO write(6,io) title 

WRITE! 6,17) 

IF (.NOT.NRMSKP) WRITE(6,51) 

IF (NDSI) WRITE (6,18) 

IF (.NOT.NOSn WRITE (6,15) 

IF (.NOT.NRMSKP) WRITE(6,52) 

DO 513 J=1,NJ 

IF !H( J).GF.1.E30) GO TO 510 

WRITE (6, 19) H(J) ,MDOT (J ) ,Q ( J ) ,KN! J ) ,LAMDA( J) ,STIFF ( J) , 

1 FORCE(J), XC!J) 

IF(J.EQ.l) WRITE(6,58) 

IF (REP(J).GE.RETURB) WRITE (6,60) 

IF(REP! J).GT.RELAM.AND,REP( J ).LT.RETUR6) WRITE (6, 59) 

IF (NRMSKP) GO TO 510 
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XCBARIJJ* XC(J)/RDIF 

FBAR(J)= FOkCE( J)/«P0-P3»/WIDTH/RDIF 
WRITE 16,20) FBAR( J) ,XCBAR( J) 

510 CONTIMJE 
C 

WRITE(6,21) 

IF {.NOT.PWRSKP) WRITE(6,53) 

IF (NOSI) WRITE (6,22) 

IF (.'JDT.NOSD WRITE (6,16) 

IF (,'^OT.PWRSKP.AND.NOSI ) WRITE (6,54) 

IF (.MOT, PWRSKP.ANO. .NOT.NOSI) WRITE (6,33) 

DO 523 J=1,NJ 

IF (H( J).GE,i.E30) GO TO 520 
WRITE(6,19) H(J), FRICT( J) ,REP(J),RER(J) 

IF (J.EQ.l) WRITE (6,58) 

IF(REP ( J) .GT.RFLAM,AND.REP( J).LT,RETURB) WRITE (6, 59) 

IF (REP( J).GE,RETURB) WRITE (6,60) 

IF(H(J).LE.H(D) WRITE(6,57) PSTAR(J), LSTAR(J) 

IF {, MOT. PWRSKP. AND. REP(J).LE, RELAM) WRITE (6,24) POWER(J), 

I HSHEAR(J) ,DELT(J),TORQUE(J) 

520 CONTINUE 
C 

IF (PRSSKP) GO TO 600 
DO 533 J=1,NJ 

IF (M3D( J,3).E0.1) WRITE (6,26) 

IF (H(J).GE,1.E30) GO TO 533 

IF (NOSI) WRITE (6,25) H(J) 

IF (.MOT. NOSI) WRITE (6,34) H(J) 

IF (J.EQ.l) WRITE (6,58) 

531 WRITE(6,27) 

IF (NOSI) WRITE (6,35) 

IF (.MOT. NOSI) WRITE (6,36) 

DO 532 1=1,101,10 
IJ = I/lO+l 

WRITE (6,23) X( I ),P( J, I) ,MACH( J, IJ) ,FF( J,IJ) ,RHO( J,IJ) , 

1 V(J,IJ),T( J,IJ ) 

532 CONTINUE 

533 CONTINUE 

HAKE PLOTS OF - POWER, CENTER OF PRESSURE, AND FORCE VS FILM 
THICKNESS AND OF PRESSURE, TEMPERATURE, DENSITY, MACH NUMBER, AND 
FRICTION PARAMETER 4FL/D VS RADIUS 

600 CALL GRAFIC(PLTSKP,NJ, TITLE) 

700 GO TO (100,110,120), INCODE 

END OF PROGRAM 

1 FORMAT (12A6) 

2 FORMAT (13) 

3 FORMAT (6F12.0) 

10 FORMAT (1H1,52HQUASI-0NE DIMENSIONAL COMPRESSIBLE FLOW SEAL PROGRA 
1M,5X,12A6) 

11 FORMAT (IHO, 12HINPUT DATA -, /,1H0, 10X,9HR1, INCHES, 22X,7HP0,PSIA, 

I 18X,16HM0LECULAR WEIGHT, 17X, 9HF=K/RE**N, /, IH ,F18. 4,F30. 3, 

"^2 F29.3,22X, 9H*» •****•*, /, IHO, 10X,9HR2, INCHES, 22 X ,7HP3, PS I A, 

3 18X,16HCP,BTU/LBM-DEG R, 17X , 10HK( LAMINAR) ,/, IH ,F18.4, F30. 3 , 

4 F28.3,F30.3,/, 1H0,6X,18HFL0W LENGTH, INCHES, 18X,5HP0/P3,24X, 

5 5HGAMMA,23X,10HN(LAMINAR),/,1H ,F18.4,F29.3,F29.3,F30.2/, 
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6 IH0,6X,17HFL0W WIDTH, INCHES, 17X.8HT0, DEG F,16X, 

7 20HVISCaSITY,LB-SeC/IN2,8X,24HUPPFR LIMIT RE (LAMINAR),/, 

8 IH ,Fi8.4,F29. 1,G3.’,4, F28.1,/, 1H0,39X,10HL0SS COEF, ,21X, 

9 9HSPEtO,RPM,19X,12HK(TURBULENT),/,lH ,F47, 2 , E32. 4,G30. 4 , / , 

X 1H0,70X,8HV, FT/SEC, 20X,12HN(TURBULENT),/,1H , F 77.2 , G32. 4/ , 

1 lH0,9iX,26HL0WER LIMIT RE ( TURBULENT ),/, IH ,F107.1) 

12 FORMAT ( IHC, 13HOUTPUT DATA - , / , IHO, 1 OX, 9HR1 , INCHE S , 22X ,7HPO , PS I A , 
18X,:bHMOLECULAR W El GHT, 1 7X, 9HCALCULATE , / , IH ,F18.4,F30.3, 
F29.3,22X,9H*****#***,/, IHO, 10X,9HR2,INCHES,22X,7HP3,PSIA, 
18X,16HCP,BTU/LBM-DEG R, lOX , 4A6, / , IH , F18, 4, F30. 3 , F28. 3 ,/ , 
IH ,92X,4A6,/, IH ,6X,18HFLOW LENGTH, INCHES, 18X, 5HP0/P3,24X, 
5HGAMMA,/, IH .F18.4,F29.3,F29, 3,16X,4A6,/1H0,6X, 

17HFLOW WIDTH, INCHES, 17X,8HT0, DEG F.16X, 

20HVI SCOSITY, LB-SEC/ IN2,18X,4HPL0T,/,1H ,F18.4,F29.3,G32.4, 
14X, 22H****»** »»»**»**»»#♦*»*,/, IHO, 10X,8HAREA, I N2,14X, 
24HGAS CONSTANT, LB-FT/LBM-R,14X,9HSPEED,RPM,14X,4A6,/,1H , 
F18.4,F29.5,G32.4,/, IH ,93X,4A6,/, IH ,39X,10HLOSS COEF., 
21X,8HV,FT/SEC,/, IH , F47. 2 , F30. 2, 16X,4A6 , /, 1H0,93X,4A6 ,/ , 
IHO, 44X,30H*** *»»*»»#*♦»#***»************, /,1H ,44X,1H*, 
28X,1H*,/,IH ,44X,30H* * - CHOKING FILM THICKNESS *,/,lH , 
44X,30H* + - TRANSITION REGION *,/,lH ,44X, 

30H* / - TURBULENT FLOW *,/,lH ,44X, 1H*,28X, 

1H*,/,1H ,44X, 30H« «•»»*«»•**•»*•••«••»»«••»•**») 

13 FORMAT (1H0,12HINPUT DATA -,/, IHO, 10X,9HR1 , METERS, 22X,7HPO, NSMA, 

1 18X.16HM0LECULAR WE I GHT , 17X, 9HF=K/RE**N, /, IH ,E 18. 3 , E30. 3 , 

2 E29. 3, 22X,9H** ****•**,/, IHO, 10X,9HR2, METERS, 22X,7HP3, NSMA, 

3 17X,18HCP, JOULES/KG-DEG K , 16X, lOHK (LAM INAR ) , / , IH ,E18.3, 

4 E30.3,E28.3,F30.3,/,1H0,6X,18HFL0W LENGTH, METERS, 10X, 

5 5HP0/P3,24X,5HGAMMA,23X, lOHN(LAMINAR) , /,1H , E18. 3, E29. 3 , 

6 F29.3,E30.2,/, 1H0,6X,17HFL0W WIDTH, METERS, 17X,8HTO, DEG K, 

7 17X,18HVISC0SITY,N-SEC/M2,9X,24HUPPER LIMIT RE (LAMINAR),/, 

8 IH ,E18.3,F29.1,G32.4,F28.1,/,1H0,39X,10HL0SS COEF.,21X, 

9 9HSPEED,RPS,19X,12HK(TURBULENT) ,/,lH ,F47.2,E32.4,G30.4,/, 

X 1H0,70X,7HV,M/SEC,21X,12HN(TURBULENT),/,1H ,F77.2,G32.4,/, 

1 1H0,91X,26HL0WER LIMIT RE ( TURBULENT ),/, IH ,F107.1) 

14 FORMAT (1H0,13H0UTPUT DATA -,/, IHO, 10X,9HR1, METERS, 22X,7HPO, NSMA, 

1 18X,16HM0LECULAR WEIGHT, 17X,9HCALCULATE,/,1H , El 8. 3 , E30. 3 , 

2 F29. 3, 22X,9H** *****•*,/, IHO, 1 OX, 9HR2, METERS, 22X,7HP3, NSMA, 

3 17X,18HCP, JOULES/KG-DEG K,9X,4A6,/,1H , E18. 3 , E30. 3 ,E28. 3 , / , 

4 IH ,92X,4A6,/, IH ,6X,18HFLOW LENGTH, METERS, 18X,5HPO/P3,24X, 

5 5HGAMMA,/,1H , E18. 3, F29.3,F29. 3,16X,4A6,/, 1H0,6X, 

6 17HFLOW WIDTH, METERS, 17X,8HTO. DEG K,17X, 

7 18HVISCOSITY,N-SEC/M2, 19X,4HPL0T,/,1H ,E18.3,F29.3,G32.4, 

8 1-VX,22H********»*«*********#»,/,1H0,11X,7HAREA,M2,12X, 

9 28HGAS CONSTANT, JOULES/KG-DEG K, 12X,9HSPEED,RPS,14X,4A6, /, 

X IH ,E18.3,E29.3,G32.4,/,1H , 93X.4A6, /, IH ,39X,10HL0SS COEF., 

1 21X,7HV,M/SEC,/, IH , F47. 2 , E30. 3, 16X, 4A6, / , IHO, 93X, 4A6 , / , 

2 1H0,44X,30H«»**«»«**««****»**»»»**»««»*»«,/.ih ,44X,1H«,28X, 

3 1H*,/1H ,44X,30H* * - CHOKING FILM THICKNESS *,/,lH ,44X, 

4 30H* + - TRANSITION REGION *,/,lH ,44X, 

^ 30H* / - TURBULENT FLOW *,/,lH ,44X, 1H*,28X,1H*,/, 

6 IH .44X,30H«***««»««***»«***»»*»»**»«**»») 

15 FORMAT (IH , 5X,6HMETERS , 7X , 6HKG/SEC , 7X, 4HSCMS, 8X ,6HNUMBER,6X , 

1 6HMETERS,7X,3HN/M, lOX, 1HN,9X, 6HMETERS) 

16 FORMAT (IH , 5X , 6HMET ERS , 8X , 6HF ACTOR , 30X, 4HNSMA , 7X,6HMETERS ) 

17' FORMAT (1H0,15H FILM THICKNESS , 3X, 6HM( DOT) ,9X, IHO, 9X, 7HKNUDSEN, 

1 5X,6HLAMBDA,4X,9HSTIFFNESS,5X,5HF0RCE,9X,2HXC) 

18 FORMAT (IH , 5X ,6H INC HES , 7X , 6HLB/M I N, 7X, 4HSCFM , 8X ,6HNUMBER, 6X , 

1 6HINCHES,6X,5HLB/IN,8X,2HLB,9X,6HINCHES) 
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19 

20 
21 

22 

23 

24 

25 

26 
27 

1 

30 

31 

32 

33 

34 

35 

36 

51 

52 

53 

54 
56 


FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 
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60 
63 


FORMAT (IH ,G15.5,7G12.3) 

FORMAT (1H+, 102X,2G12.3) 

FORMAT (1H0,15H FILM TH I CKNE SS , 3X , 8HFRI CTI ON f 4X , 5HRE I P ) , 7X f 
5HRE(R) ,8X,4HP(*>,8X,4HL(*n 

FORMAT (IH ,5X,6HINCHES,RX,6HFACTOR,30X,4HPSIA,7X,6HINCHES) 

(IH .7G12.3) 

{1H+,77X,4G12.3) 

(1H0,18H FILM THICKNESS =, G12. 5 , 2X , 6HINCHES ) 

(IHl) 

( IHOt 6X, IHX.lOX 1 1HP» lOXt 4HMACH. 8X, 5H4FL/D 1 7X » 7HDENSI TY »7X , 
1HV,11X,1HT) 

FORMAT IlHO.llHPO = P3 = 0) 

FORMAT (1H0,4HPU = , F IG. 3 , 5X, 4HP3 = , F 10. 3 t 5X, 9HP0/P3 = 0) 

FORMAT ( lH0t39HSEAL DIMENSION DATA INCONSISTANT - R1 =fF10.4,5X, 
4HR2 =,F10.4,5X,6HRDIF =,F10.4) 

FORMAT ( 1H+ f 8lXf 5HWATTS. 6X . 5HWATTS. 7X» 5HDEG Kt9Xt3HN-M) 

FORMAT (IHOflSH FILM THICKNESS = , G12. 5 t 2X , 6HMETER S) 

FORMAT (IH ,4X,6HINCHES,5X,4HPSIA,8X,6HNUMBER,17X,11HLB-SEC2/FT4, 
2X,6HFT/SEC,7X,5H0EG F) 

FORMAT (IH , 4X t 6HMET ERS» 5X t 4HNSMAf 8X» 6HNUMBER f 20X» 5HKG/M3 * 6Xf 
5HM/S£C,6X,5H0EG K) 

FORMAT (IH+t 105X,5HF0RCEf8Xf2HXC) 

(1H+, 106X,3HBAR,9Xt3HBAR) 

(1H+,81X,5HP0WER,4X,9HHEAT L0SS,4X,6HDEL(T) t7X,6HT0R0UE ) 
(lH+t82X,4HH. P.f 5X, 7HBTU/MINt 5X,5H0EG F 1 9X . 5HLB/F T ) 

( IHOf 23HCANN0T FIND H* FOR PO =fG14.5,6H, P3 *tG14.5f 
9H, GAMMA =fG14.5,18Hi AND LOSS COEF. *fG14.5»/tlH t 
22HASSUME NON-CHOKFD FLOW) 

FORMAT (1H+,51X,2G12.3) 

(1H+|1H*I 
(1H+.2H +) 

(^!_j4* 2H /) 

( IHolsSHERROR IN NUMERICAL INTEGRATION J =»I5»5Xt3HK =t 
I5t5X,4HKK =tI5) 


FORMAT 

FORMAT 

FORMAT 

FORMAT 


FORMAT 

FORMAT 

FORMAT 

FORMAT 


END 


tIBFTC STRT 


FIND CRITICAL FILM THICKNESS AND ENTRANCE CONDITIONS FOR CRITICAL 
FLOW 

PROGRAM VARIABLES 


Ml - ENTRANCE MACH NUMBER 

Cl - CONSTANT IN SEVERAL EQUATIONS - TEMPORARY STORAGE 
MTEST - NEW VALUE OF Ml IN ITERATIVE SOLUTION OF EQ. 23 
Bl - FRICTION PARAMETER AT Ml (EQ. 9) 

TT - ENTRANCE TEMPERATURE 

MJ - ENTRANCE VISCOSITY 

VV - ENTRANCE VELOCITY 

PP3 - ENTRANCE PRESSURE RATIO, Pl/PO 
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C - ENTRANCE DENSITY 

C C2 - CONSTANT IN H* EQUATION 

C EK - EXPONENT IN FRICTION FACTOR - REYNOLDS NUMBER RELATION 

C FOR LAMINAR OR TURBULENT FLOW 

C CON - CONSTANT IN FRICTION FACTOR - REYNOLDS NUMBER RELATION 

C FOR LAMINAR OR TURBULENT FLOW 

C «****» 

C XI * 

C X2 * - INTERMEDIATE VARIABLES IN FRICTION FACTOR - REYNOLDS 

C XY * NUMBER RELATION FOR TRANSITION REGION 

C X.1 • 

C »****• 

C P - PRESSURE 

C REP - ENTRANCE PRESSURE REYNOLDS NUMBER 

C FRICT - MEAN FRICTION FACTOR 

C HNEW - TRIAL VALUE FOR H* 

C N - CHARACTERISTIC OF LOG(IO) OF HNEW 

C HTOL - ROUND-OFF ERROR TOLERANCE IN CALCULATION OF H* 

C 

SUBROUTINE START 

DOUBLE PRECISION X1,X2.XX,Y2 

REAL MU,MUO,MACH,Ml,MTEST,LSTAR,LOSS 

COMMON/ARRAYS/X( 101) ,P(2C, 10 1 ) , MACHI 20, 1 1 ) , V ( 20, 1 1 ) , T ( 20 , 11 ) , 

1 RH0(20,ll),FF(20,li) ,H(20),XCBAR(20) ,LSTAR(20) ,PSTARI20) , 

2 REP(20) , FRICT ( 20 ) , POWER ( 20 ) , FORCE ( 20 ) , ST IFF ( 20 ) ,FBAR(20) 
COMMON/CONSTS/G,RDIF,XLAM,CONLAM,XTURB,CONTRB,TO,PO,P3,PTOL, 

1 RGAS,LOSS,RELAM,RETURB,MUO,TSTAR 

COMMON/CONVRT/CONV(12,2) ,IUNITS 
EQ1(G,X) = 1./ ( 1.+.5* {G-1, )*X**2)**{G/(G-1, n 
EQ2(G,X) = l./(l.+,5*(G-l. )*X**2) 

EQ3(G,X) = ( l./X**2-l.+. 5* (G+1. )*ALOG( (G+1. )*X**2/(2.+IG-l. ) 

1 *X**2)))/G 

C 

IF (G.NE.l. ) GO TO 90 
Ml = P3/P0 
GO TO 110 


ITERATE FOR Ml BASED ON FORMULA USING M2=l, P2=P3, AND L*=R2-R1 
90 Ml = 0. 

Cl = P3*S0RT(.5*(G+1. ) )/P0 
100 MTEST = C1*SQRT( EQ2I G,M1 ) )/EQl(G,Ml/LOSS) 

IF (Ml.GT.l. ) GO TO 160 
IF (ABS(Ml-MTEST).LT.l.E-7) GO TO 110 
Ml = MTEST 
GO TO 100 
110 B1 = EQ3IG,M1) 

Cl = 2.*R0IF/B1 

CALCULATE TEMPERATURE, VISCOSITY, VELOCITY, PRESSURE, DENSITY, 
AND REYNOLDS NUMBER AT ENTRANCE. 

TT = EQ2(G,M1/L0SS)*T0 
MU = MUO 

IF(G.EQ.1,4» MU=CONV( 2, IUNITS)*TT**1.5/(TT+C0NV(3,IUNITS) ) 

VV= M1*SQRT (TT*RGAS*G*CONV( 11, I UN ITS I ) 

PPO = 1. 

IF (G.NE.l.) PP0=EQ1(G,M1/L0SS) 
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c 

c 

c 

c 


P(l,l) = PPO*PO 

RRH P(l,i)*C0NV(12t IUNITS)/RGAS/TT 
C2 = 2. *VV*RRHO/MU/CONV( 7, lUNITS) 

120 REP(l) = C2*H(1) 

FOR LAMINAR OR TURBULENT FLOWt CALCULATE FRICTION FACTOR = K/RE»*N 
FOR TRANSITION FLOW, USE SPECIAL FORMULA 


IFIREPI 1) .GT.KELAM) GO TO 130 
EX = XLAM 
CON = CONLAM 
GO TO 140 

130 IF(REP(1).GT.RETURB) GO TO 135 
XI = ALOG(RtLAM) 

X2 = ALOG(RtTURB) 

XX= ALOG(REPd)) 

Y2 = AL0GIC0NTRB)-XTURB*X2-2«*( AL0G(C0NTRB/C0NLAM)+XLAM*X1- 

1 XTURB»X2 ) * ( XX* (XX*(XX-1.5*( X1+X2 )»+3.*Xl*X2)-.5*X2**2 

2 *(3. *X1-X2) )/ (X2-X1 1**3 


FRICTll)= DEXP(Y2) 

GO TO 145 
135 EX= XTURB 

CON = CONTRB 

140 FRICT(l) = C0N/REP(1 )**EX 


C 

C ITERATE FOR H(1I=H* 

C 

145 HNEW= C1*FRICT(1) 

N = ALOGIO(HNEW) 

HTOL = 10,**(N-7) 

IF (ABS(H(l)-HNEW).LT.HTnL» GO TO 150 
H(i) = .5*(H(1)+HNEW) 

GO TO 120 


C SATISFACTORY H* FOUND - STORE L*, P*» MACH NUMBER, FRICTION 

C PARAMETER 4FL/D, DENSITY. TEMPERATURE, AND VELOCITY. 

C CALCULATE T*. 

C 

150 LSTAR(l) = RDIF 
PSTARI 1 ) = P3 
MACH(1,1) = Ml 
FE(1,1) = Bi 
RH0(1,1) = RRHO 
T(l,l)= TT-C0NVI1,IUNITS ) 

V(l,l) = VV 

TSTAR = TT/EQ2(G,Ml)/.5/(G+l.) 

RETURN 


ERROR CONDITION - SUPERSONIC ENTRANCE CONDITIONS. 

160 HID = 1.E3Q 

10 FORMAT ( IHO , 61HSUPERS0iNl C ENTRANCE CONDITIONS FOR H* , ASSUME NON-C 
IHOKED FLOW) 


RETURN 


C 


END 
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ilBFTC CSTRT 


STARTING ROUTINE FOR SUPERCRITICAL FLOW 
PROGRAM VARIABLES 


KNT - COUNTER FOR NUMBER OF TIMES THROUGH THE ITERATION 
Ml - ENTRANCE MACH NUMBER 
TT - ENTRANCE TEMPERATURE 
MU - ENTRANCE VISCOSITY 
PPO - PRESSURE RATIO, Pl/PO 
P - PRESSURE 
RHO - DENSITY 
V - VELOCITY 

REP - ENTRANCE PRESSURE REYNOLDS NUMBER 

EX - EXPONENT IN FRICTION FACTOR - REYNOLDS NUMBER RELATION 
FOR LAMINAR AND TURBULENT FLOW 

CON - CONSTANT IN FRICTION FACTOR - REYNOLDS NUMBER RELATION 
FOR LAMINAR AND TURBULENT FLOW 


XI * 

X2 * - INTERMEDIATE VARIABLES IN FRICTION FACTOR - REYNOLDS 
XX * NUMBER RELATION FOR TRANSITION REGION 

Y2 » 


FRICT - MEAN FRICTION FACTOR 

FF - FRICTION PARAMETER 4FL/D BY EO. 9 
FLl - FLOW LENGTH CALCULATED FROM EQ. 9 

SM - MACH NUMBER (DEPENDENT VARIABLE) FOR CURVE FIT 
SL - FLOW LENGTH (INDEPENDENT VARIABLE) FOR CURVE FIT 
FLMAX - MAXIMUM DIFFERENCE BETWEEN NEW FLOW LENGTH AND SL 
VALUES 

II - INDEX OF ELEMENT OF SL ARRAY FOR WHICH FLMAX IS LARGEST, 
REPLACE THE II (TH) POINT BY THE NEW POINT 
SUM3 - SUM OF SQUARES OF INDEPENDENT VARIABLE 
SUM2 - SUM OF INDEPENDENT VARIABLE 

SUMl - SUM OF PRODUCT OF INDEPENDENT AND DEPENDENT VARIABLES 
SUM3 - SUM OF INDEPENDENT VARIABLE 
DET - DETERMINANT IN CURVE FIT 

A - SLOPE OF FITTED STRAIGHT LINE 
B - INTERCEPT OF FITTED STRAIGHT LINE 
FM - VALUE OF MACH NUMBER FROM STRAIGHT LINE THAT 
CORRESPONDS TO L - R2-R1 
T - TEMPERATURE 

TSTAR - TEMPERATURE AT POINT OF CHOKING 


SUBROUTINE CHOKE (J) 

REAL MU,HUO,Ml,HTEST,MACH,LSTAR,LOSS 

DOUBLE PRECISION SL, SM, DET, SUM3, SUM2.SUM1 , SUMO, A,B,X1 ,X2.XX,Y2 
DIMENSION SM(3),SL(3) 

COMMON/ARRAYS/XdOl) , P (20, 101 ) , MACHI 20, 11 ) , V(20,ll ) ,TI20, 11 ) , 

1 RH0(20,11) ,FF(20,11),H(20),XCBAR(20),LSTAR(20),PSTAR(20), 

2 REPI20) ,FRICT( 2O),POWER(20),FORCE(2O),STIFF(2O),FBAR(2O) 
C0MM0N/C0NSTS/G,RDIF,XLAM,C0NLAM,XTURB,C0NTRB,T0,P0,P3,PT0L, 

1 RGAS, LOSS, RELAM, RETURB,MUO, TSTAR 



COMMOM/CONVRT/CONV(12f 2) tIUNITS 

EQKGfX) = I,/ ( 1,+,5*(G-1* »*X**2)**<G/(G-1* ) ) 

EQ2(G,X) = l./(l.+,5*(G-l. l*X**2l 

EQ3(Gt X) = ( l./X**2*l. + .5*<G+l. )«AL3G( .5*(G+l. )*X**2/ ( 1.+.5* (G-l. > 
1 *X**2)))/G 

C 

C GUESS HI 

C 

Ml = MACH (1,1) 

KNT = 0 
C 

C CALCULATE TEMPERATURE, VISCOSITY, PRESSURE, VELOCITY, AND 

C REYNOLDS NUMBER AT ENTRANCE 

C 

100 TT = TO*EQ2(G,M1/LOSS) 

MU = MUO 

IFCG.EQ.1.4) MU=CONV( 2, IUNITS)»TT**1.5/(TT+CONV(3, lUNITS) ) 

V(J,l) = Ml*SQRT(TT*RGAS*G*CONV(ll,IUNITS) ) 

PPO =* 1. 

IF (G.NE.l.) PP0=EQI(G,M1/L0SS) 

P(J,1) “ P0*PP0 

RHO( J, 1 ) » P (J , 1 ) *CONV ( 12 , lUN ITS ) /RGAS/TT 

REPU) * V( J,1)*RH0( J,1)*2.*H( J)/MU/C0NV(7, lUNITS) 

C 

C FOR LAMINAR OR TURBULENT FLOW, CALCULATE FRICTION FACTOR = K/RE**N 

C FOR TRANSITION FLOW, USE SPECIAL FORMULA 

C 

IF(REP(J).GT, RELAM) GO TO llO 
EX = XLAM 
CON = CONLAM 
GO TO 120 

110 IF(REP( J).GT.RETURB) GO TO 115 
XI = ALOG(RELAM) 

X2 = ALOG(RETURB) 

XX- ALOG(REP(J)) 

Y2 = ALOG(CONTRB)-XTURB*X2-2.*(ALOG(CONTRB/CONLAM)+XLAM*X1- 

1 XTURB*X2)*(XX*(XX*(XX-l.5*(Xl+X2) )+3.*Xl*X2)-.5*X2**2 

2 »(3, *Xl-X2) )/ (X2-Xl)«*3 
FRICT(J)* 0EXP(Y2) 

GO TO 125 
115 EX= XTJRB 

CON - CONTRB 

120 FRICTU)- CON/REP( J) »*EX 
C 

C CALCULATE FRICTION PARAMTER 4FL/0 FOR Ml 

C 

125 FF(J,l) - EQ3(G,M1) 

FLl - FF( J,i)»H( J)/2./FRICT( J) 

C 

C CHECK FOR CHOKING LENGTH = R2-R1. 

C IF NOT EQUAL, STORE Ml AND FLl FOR CURVE FIT 

C 

IF (ABS(FLl-RDIF).LT.l.E-7) GO TO 130 
KNT » KNT+1 

IF (KNT.GT,3) GO TO 126 
SMIKNT) » Ml 
SLCKNT) » FLl 
IF (KMT,EQ.3) GO TO 128 
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HI = M1»RDIF/ABS(FL1 ) 

GO TO 100 

126 FLMAX = 0. 

DO ;.27 1 = 1,3 

IF ( ABSISLt I)“FL1 ) .LT. FLHAX ) GD T3 127 
II = I 

FLMAX = ABS(SL( I )-FL 1 ) 

127 CONTIMUE 
SL(II) = FLl 
SMim = Ml 

128 CONTIMUE 

C 00 LIMEAR CURVE FIT WITH 3 POINTS TO DETERMINE NEW VALUE OF Ml 

C CORRESPONDING TO L=R2-R1 

C 

SUM3 = 0, 

SUM2 = 0. 

SUMl = 0, 

SUMO » 0. 

DO 129 1=1,3 

SUM3 = SUM3+SLI I )**2 

SUM2 » SUM2+SLI I ) 

SUHl = SUMl+SMI I )*SL ( I ) 

SUMO * SUMO+SM(I) 

129 CONTIMUE 

OET = 3.*SUM3-SUM2*»2 
A = (SUH1*3.-SUM2*SUM0)/DET 
B = (SUM3*SUM0-SUM2»SUM1 )/DET 
FH = A-RDIF+B 

IF IKMT.GT.3) Ml = . 5 * ( M 1 + FM ) 

IF IMl.LT.O.) M1=ABS(M:) 

GO TO 100 
C 

C SATISFACTORY Ml FOUND 

C 

130 T(J,1)= TT-CONVI 1 . lUN ITS ) 

TSTAR = TT/EQ2(G,Ml)/.5/(G+i. I 
MACHIJ.l) = Ml 

C 

RETURN 

END 


$IBFTC NCSTRT 
C 

C STARTING CONDITIONS FOR SUBCRITICAL FLOW 

C 

C PROGRAM VARIABLES 

C *******.»•******* 

c 

c LM - NUMBER OF POINTS IN THE CURVE FIT OF Ml VS P2 

C CUT - FACTOR BY WHICH Ml IS CUT IN SEARCHING FOR AN Ml 

C THAT GIVES L GREATER THAN R2-R1 

C Ml - ENTRANCE MACH NUMBER 


non 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


M2 - EXIT MACH NUMBER 

DELM - MACH NUMBER INCREMENT TO GET THE ITERATION STARTED 
KNT - COUNTER FOR THE NUMBER OF TIMES THROUGH THE ITERATION 

TT - E'iTRANCE TEMPERATURE 
MJ - ENTRANCE VISCOSITY 
VV - ENTRANCE VELOCITY 
PPD ~ PRESSURE RATIOf Pl/PO 
P - PRtSSURE 
RRHD - ENTRANCE DENSITY 
REP - ENTRANCE REYNOLDS NUMBER 

EX - EXPONENT IN FRICTION FACTOR - REYNOLDS NUMBER RELATION 
FOR LAMINAR AND TURBULENT FLOW 

CON - CUNSTANT in friction FACTOR - REYNOLDS NUMBER RELATION 
FOR LAMINAR AND TURBULENT FLOW 


XI * 

X; • - KirERMEDIATE VA«IABLES IN FRICTION FACTOR - REYNOLDS 
X< • NUMBER RELATION FOR TRANSITION REGION 

Y • 


FR - MEAN FRICTION FACTOR FOR TRANSITION REGION. IF FR 
IS NOT A GOOD NUMBERf ALTER Ml BY DELM. 

FRICT - MEAN FRICTION FACTOR 

R1 - FRICTION PARAMETER 4FL/0 FROM EO, 9 
LISIAR - FICTITIOUS FLOW LENGTH FROM ENTRANCE TO POINT OF 
CHOK ING 

PSTAR - pressure AT POIiT OF CHOKING 

L2 - DISTANCE FROM E vD OF SEAL TO POINT OF CHOKING 
B2 - FRICTION PARAMETER 4FL2/0 
MTEST - NEW VALUE OF M2 FROM ITERATIVE SOLUTION OF EQ. 9 
P2 - EXIT PRESSURE 

SM - MACH NUMBER (DEPENDENT VARIABLE) FOR CURVE FIT 
S? - EXIT PRESSURE (INDEPENDENT VARIABLE) FOR CURVE FIT 
PMAX - MAXIMUM DIFFEREfiCE BETWEEN NEW P2 AND STORED P2 VALUES 
II - INDEX OF ELEMENT OF SP ARRAY FOR WHICH PMAX IS LARGESTt 
REPLACE THE IIITH) POINT BY THE NEW POINT 
SUM3 - SUM OF SQUARES OF INDEPENDENT VARIABLE 
SUM2 - SUM OF INDEPENDENT VARIABLE 

SUMl - SUM OF PRODUCT OF INDEPENDENT AND DEPENDENT VARIABLES 
SUMO - SUM OF DEPENDENT VARIABLE 
DET - DETERMINANT IN CURVE FIT 

A - SLOPE OF FITTED STRAIGHT LINE 
B - INTERCEPT OF FITTED STRAIGHT LINE 

PDIF - MINIMUM DIFFERENCE BETWEEN P3 AND STORED VALUES OF P2 FOR 
THE CONDITION WHEN P2 CANNOT BE DETERMINED WITHIN PTOL 
TOLERANCE OF P3 

TSTA^ - TEMPERATURE AT CHOKING 


SUBROUTINE NCHOKEIJ) 

DOUBLE PRECISION SM, SP , A ,B ,DET, SUM3, SUM2 , SUMl , SUMO.Xl ,X2 , XX, Y2 
REAL MU,L1STAR,M1,M2,MTEST,L2.MACH,LSTAR,L0SS,MU0 
DIMENSION SP(3),SMI3) 

COMMON/ ARRAYS/XI 101) ,P( 20, 101) ,MACH( 20, 1 1 ) , V ( 20, 11 ) ,TI20,11 ) , 
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1 RHn(20,lH ,FF(20,in ,H( 20 ) ,XCBAR(20 J.LSTAR(20), PSTAR(20|, 

2 REP (20 fFRICT(20) , POWER (201 ,F0RCE(20 ) ,STIFF(20) t FBAR( 20) 
COMMON/CONSTS/G,RDIF , XLA M .CONLAM , XTURB ,C ONTRB ,TOf P0»P3, PTOLt 

I R5AS»L0SS,RELAM,RETURBtMU0fTSTAR 

DIMENSION SP(3)tSM(3) 

EOKGtX) = lo /(!.• ■•■o 5* (G-l» ) *X**2) **(G/ (G-l» ) ) 

E02(G,X) * l,c /(I j 5*(G-1. ) *X**2) . 

EQ3(GtX) = ( X, /X**2-l» +» 5* ( G+1* ) ♦ALOG( (G+l« ) *X**2/ ( 2« ♦( G-le ) 

1 ^ *X**2) ) ) /G 

E05(GfX) = SQRT( (G + l* ) / (2. + (G-lo ) *X**2) ) /X , ,*v 

EOB(G,B, X »= 1, /SQRT(G*B+1. -«5*(G+l. ) *ALOG( ( G+ 1« 1 *X**2/ ( 2. ♦( G-1* ) X 

1 **2)>) 


C 

C 

C 

C 

C 


C 

C 

C 

C 


SET CONSTANTS AND MAKE INITIAL GUESS FOR Ml, FOR P3/ PO GREAT 
THAN P3/P0 (CRITICAL)* USE A TWO POINT FIT, FOR P3/PC LESS T 
P3/P0 (CRITICAL) f USE A THREE POINT FIT 


L IM » 3 

IF (H(l ),GE,l,E30) LIM=2 

IF (ABS(1,0-H(1) /H( J) ), LE,2.0E-2) L1M=2 

CL'T » c9B 
MTEST * 0. 

Ml = ,75 
DELM = Ml/2, 

M2 = 1 
KNT = 0 

90 IF (MTEST.GT,, 9) CUT=, 9999 
IF (KNT.GT, 50) CUT=, 9999 
IF (KNToGT,100) CUT=, 9999999 

CALCULATE TEMPERATURE, VISCOSITY, VELOCITY, PRESSURE, DENSITY, 
AND REYNOLDS NUMBER AT ENTRANCE, 


TOO TT = E32CG,MWL0SS)*TC 
Ml = MUO 

IF (G.EO.1.4) MU=2.27E-8*TT**1«5/(TT+198,6) 

VV * Ml*SQRT(TT*RGAS*32-,l74*G) 

PPO = lo 

IF (G.NE.l.) PPO=EOl (G,M1 /LOSS) 

P(J,1 ) = PPO*PO 

PPHO = P ( J,l)*l44, /RGA S/TT/32. 174 
120 REP(J)* VV*RRH0*2.*H (J) /MU/12, 


FOR LAMINAR OR TURBULENT FLCW, CALCULATE FRICTION FACTOR = K/RE**N 
FOR TRANSITION FLOW, USE SPECIAL FORMULA 

IF(REP( J ),GT,RELAM) GO TO 130 
EX * XLAM 
CON = CONLAM 
GO TO 140 

130 IF(REP(J).GT.RETURB) GO TO 135 
XI = ALOG(RELAM) 

X2 = ALOG(RETURB) 
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XX= ALOG(REP(J») 

Y2 = AL0G(C0NTRB)-XTURB*X2-2.*CAL0GIC0NTR8/C0NLAM»+XLAM*Xi- 

1 XTURB*X2)*(XX»(XX*(XX-1.5*(X1+X2) »+3.*Xl*X2)-,5*X2**2 

2 *(3. *X1-X2)>/(X2-Xl)**3 
IF (K'iT.EQ.O) GO TO 131 

FR = DEXP(Y2) 

IF (FR.NE.FRICK J) ) GO TO 131 
Ml = Ml+OELM 
GO TO 100 

131 FRICT(JI= DEXP(Y2) 

GO TO 145 
135 EX= XTURB 

CON = CONTRB 

140 FRICT(J)= CON/REP(J)**EX 

CALCULATE FRICTION PARAMETER Bl, LI*. AND P* FOR Ml, DETERMINE 
L2*. FRICTION PARAMETER B2. M2, AND P2 

145 Bl ^ EQ3(G,M1) 

LISTAR = Bl*HU)/2./FRICT( JI 
PSTARUJ * P( J,1)/EQ5(G,M1) 

L2 = LISTAR-RDIF 
IF (L2.GE.0, ) GO TO 150 
MX = »11*CUT 
DELM * DELM*CUT 
GO TO 100 

150 B2 = 2,*FRICT(J»*L2/H( J) 

160 MTEST - EQ8{G,B2,M2) 

IF(ABS(H2-MTEST).LT, l.E-7) GO TO 170 
M2 = MTEST 
GO TO 160 

170 P2 = PSTAR(J)*EQ5(G,M2) 

CHECK FOR P2*P3. IF NOT EQUAL, STORE Ml AND P2 FOR CURVE FIT. 


IF 

IABS(P2-P3J,LT. 

PTOU 

GO TO 

IF 

(DELM.EQ.O.) 

GO 

TO 

180 

IF 

(DELM.LT. l.E-7} 


GO 

TO 175 

KNT 

» KNT+1 




IF 

(KNT.GT.LIM) 

GO 

TO 

171 

SMI 

KNT) » Ml 




SPI 

KNT) » P2 




IF 

(KNT.EO.LIM) 

GO 

TO 

173 

Ml 

= Ml+DELM 





GO TO 100 

171 PMAX » 0. 

DO 172 I«1,LIM 

IF (ABS(SP(I )-P3).LT.PMAX) GO TO 172 
II = I 

PMAX » ABS(SP(I)-P3I 

172 CONTINUE 
SPUD = P2 
SM(II) = Ml 

173 CONTINUE 

DO TWO OR THREE POINT LINEAR CURVE FIT TO DETERMINE NEW VALUE OF 
Ml CORRESPONDING TO P2=P3 
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c 

c 

c 


SUM3 = O.DO 
SUM2 * O.DO 
SUMl = O.D.; 

SUMO = 0,00 
DO 17'^ I = 1,LIM 
SUM3 = SUMi+SP( I )**2 
SUM? = SUM2+SP(I) 

SUMl = SUM1+SM( I ) *SP ( I ) 

SUM; * SUMO+SM(I) 

174 CONmUE 

OET = FL0AT(LIM»*SUM3-SUM2**2 

A = ( SUMl *FLOAT ( L IM ) -SUNi.'* SUPO ) /OET 

B = (SUM3*SUM0-SUM2*SUMi )/l)ET 

Ml = A*P3+B 

iELi"''ii;NllABS(«-SHl'nl»ltS(m-S«.21..ABS(Hl-SH(3n. DELHI 
GO TO 90 

ERROR CONDITION - CANNOT DETERMINE P2 WITHIN TOLERANCE PTOL OF P3 


175 


175 


GO TO 176 
M1=SM(3) 


C 

C 

C 

C 

C 


DELM = 0.0 
PDIF = ABS(SP(1)-P3) 

Ml = SM(1) 

IF (ABS(SP(2)-P3).GT.P0IF) 

PDIF = ABS(SP(2)-P3) 

Ml — S H ( 2 ) 

IF (ABS(SP(3)-P3).LT.PDir) 

GO TO 90 

<iATISFACTDRY Ml FOUND - STORE CURRENT VALUES OF VELOCITYf 
?eIiPERATURE, DENSITY, FRICTION PARAMETER, MACH NUMBER. AND L^ 

CALCULATE T* 

180 V(J.l) = VV 

T(J.1)= TT-CONV( 1, lUNITS) 

RHO(J.l) = RRHO 
FFU.l) =• B1 
MACHU.l) = Ml 
LSTAR(J) = LiSTAR 

TSTAR = TT/E02(G,M1) /. 5/ (G+1. ) 

RETURN 


END 


$1BFTC DISTS 

PRESSURE^ TEMPERATURE, VELOCITY, MACH NUMBER, DENSITY, AND 
FRICTION PARAMETER PROFILES VS X 

PROGRAM VARIABLES 
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X - DISTANCE ACROSS THE FACE OF THE SEAL 
X - POINT COUNTER FOR VARIABLES OTHER THAN X AND P 
MX - MACH NUMBER AT X 

LX - DISTANCE FROM X TO POINT OF CHOKING 
BX - FRICTION PARAMETER AFLX/D 
MTEST - NEW VALUE OF MX IN ITERATIVE SOLUTION OF EQ. 23 
TT - TEMPERATURE AT X 
VV - VELOCITY AT X 
RRHD - DENSITY AT X 
P - PRESSURE AT X 

I - POINT INDEX FOR X AND P 
II - KMODULO 10) - WHEN 1 1 = 1, INCREASE K 
J - FILM THICKNESS INDEX 

SUBROJTINE DIST(J) 

REAL LX, MX, MACH, MTEST, LOSS, LSTAR 

COMMON/ARRAYS/X ( 101) ,P( 20, 101) ,MACH( 20,11), V(20,11),TI20,11) , 

1 RH0(20,11) ,FF( 20, 11) , H(20 ) , XCBAR( 20) ,LSTAR( 20) ,PSTAR ( 20) , 

2 REP(20),FRICT(20),POWER(20),FORC£(20),STIFF(20),FBAR(20) 
COMMON/CONSTS/G,RDIF ,XLAM,C0NLAM,XTURB,C0NTRB,T0,P0,P3,PT0L, 

1 RGAS, LOSS, RELAM, RETURB,FMU,TSTAR 

COMMON/CONVRT/CONV(12,2) ,IUNITS 
E02(G,X) » l./(l.-«-.5*(G-l, )«X»»2) 

E06(G,B,X) = l./SQRT ( G*B+l.-.5*(G-»-l« )*ALOG( (G+1, )*X**2/ 

1 (2.+(G-l. )*X**2) ) ) 

K=1 

MX= MACH(J,l) 

DO 120 1=2,101 

11= MDD(I,10) 

IF (H( J).LT.H(l) ) GO TO 90 
IF (I.NE.lOl) GO TO 90 
MX = 1. 

LX = 0. 

BX = 0. 

GO TO 110 

90 MX= MX*(1.0A+.2E-3*FL0AT(I ) ) 

IF (MX.GT.l.) MX=1. 

LX* ROIF-Xm 

IF (H( J).LT.H(l) ) LX=LSTAR(J)-Xm 
BX* 2.*FRICT( J)*LX/H( J) 

100 MTEST* EQ8(G,BX,MX) 

IF(ABS(MTEST-MX) .LT. l.E-7) GO TO 110 
MX* MTEST 
GO TO 100 

110 TT * TSTAR*,5*(G+1. )*EQ2(G,MX) 

116 VV* MX*SQRT (G*RGAS*TT*CONV( 11, lUNITS ) ) 

RRHO* RH0(J,1) *V(J,1)/VV 

P(J.I)= RRHO*RGAS*TT/CONV( 12,IUNITS) 

C 

C STORE EVERY 10-TH VALUE OF TEMPERATURE, DENSITY, VELOCITY, MACH 

C NUMBER, AND FRICTION PARAMETER 
C 

IFdl.NE.l) GO TO 120 
K* K+l 

T(J,K)* TT-C0NV(1,IUNITS) 

V(J,K)» VV 
MACHIJ.K)* MX 
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FF(J.<)= BX 
RHO(J,K)= RRHO 
120 CONTINUE 
C 

RETURN 

END 


SIBFTC PX 


INTEGRAND OF FORCE INTEGRAL 

FUNCTION FFUNC(JtX) 
COMMON/CONSTS/D(8) »P3.DD(7> 

FFUNC»PRESSU,X)-P3 

RETURN 

END 


IBFTC PXX 


INTEGRAND OF CENTER OF PRESSURE INTEGRAL 

FUNCTION XCFUNCUtX) 

COMMON/CONST S/D ( 8 ) ,P3 t DD (7 ) 

XCFUNC = (PRESS! J,X)-P3)*X 

RETURN 

END 
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ilBFTC PRSS 


C 

C INTERPOLATION ROUTINE FOR PRESSURE 

C 

C CALL VECTOR VARIABLES 

C «••••»•••••••«••«»••• 

c 

C J - FILM THICKNESS INDEX 

C XX - DISTANCE ACROSS FACE OF SEAL 

C 

FUNCTION PRESSIJ.XX) 

REAL MACHtLSTAR 

COMMON/ARRAYS/X( 101) ,P(20, 10 1 ) , MACHI 20, 1 1 ) . V (20, 1 1 ) , T ( 20 ,11 ) , 

1 RHO(20,ll),FF(2C,li ) ,H(20),XCBAR(20),LSTAR120) ,PSTAR(20) , 

2 REP(20),FRICT{20),POWER(20),FORCE(20),STIFF(20»,FBAR(20) 

C 

DX = X(2)-Xm 
DO iOO 1=1,101 
II = I 

IF (XX-X(in 110,120,100 
100 CONTINUE 
110 IF (II.LT.2) 11=2 

IF (II.GT.lOO) 11=100 

PRESS = (P(j,ii-i)*(xx-x{iin*(xx-x(ii+i))-p(j,in* 

1 (XX-X( I I-l ) )*( XX-X ( I I+l) )*2.*P( J,II+1)»(XX-X(II-1) )♦ 

2 (XX-XI II n )/2./DX*»2 
RETURN 

120 PRESS = P(J,II) 

RETURN 

END 


MBFTC OERIV 


LAGRANGE NUMERICAL DIFFERENTIATION OVER MAXIMUM OF 5 POINTS 
TO DETERMINE AXIAL FILM STIFFNESS 


CALL VECTOR VARIABLES 


MMAX - NUMBER OF FILM THICKNESSES 
PROGRAM VARIABLES 


XX - FILM THICKNESSES 
YY - FORCE 

DDY - AXIAL FILM STIFFNESS 

X - WORKING ARRAY FOR FILM THICKNESSES 
Y - WORKING ARRAY FOR FORCE 
DY - WORKING ARRAY FOR AXIAL FILM STIFFNESS 
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C 

C IND - INDEX OF TYPE OF FLOW 

C IND = 1 MEANS CHOKED LAMINAR FLOW 

C IND = 2 MEANS CHOKED TURBULENT FLOW 

C IND = 3 MEANS NON-CHOKED LAMINAR FLOW 

C IND = 4 MEANS NON-CHOKED TURBULENT FLOW 

C 

C H - POINT INDEX 

C MM - POINT COUNTER FOR EACH TYPE OF FLOW 

C 

C #***»» 

C N » 

C 1ST * 

C I N » 

C K » - AUXILLIARY INDEXES IN SETTING UP MATRIX A 

C I I * 

C JJ * 

C J » 

c ***»*» 

C 

C A - MATRIX OF DIFFERENCES DEFINED BY EQ,- 

C 

c «**«»« 

C SI * 

C S2 » - SUMS AND PRODUCTS DEFINED BY EQS. - 

C PI * 

C P2 » 

C «««•»» 

C 

C KY - DERIVATIVE (4XIAL FILM STIFFNESS) ROUNDED OFF TO 

C NEAREST WWH.E NUMBER 

C 

SUBROUTINE STFNSS(MAX) 

DIMENSION X(20), Y(20), DY ( 20) , A ( 5t 5 ) 

COMMOM/CONSTS/Dl ( 12 )» RELAM , RETURB. FMU, TSTAR 
COMMO>) / ARRAYS/D2 ( 322 1 ), X X( 20 ), 03(60 ), REP (20) ,04(40) ,YY( 20) , 

1 DDY(20),05(20) 

ELIMINATE INVALID POINTS AND ARRANGE VALID POINTS IN ASCENDING 
ORDER, IF THERE ARE LESS THAN 2 VALID POINTS, NO DIFFERENTIATION 
IS POSSIBLE, 

DO 90 M=1,MAX 
90 CDY(M) * 0, 

DO 400 IND=1,4 
PM = 3 

100 00 113 M=1,MAX 

IF (XX(M),GE.1.E30) GO TO 110 
GO TO (1,2, 3, 4), IND 

1 IF (XX(M).GT,XX(1),AN0,REP{M).LE,RELAM) GO TO 105 
GO TO llO 

2 IF (XX(M),GT.XX(1).AND.REP(M),GE,RETURB) GO TO 105 
GO TO 110 

3 IF (XX(M),LE,XX( 1),AND.REP(M),L£. RELAM) GO TO 105 
> GO TO 110 

4 IF (XX(M).LE,XX(1),AN0,REP(M).G£.RETURB) GO TO 105 
GO TO 110 

105 PM = MM41 
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XCMM) » XX(M) 

Y(MM) - YY(M) 

110 CONTIMUE 

IF (MH.LT.2) GO TO 400 
130 CALL S0RTXY(X,Y,MM) 

SET UP MATRIX OF X DIFFERENCES FOR EACH POINT X(K» 

200 N = MIN0(MM,5) 

DO 250 K=1,MM 
1ST * MAX0(K-2.1) 

1ST » MINO(MM-N+ltIST) 

IN = IST+N-1 

DO 211 II=IST,IN 

I = II-IST+i 

DO 210 JJ=IST,IN 

J = JJ-IST+1 

A(I,J) * X( II )-X( JJ ) 

210 CONTINUE 

211 CONTINUE 

FORM SUMS AND PRODUCTS FOR DERIVATIVE FORMULA 

220 SI = 0. 

S2 = 0. 

P2 = 1. 

DO 231 II=IST,IN 
IF (lI.EQ.KI GO TO 231 
I = II-IST+1 
PI = X(in-X(K) 

S2 = 52-1. /PI 
P2 = P2*P1 
DO 230 JJ=1,N 

IF (I.NE.JJ) P1=P1*A(I,JJ) 

230 CONTIMUE 

SI = S1*Y(II)/P1 

231 CONTINUE 

IF (fN/2)*2.NE.N» S2=-S2 

DERIVATIVE 

KY = S2*Y(KJ+P2*S1 
DY(K) » KY 
250 CONTINUE 

PUT CALCULATED DERIVATIVES IN ORDER TO CORRESPOND TO INPUT XX 
ARRAY 

300 DO 320 M^^ltMAX 
DO 310 II^ltMM 

IF (XX(M).NE.Xnin GO TO 310 
IF I(N/2)*2.NE.N) GO TO 311 
COY(M) = -DYUI) 

GO TO 320 

311 CDY(M» = DY(II) 

GO TO 320 
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310 CONTINUE 
320 CONTINUE 
400 CONTINUE 

RETURN 

END 


SIBFTC SRT 


ARRANSE POINTS (X,Y) IN ORDER OF ASCENDING X 
CALL VECTOR VARIABLES 


X - INDEPENDENT VARIALBE 

y - DEPENDENT VARIABLE 

N - NUMBER OF POINTS TO BE SORTED 

PROGRAM VARIABLES 


T - TEMPORARY STORAGE 
I - POINT INDEX 
J - POINT INDEX 
NN - N-1 


SUBROUTINE SORTXXIX, Y,N) 
DIMENSION X(IOO), Y(IOO) 

NN= N-l 

DO 110 I>1,NN 
II = I 

DO 100 J*2,N 

IFIXU) .GE. xm) GO TO 100 
T= X(J) 

xu):= xm 
xm= T 
T= YCJ) 

YIJ)= YII) 

Y(I)= T 
100 CONTINUE 
110 CONTINUE 

RETURN 

END 
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tlBFTC SMPS^2 


C 

C NUMERICAL INTEGRATION BY SIMPSONS RULE 

C 

C INTEGRAL OF Y OX FROM XMIN TO XMAX » ( H/3 ) * ( YO+4Y1+Y2 ) 

C WHERE YO = Y EVALUATED AT XMIN 

C Yl = Y EVALUATED ATA ( XMIN+XMAX ) /2 

C Y2 = Y EVALUATED AT XMAX 

C H = (XMAX-XMINI/2 (STEP SIZE) 

C 

C CALL VECTOR VARIABLES 

C ********************* 

c 

c J - INDEX OF FILM THICKNESS 

C XMAX - LOWER LIMIT OF INTEGRATION 

C XMAX - UPPER LIMIT OF INTEGRATION 

C FUNC2 - FUNCTION SUBPROGRAM TO EVALUATE Y 

C KER - NUMERICAL CONSTANT TO INDICATE IF INTEGRATION IS ACCURATE 

C 

C PROGRAM VARIABLES 

C **#**••»•*»**•*** 

c 

c V - INDEPENDENT VARIABLE IN INTEGRATION 

C H - STEP SIZE 

C A - YO 

C B - Yl 

C C - Y2 

C P - 3*(VALUE OF INTEGRAND) 

C E» - DIFFERENCE BETWEEN ANSWERS USING DIFFERENT STEP SIZES 

C NE* 

C **• 

C ANS - ACCUMULSTED ANSWER FOR MANY SUBINTERVALS 

C N - SUBINTERVAL COUNTER (N GE 200 MEANS INTEGRAL IS 

C INACCURATE) 

C T - ERROR TOLERANCE 

FRAC - FRACTION OF ERROR TOLERANCE APPLICABLE TO N(TH) 
SUBDIVISION 


test* - TEST VALUE FOR ERROR IN INTEGRAL 
NTEST* 


Q - TEST VALUE OF FINAL ANSWER 
FUNCTION SIMPS2U,XMIN,XMAX,FUNC2fKER) 

DIMENSION V( 200).H( 200). A( 200), 8(200) fC (200) »P(200) »E (200) ,NE(200) 
EQUIVALENCE ( E , NE ),( TEST .NTEST ) 

DEFINE STARTING VALUES 


T=3.0E-5 

V(1)=XMIN 

H(1)=3,5*(XMAX-XMIN) 
A(1)=FUNC2( J.XMIN) 

B(l)*=FUNC2( J,XMIN+H( 1) ) 
C(1)=FUNC2( J.XMAX) 
P(l)=H(l)*(A(l)+A,0*B(l)+C(l)) 
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Ed )=P(1) 

ANS=P(1) 

N=1 

FRAC=2.0*T 

BEGIN INTEGRATION USING MORE SUBINTERVALS WHERE VALUE OF INTEGRAND 
IS CHANGING RAPIDLY 

1 FRAC=3.5*FRAC 

2 TEST=ABS(FRAC*ANS) 

K=N 

DO 7 1=1, K 

IF INFEST.GT.IABSINEI I) ) I GO TO 7 
N = N«^l 

V{N)=V( I)+H( I ) 

H(N)=0.5*H( I ) 

A(N)=Bd ) 

B(N)=FUNC2{ J,V(N)+H(N) ) 

C(N)=C ( 1 ) 

P(N)=H(N)*(A(N)+4,0*B(N)+C(N» ) 

0=P(I) 

H( I )=H(N) 

B(I )=FUNC2( J,V( n+H( I ) ) 

C(I)=A(N> 

Pin=H(I)*(A{I)+4.0*B( n+Cd )) 

Q=Pd)+P(N)-Q 
ANS=AVS+0 
E(H=3 
E{N»=3 

IF (N.GE.2CO) GO TO 13 
7 CONTINUE 

IF (N.GT.K) GO TO 2 
0 = 0. D 
DO U 1=1, N 

11 C=Q+E(I) 

12 IF (ABSIQ)-T*ABS( ANS ) ) 14,14,1 

13 KER=KER+2 
C 

C ACCUMULATE FINAL ANSWER 

C 

14 ANS=0.0 

DO 16 1=1, N 
16 ANS = ANS+Pm 

SIMPS2=(ANS+0/30.0)/3.0 

C 

RETURN 

END 
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SIBFTC GRAFS 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


MICROFILM PLOTTING ROUTINE 


CALL VECTOR VARIABLES 


PTLSKP 

NJ 

TITLE 


ARRAY OF LOGICAL VARIABLES THAT TELL 
NUMBER OF FILM THICKNESSES 
ALPHA-NUMERIC IDENTIFICATION OF DATA 


WHICH PLOTS TO OMIT 
= PLOT TITLE 


PROGRAM VARIABLES 

•****»»**»*»»*•*• 


HLES 

PWRLES 

XCLES 

FLE3 


FILM THICKNESS LEGEND FOR 
POWER LEGEND FOR ORDINATE 
CENTER OF PRESSURE LEGEND 
FORCE LEGEND FOR ORDINATE 


ABSCISSA OF FIRST 3 PLOTS 
OF FIRST PLOT 

FOR ORDINATE OF SECOND PLOT 
OF THIRD PLOT 


XLE3 

PLE3 

TLES 

RLE3 

MLES 

BLEG 


X LEGEND FOR ABSCISSA OF LAST 5 PLOTS 
PRESSURE LEGEND FOR ORDINATE OF FOURTH PLOT 
TEMPERATURE LEGEND FOR ORDINATE OF FIFTH PLOT 
DENSITY LEGEND FOR ORDINATE OF SIXTH PLOT 
MACH NUMBER LEGFND FOR ORDINATE OF SEVENTH PLOT 
FRICTION PARAMETER LEGEND FOR ORDINATE OF EIGHT PLOT 


NORMALIZING PARAMETERS 


C 

C 

C 


HMAX - MAXIMUM FILM THICKNESS 
PWRMAX - MAXIMUM POWER 
RDIF - FLOW LENGTH 

P3 - SUPPLY PRESSURE 

TO - TOTAL TEMPERATURE 

RT - DENSITY OF SUPPLY RESERVOIR 

XX - AUXILLIARY ARRAY FOR INDEPENDENT VARIABLE IN EACH PLOT 
YY - AUXILLIARY ARRAY FOR DEPENDENT VARIABLE IN EACH PLOT 
IJ - POINT COUNTER 

SUBROUTINE GRAF IC (PLTSKP ,NJ t TITLE ) 

LOGICAL PLTSKP 

REAL MACH.LSTAR 

DIMENSION PLTSKP ( 8 ) t T ITLE( 12 ) f XX ( 25 ) f YY( 200 ) , , m 

D I MENS ION HLEG( 10 ) f PWRLEGI 10 )»XCLEG( 101 »XLEG( 10) »PLEG( 10) t 

1 TLEG(10),RLEG( 10) tMLEG(10)f BLEG(10)fFLEG(10) 

rnMMnM/ARRAYS/XIlOl) tP(20*l01) » MAC HI 20 1 11)»V120»11) tT(20»ll# i 

C0MM0V/ARRAYS/X(lOl).PI/^.^^^^^^^^^^^^^^^^^l^^^^^j20),PSTAR(20), 

2 REPI20) ,FRICT( 20 ) t POWER! 20 ) .FORCE ( 20 ) i STIFF ( 20) .FBARI20) 

COMMON/CONSTS/G,RDIF,XLAM.CONLAM.XTURB.CONTRB.TO,PO,P3,PTOL, 

1 rGAS, LOSS, RELAM, RETURB.FMUO.TSTAR 

C0MM0'I/C0NVRT/C0NV(12,2) , I UN ITS 


ALL PLOTS ARE DIMENSIONLESS 


HMAX » HID 
PWRMAX = POWER! 1) 
DO 100 J=2,NJ 
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^ HMAX » AMAXi.{HMAX,H( J)) 

PWRMAX = AMAX1(PWRMAX,P0WER( J) ) 
ICO CONTHUE 

DO UD J = 1,NJ 
110 XX(J> = H(J)/HMAX 
C 

IF(PLrSKP(l) ) GO TO 210 
DO ?0D J=1,NJ 

200 VY(J) = POWERC J )/PWRMAX 
C 


C**** CALL PLOT ROUTINE TO PLOT POWER VS FILM THICKNESS ************ 

Q 

210 IF (PLTSKPI2)) GO TO 220 
C 

C***» CALL PLOT ROUTINE TO PLOT CENTER OF PRESSURE VS FILM THICKNESS 

220 IF (PLTSKPI3)) GO TO 300 
C 

C**** CALL PLOT ROUTINE TO PLOT FORCE VS FILM THICKNESS *****•♦•*»** 
C 

C MULTIPLE PLOTS OF PARAMETERS THAT VARY WITH X 

C 


300 IJ = D 

DO 301 J=1,NJ 
DO 301 1=1,101,10 
IJ = IJ+1 

301 XX(IJ) = Xm/ROIF 

IF (PLTSKPI4)) GO TO 310 
IJ = 3 

DO 302 J=1,NJ 
DO 302 1=1,101,10 
IJ = IJ+1 

302 YY(IJ) = P(J,1)/P0 
C 

C**«* CALL PLOT ROUTINE TO PLOT PRESSURE VS DISTANCE 
C 


310 IF (PLTSKPISn GO TO 320 
IJ = 0 

DO 311 J=1,NJ 
DO 311 1=1,11 
IJ = IJ+I 

311 YY(IJ) = CT(J,n+CONVIl, lUNITS) )/T0 
C 

C***« CALL PLOT ROUTINE TO PLOT TEMPERATURE VS DISTANCE 

0 

320 IF (PLTSKPI6)) GO TO 330 
IJ = 0 

RO = P0*C0NV(12,IUNITS)/RGAS/T0 
DO 321 J=1,NJ 
DO 321 1=1,11 
IJ = lJ+1 

321 YY(IJ) - RHOIJ,I)/RO 
C 

C***» CALL PLOT ROUTINE TO PLOT DENSITY VS DISTANCE *** 
C 

330* IF tPLTSKPIT)) GO TO 3A0 
IJ = 3 

DO 331 J=1,NJ 


»«« 
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00 ^31 
IJ = IJ+1 

331 YY(IJ) = MACH(J.I) 

C**** CALL PLOT ROUTINE TO PLOT MACH NUMBER VS DISTANCE ***••* 
C 

340 IF (PLTSKP(8)) RETURN 
IJ = 3 

DO 341 J=lfNJ 
DO 341 1=1 »il 
IJ = IJ+1 

341 YY(IJ) = FF(J,I) 

C***» CALL PLOT ROUTINE TO PLOT FRICTION PARAMETER VS OISTNACE 
C 

RETURN 

END 


*I8FTC BDATA 


C 

C 

C 


C 


UNITS CONVERSION DATA 


BLOCK DATA 

COMMON/CONVRT/CONVI 12f 2) , I UN ITS 
DATA (CONV( If l)f 1=1. 12)/0.t .14591E-5. 110. 
1 .5051554E-2. 6*1. / 

DATA ( CONV( I .2) . I=lf 12)/ 460.. 1.57639E-10. 

1 13.405833. 13.083. 1728.. .021818182. 

2 4.4756636 / 


33333. 

198.6. 

42.42. 


2 * 1 .. 

720, . 
3.3E4, 


32.174 


END 
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Comment: Use method of iteration to solve flow equations for entrance Mach number, friction factor, and 
pressure Reynolds number when flow is exactly choked (i.e., M2 -1, P2 • P3. and L • R2 - Rl) 
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CHOKE 
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SAMPLE PROBLEM 


OUASI-ONE DIMENSIONAL COMPRESSIBLE FLOW SEAL PROGRAM 

SAMPLE PROBLEM 

INTERNATIONAL UNITS 

INPUT DATA - 




Rl.METERS 

0.829E-01 

PO.NSMA 

0.44BE+06 

MOLECULAR HEIGHT 
0.290E-01 

F»K/RE*«N 

R2 •METERS 
0.B42E-O1 

P3*NSMA 

0.103E4^06 

CP, JOULES/KG-OEG K 
O.IOOE«04 

KILAMINAR) 

24.00C 

FLOW LENGTH. METERS 

0. 

PD/P3 

0. 

GAMMA 

1.400 

NILAHINARt 

0.10E401 

FLOW UlOTHiMETERS 

0. 

TOtDEG K 
310.9 

VISCOSITV,N-SEC/MZ 

0 

UPPER LIMIT RE (LAMINAR) 
2300 .C 


LOSS COEF. 
1.00 

SPEED,RPS 

0. 

K(TURBULENT) 

0.7900E-C1 



V, M/SEC 
60.96 

NITURBULENT ) 
0.2500 




LOWER LIMIT RE (TURBULENT) 
3000. C 

OUTPUT DATA - 




Rl.METERS 

0.929E-01 

PO.NSMA 

0.448E«06 

MOLECULAR WEIGHT 
0.029 

CALCULATE 

R2, METERS 
0.842E-01 

P3.NSMA 

0.103E406 

CP. JOULE S/KG-DEG K 
0.100E*04 

POWER 

DIMENSIONLESS CUANTITIES 

FLOW LENGTH. METERS 
0.127E-02 

P0/P3 

4.333 

GAMMA 

1.400 

PRESSURE OISTPIBUTIONS 

FLOW WIDTH. METERS 
0.525E+00 

TO. DEG K 
310.928 

VISCOSITY. N-SEC/H2 
0.1899E-04 

PLOT 

A RE A. M2 
0.667E-03 

GAS CONSTANT. JOULES/KG-OEG K 
0.287E^06 

SPEED, RPS 
116, 10 

POWER XC(BAR) 

FORCE PRESSURE 


LOSS COEF. 
1.00 

V.H/SEC 

0.610E*02 

TEMPERATURE DENSITY 


MACH NUMBER AFL/C 


* - CHOKING FILM THICKNESS 

♦ - TRANSITION REGION 
/ - TURBULENT FLOW 


QUASI -ONE DIMENSIONAL COMPRESSIBLE FLOW 


H THICKNESS 

M(DQT) 

0 

METERS 

KG/SEC 

SCMS 

0.73690E-04 

0.549E-03 

0.277E-05 

0.25400E-05 

0.317E-07 

0.160E-09 

0.50800E-05 

0.253E-06 

C.128E-08 

0.76200E-05 

0.855E-06 

0.432E-08 

0.10160E-04 

0.203E-05 

0.102E-07 

0.12700E-04 

0.396E-D5 

0.200E-07 

0.15240E-04 

0.683E-05 

0.345E-07 

0.177B0E-04 

0,rOBE-04 

0.548E-07 

0.20320E-04 

0.162E-04 

0.81.6E-07 

0.22860E-04 

0.230E-04 

0. 116B-06 

0.254O0E-04 

0.314E-04 

C.159E-06 

0.27940E-04 

0.416E-04 

0.210E-06 

0.30480E-04 

0.537E-04 

0.272E-06 

0. 3302 OE -04 

0.679E-04 

0. 343E-06 

0.35560E-04 

0.841E-04 

0.425E-06 

0.38100E-04 

0.1026-03 

C.51BE-06 

0.40640E-04 

O.123E-03 

C.621E-06 

0.43180E-04 

0.1466-03 

0.735E-D6 

0.45720E-04 

0.170E-03 

0.860E-06 


SEAL PROGRAM 

SAMPLE PROBLEM 

KNUDSEN 

LAMBDA 

STIFFNESS 

NUMBER 

METERS 

N/M 

0.266E-01 

0.196E-05 

-0. 7436^05 

0.855 

0.217E-05 

-0.4126*^04 

0.427 

0.217E-05 

-0.177E*04 

0.285 

0.217E-05 

-589.0 

0.214 

0.217E-05 

-83S.0 

0.171 

O.217E-05 

-0.209E+04 

0.142 

0.217E-05 

-0.346E*-04 

0.122 

0.217E-05 

-0.524E*04 

0.107 

0.217E-05 

-0.779E+04 

0.947E-01 

0.216E-05 

-0.109E*05 

0.851E-01 

0.216E-05 

-0. 147E*05 

0.773E-01 

0.216E-05 

-0.191E+05 

0.707E-01 

0.215E-05 

-0.239E+05 

0.651E-01 

0.215E-05 

-0.291E^05 

0.602E-01 

0.214E-05 

-0.3466*05 

0.560E-01 

0.213E-05 

-0.401E+05 

0.523E-01 

0.213E-05 

-0.454E*^05 

0.490E-01 

0.212E-05 

-0.504E+05 

0.461E-01 

0.211E-05 

-0.551E*05 


INTERNATIONAL UNITS 


FORCE XC FCRCE 

N METERS EAR 

141*7 0.504E-03 C.616 

130.9 0*473E-03 0.604 

13B.9 0.474E-03 C.604 

130.9 0.474E-03 0.6C4 

138.9 0.474E-03 C.60A 

138.9 0.474E-03 C.604 

138.9 0.474E-03 C,60^ 

130.9 0.474E-03 0.604 

138.9 0.474E-03 C.60A 

138.9 0.474E-03 0.604 

139.0 0.474E-03 C.605 

139.0 0.475E-03 0.605 

139.1 0.475E-03 C.605 

139.1 0.476E-03 0*605 

139.2 0.476E-03 C.606 

139.3 0.477E-03 C.606 

139.4 O.470E-C3 C.607 

139.6 0.479E-C3 C.607 

139.7 0.481E-03 0.608 


XC 
BAR 
0.397 
0.373 
0. 373 
0.373 
C.373 
0.373 
0.373 
0.373 
0.373 
0.373 
0.373 
0. 374 
0.374 
0.375 
0.375 
0.376 
0.377 
0.378 
0.37S 
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FILM THICKNESS 

FRICTION 

REiP» 

RE<R) 

METERS 

FACTOR 



• 0.73690E-04 

0.216 

111.3 

1*812 

0.25400E-05 

0.378E404 

0.635E-02 

0.594E-01 

0.50800E-05 

472.2 

0.508E-01 

0.119 

0.76200E-05 

139.9 

0*172 

0*178 

0.10160E-04 

59.04 

0*406 

0.238 

0.12700E-04 

30.24 

0*794 

0.297 

0.15240E-04 

17.51 

1*371 

0.357 

0.1778DE-04 

11*04 

2*174 

0.416 

0.20320E-04 

7*405 

3.241 

0*476 

0.22860E-04 

5.212 

4.605 

0.535 

0.25400E-04 

3.811 

6.298 

0.595 

0.27940E-04 

2.874 

8*351 

0.655 

0. 30480E-04 

2.225 

10.79 

0*715 

0.33020E-04 

1*761 

13.63 

0.776 

0.35560E-04 

1.421 

16,89 

0*837 

0.38100E-04 

1*166 

20.58 

0.898 

0.40640E-04 

0.971 

24.71 

0*960 

0.43180E-04 

0.820 

29.26 

1*022 

0.45720E-04 

0.701 

34.24 

1.085 


P(*» 

LC*I 

POWER 

HEAT LOSS 

NSMA 

METERS 

WATTS 

WATTS 

0.103E406 

0.127E-02 

0.627 

0.627 

173.1 

0.134E-02 

18*53 

18.53 

692.4 

0.134E-02 

9.263 

9.263 

0.156E404 

0.134E-02 

6.175 

6.175 

0.277E404 

0.134E-02 

4.632 

4.632 

0.432E404 

0.134E-02 

3*705 

3* 705 

0.622E404 

0.134E-02 

3*088 

3.088 

0.846E404 

0.134E-02 

2.646 

2*646 

0.110E405 

0.134B-02 

2.315 

2.315 

0.139E405 

0.133E-02 

2.058 

2.058 

0.172E405 

0.133E-02 

1*852 

1.852 

0.207E405 

0.133E-02 

1*683 

1.683 

0.245E405 

0.132E-02 

1*542 

1.542 

0.285E405 

0.132E-02 

1.423 

1.423 

0.328E405 

0.132E-02 

1*321 

1.321 

0.373E405 

0.131E-02 

1.232 

1,232 

0.420E405 

0.131E-02 

1.154 

1. 154 

0.46BE40S 

0.130E-02 

1.086 

1.086 

0.517E405 

0.130E-02 

1.025 

1.025 


OEUT) 

DEC K 

I. 137 
0.S82E«04 
0.364E405 
0.719E«04 
0.227E404 

932.0 
M9.7 

243.0 
142.6 
89.22 
58.71 
40*25 
26*56 
20*86 
15*63 

II. 97 
9*343 
7.422 
5*988 


FILM THICKNESS - 0.73690E-04 METERS 





X 

P 

MACH 

4FL/0 

DENSITY 

V 

T 

METERS 

NSMA 

NUMBER 


KG/M3 

M/SEC 

DEG K 

0 

0.427E*06 

0.263 

7.433 

0.485E-02 

0.292E*04 

306.7 

0.127E-03 

0.410E«^06 

0.275 

6.690 

0o466E-02 

0.305E«04 

306.3 

0.254E-03 

0.391E406 

0.287 

5.947 

0.446E-02 

0.318E^04 

305.9 

0.381E-03 

0.372E^06 

0.3C2 

5.203 

0.424E-02 

0.335E^04 

305.4 

0.508E-03 

0.351E^06 

0.320 

4.460 

0.401E-02 

0.354E^04 

304*7 

0.635E-03 

0.328E406 

0.341 

3.717 

0.376E-02 

0.377E^04 

303*9 

0.762E-03 

0.304E406 

0*368 

2.973 

0.349E-02 

0.406E4-04 

302.7 

0.889E-t>3 

0.276E4-06 

0.4C4 

2.230 

0.319E-02 

0.445E+04 

301.1 

0. 102E-02 

0*243E«06 

0.457 

1.487 

0.284E-02 

0.500E«04 

298.5 

0.114E-02 

0.201E406 

0*547 

0.743 

0.239E-02 

0. 594E4^04 

293.4 

0.1Z7E-0Z 

0. I03E^06 

1.000 

0 

0.139E-02 

0«102E^05 

259.1 

FILM THICKNESS - 0.25400E-C5 METERS 





X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

METERS 

NSMA 

NUMBER 


KG/M3 

M/SEC 

DEG K 

0 

0.448E+06 

0.423E-03 

0.399E+07 

0.502E-02 

4.730 

310,9 

0.127E-03 

0.426E*06 

0.445E-03 

0.361E-^07 

0.47BE-02 

4.971 

310*9 

0.254E-03 

0.404E406 

0.470E-03 

0.323E+07 

0.452E-02 

5.253 

310*9 

0.3B1E-03 

0.379E+06 

0.500E-03 

0.286E^07 

0.425E-02 

5.590 

310*9 

0.508E-03 

0.353E406 

0.537E-03 

0.248E^07 

0.396E-02 

6.000 

310.9 

0.635E-03 

0.325E406 

0.583B-03 

0.210E+07 

0.364E-02 

6.518 

310.9 

0.762E-03 

0.295E«06 

0.644E-03 

0.172E+07 

0.330E-02 

7.197 

310,9 

0. 809E-t)3 

0.260E«06 

0.729E-03 

0.135E+07 

0.292E-02 

8.144 

310*9 

0.102E-02 

0.221E«^06 

0.859E-03 

0.968E+06 

0.247E-02 

9.603 

310.9 

0. 114E-02 

0.I72E-^06 

O.llOE-02 

0.590E+06 

0.193E-02 

12.30 

310.9 

0.127E-02 

0,103E+06 

0.183E-02 

0.2I2E>06 

0.116E-02 

20.51 

310.9 

FILM THICKNESS « O.5O0OOE-O5 METERS 





X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

METERS 

NSMA 

NUMBER 


KG/H3 

M/SEC 

DEG K 

0 

0.44BE>06 

0.169E-02 

0.249E+06 

0.502E-02 

18.92 

310.9 

0.127E-03 

0.426E«06 

0.178E-02 

0.226E+06 

0.478E-02 

19.88 

310.9 

0.254E-03 

0.404E«06 

0.188E-02 

0.202E«06 

0.452E-02 

21.01 

310.9 

0.381E-03 

0.379E^06 

0.200E-02 

0.179E+06 

0.425E-02 

22.36 

310.9 

0.508E-03 

0,353E^06 

0.215E-02 

0.I55E4-06 

0.396E-02 

24.00 

310.9 

0.635E-03 

0.325E+06 

0.233E-02 

0.131E^06 

0.364E-O2 

26.07 

310.9 

0.762E-03 

0.295E+06 

0.257E-02 

0.108E<t-06 

0.33OE-O2 

26.78 

310.9 

0. B89E-03 

0.260E406 

0.291E-02 

0.841E+05 

0.292E-02 

32.57 

310.9 

0.102E-02 

0.221E4-06 

0.344E-02 

0.605E^05 

0.247E-02 

38.41 

310.9 

O.114E-02 

0.172E+06 

0.440E-02 

0.369E^05 

0.193E-02 

49.19 

310.9 

0.127E-02 

0.103E+06 

0.733E-02 

0.133E+05 

0.116E-02 

81.99 

310.9 


FILM THICKNESS * 0.76200E-05 METERS 





X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

METERS 

NSMA 

NUMBER 


KG/ M3 

M/SEC 

DEC K 

0 

0.448E^06 

0.381E-02 

0.493E+05 

0.502E-02 

42.56 

310.9 

0,127E-03 

0.426E^06 

0.400E-02 

0. 4466^05 

0.478E-02 

44.73 

310.9 

O.254E-03 

0.404E+06 

0.423E-02 

0. 3996*05 

0.452E-02 

47.27 

310.9 

O.301E-O3 

0.379E+06 

0.450E-02 

0.353E+05 

0.425E-02 

50.30 

310.9 

0.508E-03 

0.353E+06 

0.483E-02 

0.306E+05 

0.396E-02 

53,99 

310*9 

0.635E-03 

0.325E^06 

0.525E-02 

0.2596405 

0.364E-02 

58.65 

310.9 

0.762E-03 

0.295E+06 

0.5796-02 

0.213E405 

0.330E-02 

64.76 

310.9 

0. 809E-O3 

0.260E-»-06 

0.656E-02 

0. 1666*05 

0.292E-02 

73.28 

310.9 

0.102E-02 

0.221E406 

0.773E-02 

0.119E*05 

0.247E-02 

86.41 

310*9 

0.114E-02 

0.172E+06 

0.990E-02 

0.726E*04 

0.193E-02 

110.7 

310.9 

0. 127E-02 

0.103E406 

0.165E-01 

0.262E*04 

0. 116E-02 

184.4 

310.9 


TORQUE 

0*5406-02 
0.160 
0.798E-01 
0.532E-01 
0.399E-01 
C.319E-01 
0*266 E-01 
0*228E-01 
0.1 99 E-01 
0.177E-01 
0.159E-01 
0.145E-01 
0.133E-01 
0.123E-0I 
C.114E-01 
Q.106E-01 
0.994E-02 
0,935E-02 
0.682E-02 
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FILM THICKNESS » 0,10160E-04 METERS 


X 

p 

MACH 

4FL/D 

DENSITY 

V 

METERS 

NSMA 

NUMBER 


KC/M3 

M/SEC 

0 

0.448E406 

0.677E‘02 

0.156E+05 

0.502E-02 

75.65 

0.127E-03 

0.426E406 

0.711E-02 

0.141E^05 

0.478E-02 

79.51 

0.254E-03 

0.403E+06 

0.752E*0Z 

0.126E+05 

0.452E-02 

84.02 

0. 381E-03 

0.379E*06 

0.800E-02 

0.112E*05 

0.425E-02 

89.41 

0. 508E-03 

0.353E^06 

0.859E-02 

0.966E+04 

0.396E-02 

95.97 

0. 635E-03 

0.325E^06 

0.933E-02 

0.820E*04 

0.364E-02 

104. 2 

0.762E-03 

0.295E^06 

0.1C3E-01 

0.673E+04 

0.330E-02 

115.1 

0.aa9E-03 

0.260E+06 

0.117E-01 

0.525E+04 

0.292E-02 

130.3 

0. 102E-02 

0.221E+06 

0.137E-01 

0.378E+04 

0.247E-02 

153.6 

0.114E-02 

0.1T2E+06 

0.176E-01 

0.230e^04 

0. 193E-02 

196.7 

0.127E-02 

0.103E+06 

0.293E-01 

824.0 

0.116E-02 

327.8 

MLM THICKNESS - 0.12700E-04 METERS 





p 

MACH 

4FL/0 

DENSITY 

V 

METERS 

NSMA 

NUMBER 


KG/M3 

M/SEC 

0 

0.448E+06 

0.106E-01 

0.638E+04 

0.502E-02 

118.2 

0. 127E-03 

0.426E^06 

O.lllE-01 

0.578E+04 

0.478E-02 

124.2 

0.254E-03 

0.403E+06 

0.117E-01 

0.517E+04 

0.452E-02 

131.2 

0.381E-03 

0.379e+06 

0.125E-01 

0.457E+04 

0.425E-02 

139,7 

0.508E-03 

0.353E+06 

0.134E-01 

0.396E+04 

0.396E-02 

149,9 

0.635E-03 

0.325E+06 

0.146E-01 

0.336E+04 

O.364E-02 

162, 8 

0.762E-03 

0.295E+06 

0.161E-01 

0.275E+04 

0.330E-02 

179.8 

0. 889E-03 

0.260E«06 

0.182E-01 

0.215E+04 

0.292E-02 

203,4 

0. 102E-02 

0.221E+06 

0.215E-01 

0.154E+04 

0.247E-02 

239,8 

0. 114E-02 

0.172E^06 

C.275E-01 

939.5 

0.193E-02 

307, 1 

0.127E-02 

0.i03E*06 

0.458E-01 

334.7 

0.116E-02 

511,8 


FILM THICKNESS « 0.15240E 

-04 METERS 





p 

MACH 

4FL/D 

DENSITY 

V 

METERS 

NSMA 

NUMBER 


KG/M3 

H/SEC 

0 

0.448E+06 

0,152E-01 

0.3086^-04 

0. 502E-02 

170. 1 

0. 127E-03 

0.426E*06 

0.160E-01 

0.279E>04 

0-478E-02 

178. 7 

0.254E-03 

0.403E+06 

0.169E-01 

0.249E*04 

0.452E-02 

188. 9 

0.381E-03 

0.379E+06 

0.180E-01 

0.220E+04 

0.425E-02 

201. 0 

0. 508E-03 

0.353E+06 

0.193E-01 

0.191E+04 

0.396E-02 

215.7 

0.635E-03 

0.325E^06 

0.2I0E-01 

0.162E+04 

0.364E-02 

234. 3 

0.762E-03 

0.295E+06 

0.231E-01 

0.133E+04 

0.330E-02 

258.7 

0.889E-03 

0,260E*06 

0.262E-01 

0.103E+04 

0.292E-02 

292.7 

0. 102E-02 

0.221E+06 

0,309E-01 

742,9 

0.247E-02 

345.1 

0.114E-02 

0.172E^06 

0.395E-0i 

451.1 

0.193E-02 

441,8 

0.127E-02 

0.103E+06 

0,659E-01 

159.3 

0.116E-02 

736. 3 

FILM THICKNESS « 0.17780E-04 METERS 




X 

p 

MACH 

4FL/D 

DENSITY 

V 

METERS 

NSMA 

NUMBER 


KG/M3 

H/SEC 

0 

0.448E4'06 

0,2t7E-01 

O.I66E+04 

0.502E-02 

231.3 

0, 127E-03 

0.426E+06 

0,2i7E-01 

0.150E4^04 

0.478E-02 

243.1 

0, 254E-03 

0,403E*06 

0,230E-01 

0.135E+04 

0,452E-02 

256.6 

0.381E-03 

0.379E+06 

0,244E-01 

0,ll9E+04 

0.425E-02 

273 . 3 

0,508E-03 

0,353E«-06 

0,262E-01 

0.103E+04 

0.396E-02 

293.3 

0,635E^03 

0.325E*06 

0.265E-01 

872.7 

0.365E-02 

318.6 

0,762E-03 

0.295E^06 

0.315E-01 

715.0 

0.330E-02 

351.7 

0.889E-03 

0.260E»06 

0,356E-01 

557.4 

0.292E-02 

397.9 

0.102E-02 

0,221E^06 

0.420E-01 

399.7 

O.240E-O2 

469.0 

0, 114E-02 

0,l73Ef06 

0.537E-01 

242.0 

0.193E-02 

600. 2 

0. 127E-02 

0, 1C3E + 06 

0.896E-01 

84.34 

0.116E-02 

0. 100E>04 

FILM THICKNESS - 0.20320E-04 METERS 




X 

p 

MACH 

4FL/0 

DENSITY 

V 

METERS 

NSMA 

NUMBER 


KG/M3 

M/SEC 

0 

0.448E«06 

0,270E-0l 

973.8 

0.502E-02 

301.7 

0* 127E-03 

0.426E+06 

0,2B4E-0l 

881.2 

0.478E-02 

317.0 

0.254E-03 

0.403E+06 

0.3C0E-01 

788,7 

0.452E-02 

335.0 

0,381E-03 

0.379E+06 

0.319E-01 

696.1 

0.425E-02 

356. 4 

0. 508E-03 

0.353E+06 

0.342E-01 

603.5 

0.396E-02 

382.5 

0.635E-03 

0,325E^06 

0.372E-01 

511.0 

0.365E-02 

415.4 

0.76ZE-03 

0.295E^06 

0.410E-01 

410,4 

0.330E-02 

458.5 

0. 889E-03 

0.260E+06 

0.464E-01 

325.8 

0.292E-02 

518.6 

0, 102E-02 

0.221E^06 

0.547E-01 

233.3 

0.248E-02 

611.1 

0, 114E-0Z 

0.173E^06 

0.700E-01 

140.7 

0.194E-02 

781,9 

O.127E-02 

0,i03E406 

0.117 

48,16 

0, 116E-02 

0* 1 30E'«-04 


T 

OEG K 
310.9 
310.9 
310.9 
310.9 
310.9 
310.9 
310.9 
310.9 
310.9 
310.9 
310.9 


T 

OEG K 
310.9 
310.9 
310.9 
310.9 
310.9 
310.9 
310.9 
310.9 
310.9 
310.9 
310.8 


T 

DEG K 

310.9 

310.9 

310.9 

310.9 

310.9 

310.9 

310.9 

310.9 

310.9 

310.8 

310.7 


T 

OEG K 

310.9 

310.9 

310.9 

310.9 

310.9 

310.9 

310.9 

310.8 

310.8 

310.7 
310.4 


T 

DEG K 

310.9 

310.9 

310.9 

310.9 

310.9 

310.8 

310.8 

310.8 

310.7 
310.6 
310.1 



FILM THICKNESS - 0.22860E-04 METERS 


X 

P 

MACH 

4FL/0 

DENSITY 

V 

T 

METERS 

NSMA 

number 


KC/M3 

M/SEC 

DEG K 

0 

0,44BE^06 

0,341E“01 

608,2 

0.502E-02 

381.1 

310.9 

O.127E-0J 

0.426Ei^06 

0,358E-0l 

550.2 

0.478E-02 

400.4 

310.8 

O.254E-0J 

0.403E«06 

0,379E-01 

492.3 

0.452E-02 

423.1 

310.6 

0*3B1E-03 

©•379E+06 

0.403E-01 

434.4 

0.425E-02 

450.1 

310.8 

0.508E-03 

0.353Ei^06 

C.432E-01 

376.5 

0.396E-02 

483.0 

310.6 

0*635E-03 

0*325E+06 

0.4S9E-01 

318.6 

0.36SE-02 

524.4 

310.6 

0.762E-03 

0.295E+06 

0.518E-01 

260.7 

0.330E-02 

578.7 

310.6 

0.eB9E-03 

0.26IE406 

0.586E-OI 

202.8 

0.292E-02 

654.5 

310.7 

0.1O2E-0Z 

0.22IE406 

0.690E-01 

144.9 

0.248E-02 

770.9 

310.6 

O.U4E-02 

0.i73E^06 

0.883E-01 

B6.96 

0.194E-02 

985.6 

310.4 

0.127E-02 

0,103E^06 

0,147 

29.07 

0.116E-02 

0.164E404 

309.6 


FILM THICKNESS - 0.25400E-0A METERS 


X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

METERS 

.NSMA 

NUMBER 


KG/M3 

M/SEC 

DEG K 

0 

0.446E«^06 

0.420E>01 

399.3 

O.S02E-02 

469.2 

310.8 

0.127E-03 

0.426E^06 

0.441E-01 

361.2 

0.477E-02 

492.9 

310.6 

0.254E-03 

0.403E+06 

0.466E-01 

323.1 

0.452E-02 

520.6 

310.6 

0.381E-03 

0.379E+06 

0.496E-01 

285.0 

0.425E-02 

553.9 

310.8 

0.508E-03 

0,353E^06 

0.532E-01 

246.9 

0.396E<-02 

594.4 

310.6 

0.635E-03 

0.32SE^06 

0.577E-01 

206.8 

0.36SE-02 

645.2 

310.7 

0.762E-03 

0.295E^06 

0.637E-01 

170.7 

0.331E-02 

711.9 

310.7 

0.889E>03 

0.261E«06 

0.720E-01 

132.6 

0.292E>02 

804.6 

310.6 

0.102E-02 

0.221E406 

0.848E-01 

94.49 

0.248E-02 

947.5 

310.5 

0.114E-02 

0.173E+06 

0.108 

56.39 

0.194E-02 

0.I21E«04 

310.2 

0.127E-02 

0.103E406 

0.181 

18.28 

0.117E-02 

0.202E«04 

306.9 


FILM THICKNESS « 0,279A0E-04 METERS 


X 

P 

MACH 

4FL/0 

DENSITY 

V 

T 

METERS 

NSMA 

NUMBER 


KG/M3 

M/SEC 

DEG K 

0 

0.447E^06 

0.506E-01 

273.1 

0.502E-02 

565.6 

310.8 

0.127E-03 

0.426Et06 

0.532E-01 

247.0 

O.477E-02 

594.2 

310.6 

0.254E-03 

0.403E^06 

0.562E>01 

220.9 

0.452E-02 

627.7 

310.7 

0.381E-03 

0.379E+06 

0.597E-01 

194.8 

0.425E-02 

667.5 

310.7 

0.508E-03 

0.353E^06 

0.641E-01 

168.6 

0.396E-02 

716,1 

310,7 

0.635E-03 

0.325E«06 

0.696E>01 

142.5 

0.365E-02 

777.2 

310.6 

0.762E-03 

0.295E*06 

0.767E-01 

116.4 

0.331E>02 

857.2 

310.6 

0.889E-03 

0.261E«06 

0.867E-01 

90.24 

0.293E-02 

968.6 

310.5 

0.102E-02 

0.222Ef06 

0.102 

64.11 

0.249E-02 

0.114E^04 

310.3 

0.114E-OZ 

0.173E+06 

0.130 

37.99 

0.195E-02 

0.145E^04 

309.9 

0.127E-02 

0.103E«06 

0.218 

11.86 

0.U7E-02 

0.242E4^04 

308,0 


FILM THICKNESS » 0.30480E-04 METERS 


X 

P 

MACH 

4FL/0 

DENSITY 

V 

T 

METERS 

NSMA 

NUMBER 


KC/M3 

M/sec 

DEG K 

0 

0.447E406 

0.6C0E-01 

193.3 

0.5016-02 

670.0 

310.7 

0.127E-03 

0.42664^06 

0.630E-01 

174.7 

0.477E-02 

703.8 

310.7 

0.254E-03 

0.403E+06 

0.665E-01 

156.2 

0.452E-02 

743.3 

310.7 

0.381E-03 

0,379E*06 

0.7076-01 

137.7 

0.4256-02 

790.3 

310.6 

0.508E-03 

0.353E406 

0.759E-01 

119.1 

0.396E-02 

847.6 

310.6 

0.635E-03 

0.3266*06 

0.823E-Q1 

100.6 

0.365E-02 

919.5 

310.5 

0. 7626-03 

0.295E*06 

0.908E-01 

82.04 

0.3316-02 

0.10IE*04 

310.4 

0.889E-03 

0.261E*06 

0.103 

63.50 

0. 2936-02 

0.1146*04 

310.3 

0.102E-02 

0.222E*06 

0.121 

44.95 

0. 2506- 02 

0.135E*04 

310.0 

0.114E-02 

0.174E*06 

0.154 

26.41 

0.196E-02 

0.172E*04 

309.5 

0.127E-02 

0.1036*06 

0.258 

7.872 

0.117E-02 

0.2866*04 

306.9 

=ILM THICKNESS - 0.33020E-04 METERS 





X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

METERS 

NSMA 

NUMBER 


KG/M3 

H/SEC 

DEG X 

0 

0.447E*06 

0.7C0E-C1 

140.8 

0. 5016-02 

781.7 

310.6 

0,1276-03 

0.425E*06 

0.735E-01 

127.2 

0. 4776-02 

821.0 

310.6 

0.254E-03 

0.403E*06 

0.776E-01 

113.7 

0.452E-02 

866.8 

310.6 

0.381E-03 

0.3796*06 

0.8256-01 

100. 1 

0. 4256-02 

921.4 

310.5 

0.508E-03 

0.3536*06 

0.884E-01 

86.60 

0.396E-02 

987.9 

310.4 

0.635E-03 

0.3266*06 

0.959E-01 

73.05 

0.366E-02 

0.1076*04 

310.4 

0.762E-03 

0.2956*06 

0.106 

59.50 

0. 3326-02 

0.118E*04 

310.2 

0.889E-03 

0.262E*06 

0.119 

45.96 

0.2946-02 

0«133E*04 

310.0 

0.1O2E-0Z 

0.2236*06 

0,140 

32.41 

0.250E-02 

0.156E*04 

309.7 

0. 114E-0Z 

0.1746*06 

0.179 

18.86 

0.1976-02 

0.199E*04 

309.0 

0.127E-02 

0.1036*06 

0,3v'^0 

5.317 

0.1166-02 

0,332E*04 

305.4 



FILM THiC<NESS - 0.35560E-04 METERS 


X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

METERS 

NSMA 

NUMBER 


KG/H3 

H/SEC 

DEG K 

0 

0.446E+06 

0.806E-01 

105.1 

0.501E-02 

900.1 

310.5 

0.127E-03 

0.425E406 

0.846E-01 

94.98 

0.477E-02 

945.1 

310.5 

0.254E-03 

0.402E^06 

O.0^#3E-Ol 

84.63 

0.452E-02 

997.6 

310.4 

0* 381E-03 

0.379E+06 

0.949E-01 

74.68 

0. 4256-02 

0.106E*04 

310.4 

0.508E-C3 

0.353E+06 

0.102 

64.53 

0.3976-02 

0.114E*04 

310.3 

0*635E-03 

0.326E406 

O.llO 

54.38 

0.366E-02 

0.123E*04 

310.2 

0.762B-03 

0,296E*06 

0.121 

44.24 

0.3326-02 

0.136E*04 

310.0 

0.869E-03 

0.262E4^06 

0.137 

34.09 

0.2956-02 

0.153E*04 

309.8 

0.102E-02 

0.2Z3E^06 

0.161 

23.94 

0.25IE-02 

0.179£*04 

309.3 

0.114E-02 

0.175E+06 

0.205 

13.79 

0.198E-0Z 

0. 228E*04 

308.3 

0.127E-02 

O.IC3E+06 

0.34K 

3.635 

0.119E-02 

0. 380E*04 

303.7 

FILM THICKNESS - 0.38100E-04 METERS 





X 

P 

MACH 

4FL/D 

OENSITV 

V 

T 

METERS 

NSMA 

NUMBER 


KG/M3 

M/SEC 

DEG K 

0 

0.446E+06 

0.917E-01 

80.25 

0.500E-02 

0,102E*04 

310,4 

0.127E-83 

0,424E+06 

0.963E-01 

72.47 

0.476E-02 

0.108E*04 

310.4 

0.254E-Q3 

0-402E+06 

0.102 

64.70 

0.451E-C2 

0.11 3E*04 

310.3 

0*381E-03 

0.378E+06 

0.108 

56.92 

0.425E-02 

0. 1216^04 

310.2 

0.508E-03 

0,353E+06 

0.116 

49.15 

0.397E-02 

0,IZ9E*04 

310.1 

0*635E-03 

0,326E+06 

0.125 

41.38 

0.366E-02 

0. 140E*04 

310.0 

0*762E~03 

0* 296E+06 

0.138 

33.60 

0.333E-02 

0. 154E+04 

309.7 

0.889E-03 

0«263E+06 

0.155 

25.33 

0.296E-02 

0, 173E*04 

309.4 

0. 102E- 02 

0«224E-i-06 

0.182 

18.05 

0.252E-02 

0. 203E*04 

308.9 

0, 114E- 02 

0.176E+06 

0.231 

10.28 

0.199E-02 

0.257E+04 

307.6 

0. 127E-02 

0.103E+06 

0.390 

2.5H6 

0.119E-02 

0.429E*04 

301.8 

FILM THICKNESS - 0*40640E-04 METERS 





X 

K 

M«CM 

4FL/C 

OENSITV 

V 

T 

METERS 

NSMA 

NUMBER 


KG/ M3 

M/SEC 

K 

0 

0.445E^06 

0-103 

62.45 

0.499E-02 

0.1156*04 

310.3 

O.127E-03 

0.4246*06 

0.108 

56. 38 

0.^76E-OZ 

0,121E*04 

310.2 

0.254E-03 

0.4026*06 

0. Il4 

50.31 

0.451E-02 

0,128E*04 

310.1 

O* 381E-03 

0.378E+06 

0.121 

44.24 

0.425E-02 

0.136E*04 

310.0 

0.508E-03 

0.353E+06 

0.130 

38.16 

0.397E-02 

0.145E+04 

309.9 

0*635E-03 

0.3266*06 

0.141 

32.09 

0.3676-02 

0. 157E*04 

309,7 

0.762E-03 

0.296E*06 

0.155 

26.02 

0.334E-02 

0.173E*04 

309.4 

0,889E-03 

0,2636*06 

0^174 

19,95 

0.2976-02 

0, 19'.E*04 

309.0 

0. 102E-02 

0.225E+06 

0.204 

13.88 

0.254E-02 

0,227E*04 

308.4 

0. 114E-02 

0.1776*06 

0.253 

7.307 

0.201E-02 

0.23TE*04 

306.8 

0.127E-02 

0.103E*06 

0.437 

1.736 

0.120E-02 

0.479E*04 

299.5 

FILM THICKNESS « 0.43X80E-04 METERS 





X 

p 

MACH 

4FL/D 

DENSITY 

V 

T 

METERS 

NSMA 

NUMBER 


KG/H3 

M/SEC 

DEG K 

0 

0.444E*06 

0.!15 

49.45 

0.499E-02 

0-129E*04 

310.1 

0*127E“03 

0.423E*06 

0.121 

44.62 

0.476E-02 

0.135E*04 

310.0 

0.254E-03 

0.401E*06 

6.128 

39.80 

0.451E-02 

0.142E*04 

309.9 

0.381E-03 

0.378E*06 

0.135 

34.98 

0.425E-02 

0.151E*04 

309.8 

0.508E-03 

0.353E+06 

0,145 

30«15 

0.397E-02 

0. I62E*04 

309.6 

0.635E-03 

0.326E+06 

0.157 

25.33 

0.367E-02 

0.175E*04 

309.4 

0.76ZE-03 

0,297E*06 

0.172 

20., 50 

0.335E-02 

0.192E*04 

309.1 

0. 8896-03 

0.264E+06 

0.194 

15.68 

0.298E-02 

0.216E*04 

306.6 

0. i02E-02 

0.226E+06 

0.226 

iO.85 

0.255E-02 

0.251E+04 

307.8 

0, U4E-02 

0.178E*06 

0.286 

6.028 

0.203E-02 

0.317E*04 

305.9 

0.127E-02 

0.1036*06 

0.485 

1.203 

0.121E-02 

0.529E*04 

297.0 


FILM THICKNESS * 0.45720E 

-04 METERS 





X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

METERS 

NSMA 

NUMBER 


K6/H3 

M/SEC 

DEC K 

0 

0. 4436+06 

0,l?8 

39,77 

0.498E-02 

0.14ZE+04 

309.9 

0.127E-03 

0.422E+06 

0, 13-i 

35.66 

0.475E-02 

0.149E+04 

309.3 

0.254E-03 

0.401E+06 

0.141 

31.98 

0.451E-02 

0, 157E+04 

309.7 

0. 381E-03 

0.378E+06 

0.150 

28,09 

0.425E-02 

0.167E+04 

309.5 

0.50BE-03 

0.353E+06 

0.160 

24.20 

0.398E-02 

O.lTBE+04 

309.3 

0.635E-05 

0.326E+06 

0.l7i 

20.30 

0.368E-02 

0.193E+04 

309. 1 

0.762E-03 

0.297E+06 

0.19C 

16,41 

0.335E-02 

0.211E+04 

308.7 

0*Se9E-05 

0.26SE+C6 

0.213 

12.51 

0.299E-02 

0.237E+04 

308« 1 

0.102E-OZ 

0.227E+06 

0.248 

8.62C 

0.257E-02 

0.276E+04 

307.1 

0.114E-02 

0.179E+06 

0.313 

4.726 

0.205E-02 

0.346E+04 

305.0 

0.127E-02 

0.103E+06 

0.533 

0.832 

0.122E-02 

0.579E+04 

294.2 



OUASI-CNE DIMENSIONAL COMPRESSIBLE FLOW SEAL PROGRAM 
INPUT DATA - 


SAMPLE PROBLEM 


U* CUSTOMARY UMTS 


Rlt INCHES 
3.2650 

?0,PSIA 

65.000 

MOLECULAR WEIGHT 
28.966 

F»K/RE**N 

R2, INCHES 
3.3150 

P3,PSIA 

15.000 

CP,BTU/LBM-D6G R 
0.240 

K(LAMINIR) 

24.00C 

FLOW LENGTH, INCHES 

0. 

P0/P3 

0. 

GAMMA 

1.400 

NILAMINARI 

1.00 

FLOW WIDTH, INCHES 

0. 

TO, DEG F 
100.0 

VISCOSITY. 13-SEC/ I N2 
0 

UPPER LIMIT RE (LAMINARJ 
2300. C 


LOSS COEF. 
I. 00 

SPEED, RPM 

0. 

K( TURBULENT) 
0.7900E>01 



V.FT/SEC 

200.00 

NI TURBULENT I 
0.2500 




LOWER LIMIT RE (TURBULENT! 
3000. C 

r DATA - 




Rl, INCHES 
3.2650 

PO.PSIA 

65.000 

MOLECULAR WEIGHT 
28.966 

CALCULATE 

R2, INCHES 
3.3150 

P3,PSIA 

15,000 

CP,BTU/L6M-DEG R 
0.240 

POWER 

DIMENSIONLESS CJANTITIES 

FLOW LENGTH, INCHES 
D.0500 

PD/r3 
4. 3 " 5 

GAMMA 

1.400 

PRESSURE DISTRieUTiONS 

FLOW WIDTH, INCHES 
20.6717 

TO, DEC F 
100.000 

VI SCOSITY, LB-SEC/ IN2 
0.2754E-08 

PLOT 

AREA, IN2 
1.0336 

GAS CONSTANT, IB-FT/L9M'-R 
53.35221 

SPEED, RPM 
6966.1 

POWER XC(BAR) 

FORCE PRESSURE 


LOSS COEF. 
loOO 

V, FT/SEC 
200.00 

TEMPERATURE DENSITY 


MACH NUMBER AFL /D 


• 

• * “ CHOKING FILM THICKNESS 

• ♦ - TRANSITION REGION 
» / - TURBULENT FLOW 


QUASI-CNE dimensional COMPRESSIBLE FLOW SEAL PROGRAM 


SAMPLE PROBLEM - U. 


FILM THICKNESS 
I NCHES 

• O.51595E-03 

0. lOOOOE-03 
0,20000E-03 
0. 300DOE-03 
0.40000E-03 
0.50000E-03 
0.60000E-03 
0,70000E-03 
0.80000E-03 
0,9000DE-03 
0* lOOOOE-02 
+ 0,11000E-02 

♦ 0.12000E-02 

♦ 0.13000E-02 

/ 0.1A030E-02 

/ 0.15000E-02 

/ 0.16000E-02 

/ 0. 17000E-02 

/ 0.18000E-02 


MIDOT J 


a 

KNUDSEN 

LAMBDA 

STIFFNESS 

FORCE 

LB/MIN 


SCFM 

NUMBER 

INCHES 

LB/IN 

LB 

0.408 


5.342 

0.473E-02 

0.244E-05 

-0.258E+04 

31.84 

O.AlflE 

02 

0.547E-0I 

0.270E-01 

0.270E-05 

-30.00 

31.23 

0-330E-D1 

0.432 

0.134E-01 

0.269E-05 

-679. D 

31.26 

0. 106 


1.386 

O.079E-O2 

0.264E-05 

-0.157E»04 

31.37 

0.226 


2*961 

0.638E-02 

0.255E-05 

-0.222E>04 

31.56 

0.361 


4,990 

0.491E-02 

0.246E-05 

-0.254E+04 

31.80 

0.557 


7.286 

0.345E-02 

0.207E-05 

-383.0 

31.97 

0.745 


9,742 

0.257E-02 

0.180E-05 

0.154E^04 

31.91 

0.939 


12.28 

0.203E-02 

0.162E-05 

0.272E+04 

31.69 

1. 136 


I4.B6 

0.167E-02 

0. 150E-05 

0.337E*^04 

31.38 

1. 335 


17r46 

0.l4lt-02 

O.141E-05 

0.369E+04 

31.02 

1.520 


19.68 

0.123E-02 

0.136E-05 

0 

30.74 

1*680 


21.97 

0.112E-02 

0.134E-05 

0 

30.61 

1.841 


24,09 

0.102E-02 

0.132E-05 

0 

30.49 

2.011 


26.31 

0.928E-0J 

0. 130E-05 

0.172E+04 

30. 33 

2,185 


28.56 

0.B53E-03 

0. 128E-D5 

0. 166E<’04 

30.16 

2.359 


30.8 7 

0.788E-03 

0.126E-05 

0.160E+04 

30.00 

2.535 


33.16 

0.733E-03 

0.125E-05 

0,154E>04 

29.34 

2.710 


35.46 

0.684E-03 

0.123E-05 

0.149E*04 

29,69 


CUSTCHARY 

UNITS 

XC 

FCRCE 

INCHES 

EAR 

O.190E-QI 

C.616 

0.186E-01 

r .(>04 

0.187E-01 

C .60 5 

0.189E-01 

0 fSC7 

0.192E-01 

C .61 1 

0.197E-CI 

C .^15 

0.203E-01 

0 ol9 

0.207E-01 

C.617 

0.211E-QI 

0 613 

0.214E-01 

0 .oC7 

0.217E-QI 

0-6CC 

O.210E-OI 

: ,595 

0.219E-01 

^ .^92 

0.220E-01 

: .590 

0.220E-0I 

0 .58 7 

0.221E- 01 

f .584 

0.222E-01 

/ ,5 8C 

0.223E-01 

C.577 

0.223E-01 

C.575 


XC 
BAR 
0.397 
0,373 
0.37A 
0.377 
0.305 
0.395 
0.A06 
0.A15 
0.42Z 
O.A26 
0.A33 
0.A37 
0.A38 
0.43 9 
0.441 
a. 442 
e .444 
Oc i-^5 

0.4^ 7 
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FILM THICKNESS 

FRICTION 

REtP) 

RE(R) 

PI«I 

U*l 

POWER 

HEAT LOSS 

DEL ITI 

TORQUE 


INCHES 

FACTOR 



PSIA 

INCHES 

H*P* 

BTU/NIN 

0E6 F 

L8/FT 

• 

0.51593E-03 

0«38AE>01 

625*8 

322*2 

15*00 

0*500E*01 

0.473E-02 

0*201 

2.C48 

0*224E<>01 


O«I0000E-03 

3*782 

6*345 

59.41 

0*793 

0.S27E-01 

0«248E*0I 

1*054 

0*105E^04 

0*118 


0.20000E-03 

O.A79 

50.12 

119*1 

3*131 

0.522E-01 

0.124E-01 

0*526 

66.61 

0*588E-01 


0,30000E-03 

0.1A9 

161.0 

180.3 

6*696 

0.SI3E-0I 

0.825E>02 

0*350 

13.76 

0*391 E*01 


0.40000E-03 

0.695E-01 

345.1 

244.3 

10*73 

0.504E-01 

0*615E>02 

0.261 

4*805 

0*291E-01 


0.50000E-03 

0.411E-01 

584.1 

311.3 

14.46 

0*500E-01 

0.489E-02 

0*207 

2*<65 

0.232E-01 


0.60000E-03 

0.260E-01 

857.4 

380.0 



0.404E-02 

0*171 

1*283 

0*191E>01 


0.70000E-03 

0.208E>01 

0.115E404 

448.9 



0*344E-02 

0*146 

0.816 

Q.163E-01 


0«80000E*03 

0«164E>01 

0.146E«^04 

517.4 



6*299E-02 

0*127 

0*563 

0.142E-01 


0«90000E’03 

0.135E-01 

0.178E«04 

585.4 



0.264E-02 

0*112 

0*411 

0.125E-01 


O.lOOOOE-02 

0.114E-01 

0.210E«04 

652. B 



0.236E-02 

0*100 

0*313 

G*112E-01 

♦ 

0«11000E>02 

0.105E-01 

0*240E^04 

719.2 







♦ 

0.I2000E-02 

0.106E-01 

0*265E«04 

785.0 








0«1300DE*02 

0.107E*01 

0.291E«^04 

850.7 







/ 

0.14000E-02 

O.IOSE-Ol 

0*319E^04 

916*5 







/ 

0.15000E-02 

0«103E>01 

0.347E^04 

982.1 







/ 

0.I6000E>02 

O.lOlE-01 

0.375E^04 

0*105E^04 







/ 

0.17000E-02 

0.991E-02 

C.404E4>04 

0*111E^04 







/ 

0,16000E-02 

0.974E-02 

0*433E^04 

0.I18E404 








* FILH THICHNESS • 0.51595E-03 INCHES 


X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


LB-SEC2/FT4 

FT/SEC 

DEG 1 

0 

61.94 

0.263 

7.433 

0.941E>02 

303.5 

92*33 

0.500E-02 

59.41 

0.275 

6.690 

0.903E-02 

316*1 

91.68 

O.IODE-Dl 

56.74 

0.287 

5,947 

0.864E-02 

330.5 

90,91 

0, 150E-01 

53.91 

0.3G2 

5.203 

0.822E-02 

347.2 

89*96 

O.2OOE-0I 

50.89 

0.320 

4.460 

0.77BE>02 

367.1 

88*79 

0.250E-01 

47. 62 

0.341 

3.717 

0.730E-02 

391.2 

87.26 

0.300E-01 

44.03 

0.368 

2.973 

0.677E-02 

421.5 

85.21 

0.350E-01 

39.98 

0.404 

2.230 

0.619E>02 

461.6 

82*26 

0.400E-01 

35.24 

0*457 

1.487 

0.550E-02 

519.2 

77.56 

0.450E-01 

29.16 

0.547 

0.743 

0.463E-02 

616.7 

68.35 

0.500E-01 

15.00 

1.000 

0 

0.270E-02 

0.106E404 

6,667 


FILH THICKNESS - O.lOOOOE-03 INCHES 


X 

P 

MACH 

4FL/0 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


LB-SEC2/FT4 

FT/SEC 

DEG F 

0 

64.99 

0.134E-01 

0.399E+04 

0.974E-02 

15,51 

99,98 

0.500E-02 

61.84 

0.140E-01 

0.36IE4-04 

0.926E-02 

16.30 

99.98 

O.lOOE'Ol 

58.52 

0.14BE-01 

0.323E4^04 

0.877E-02 

17.22 

99.98 

O.lSOE-01 

55.00 

0.1S8E-01 

0,^B6E+04 

0.824E-02 

18,32 

99,97 

0.200E-01 

51.23 

0.170E-01 

0.248E^04 

0.768E-02 

19,67 

99.97 

0.250E-01 

47.17 

0.I84E-01 

0,210E^04 

0.707E-02 

21.36 

99.96 

0.300E-01 

42,72 

0.203E-01 

0.172E«04 

0.640E-02 

23.59 

99,95 

0.350E-01 

37.76 

0.230E-OI 

0.134E^04 

0.566E-02 

26.69 

99.94 

0.400E-01 

32.03 

0.27IE-01 

964.1 

0.480E-02 

31.46 

99.92 

0.450E-01 

25.01 

0.347B-01 

585.9 

0.375E-02 

40.28 

99.86 

0.500E-01 

15.00 

0.579E-01 

207.7 

0.225E-O2 

67.14 

99*62 


FILM THICKNESS - 0.20000E-03 INCHES 


X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


LB-SEC2/FT4 

FT/SEC 

DEG 1 

0 

64.87 

0.529E-01 

250.1 

0.972E-02 

61.29 

99.69 

0.500E-02 

61.75 

0.555E-01 

226,2 

0.926E>02 

64.39 

99,65 

0. lOOE-OI 

58.46 

0.586E-01 

202.2 

0.876E-02 

68.01 

99.62 

0.150E>01 

54.96 

0.624E-01 

178.3 

0.824E-02 

72.33 

99.56 

0*200E-01 

51.23 

0.669E-01 

154,4 

0.768E-02 

77.58 

99,50 

0.250E-01 

47.21 

0.726E-01 

130.4 

0.708E-02 

84.19 

99.41 

0.300E-01 

42.79 

0.801E-01 

106.5 

0.642E>02 

92.85 

99.28 

0.350E-01 

37.86 

0,905E-01 

82.53 

0.568E-02 

104.9 

99.08 

0*400E-0I 

32.17 

0.106 

58,59 

0.483E-02 

123.4 

98,73 

0.450E-01 

25.17 

0.136 

34.64 

0.378E-02 

157.5 

97.94 

0.500E-Oi 

15,00 

0.227 

10.70 

0.227E-02 

262.5 

94,26 


FILM THICKNESS « 0.30000E-03 INCHES 


X 

P 

MACH 

4FL/0 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


LB-SEC2/FT4 

FT/SEC 

DEG 1 

0 

64,41 

0.114 

50.94 

0. 967E-02 

131.8 

98.56 

0.500E-02 

61.39 

0.119 

45.97 

0.922E-02 

136.2 

98,41 

O.IOOE-Ol 

58.21 

0.126 

41.00 

0.B75E-02 

145.7 

98.23 

0.150E-01 

54.83 

0.134 

36.04 

0.824E-02 

156,7 

98.01 

0.200E-01 

51.22 

0.143 

31,07 

0.770E-02 

165.5 

97.72 

0.250E-01 

47.32 

0.155 

26*10 

0.712E-02 

179,0 

97.33 

0.300E-01 

43.04 

0.170 

21.13 

0.648E-02 

196.6 

96.78 

0.350E-01 

38.25 

0.191 

16.17 

0.577E-02 

220.6 

95.94 

0.400E-01 

32.70 

0.223 

11.20 

0. 495^-02 

257.7 

94.47 

0.450E-01 

25.79 

0.282 

6.23 b 

0.392E-Q2 

3». 8 

91.22 

' 500E-01 

15.00 

0.478 

1.263 

0.235E-02 

542.5 

75.51 
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FILM THICKNESS - 0.40000E-03 INCHES 


X 

P 

MACH 

INCHES 

PSIA 

NUMBER 

0 

63.47 

0.184 

0.500E-02 

60.65 

0.193 

O.IOOE-Ol 

57.67 

0,203 

0.150E-01 

54.51 

0.215 

0.200E-01 

51.13 

0.229 

0.250E-01 

47.48 

0.246 

0.300E-01 

43.47 

0.268 

0.350E-01 

38.97 

0.299 

0.400E-01 

33.71 

0.344 

0.450E-01 

27,05 

0.427 

0.500E-01 

15.00 

0.743 

FILM THICKNESS - 0. 

50000E-03 INCHES 

X 

P 

MACH 

INCHES 

PSIA 

NUMBER 

0 

62.17 

0.253 

0.500E-02 

59,59 

0.264 

O.IOOE-Ol 

56.88 

0.276 

0.150E-01 

54.00 

0,291 

0.200E-01 

50,93 

0.308 

0.250E-01 

47.60 

0.329 

O.3COE-01 

43,95 

0.356 

0.3506-01 

39,84 

0.392 

0.400E-01 

35.02 

0.444 

0.450E-01 

28.84 

0.534 

0.500E-01 

15.00 

0.969 


4FL/0 

17.52 

15.78 

14.05 

12.31 

10.57 

8.B30 

7.091 

5.352 

3.614 

1.875 

0.136 


4FL/D 

8.Z16 
7.397 
6.575 
5.753 
4.932 
4.110 
3. 288 
2.466 
1.645 
0.823 
0.121E-02 


DENSITY V T 

L8-SEC2/FT4 FT/SEC DEC F 

0.957E-02 213.3 96.21 

0.9I5E-02 223.1 95.86 

0.871E-02 234.5 95.43 

0.824E-02 247.8 94.89 

0.774E-02 263.8 94.21 

0.720E-02 283.7 93.30 

0.661E-02 309.1 92.05 

0.594E-02 343.7 90.17 

0.517E-02 395.0 87.02 

0.420E-02 486.2 80.33 

0.250E-02 818.3 44.27 


DENSITY V T 

LB-SEC2/FT4 FT/ SEC DEG F 

0.943E-02 291.9 92.91 

0.905E-02 304.1 92.30 

0.865E-0Z 318.2 91.57 

0.823E-02 334.6 90.68 

0.777E-02 354,1 89.57 

0.729E-02 377,8 88.12 

0.675E-02 407.8 86.16 

0,615E-‘02 447,5 83,33 

0.545E-02 504.8 78.79 

0.457E-02 602.7 69,77 

0,267B-02 0,103E^04 11.46 


FILM THICKNESS - 0.60000E-03 INCHES 


X 

P 

INCHES 

PSIA 

0 

60.69 

0.500E-02 

56.38 

O.IOOE-Ol 

55,94 

0.150E-01 

53.35 

0.200E-01 

50.58 

0.250E-01 

47.59 

0.300E-01 

44.31 

0.350E-01 

40*61 

0.400E-0I 

36.26 

0.450E-01 

30.6B 

0,500E-DI 

17,59 


MACH 

NUMBER 

4FL/D 

0.314 

4,666 

0,327 

4.199 

0,341 

3.732 

0.3^7 

3.266 

0.376 

2.799 

0.399 

2,333 

0.427 

1,866 

0.465 

1.400 

0.518 

0.933 

0.606 

0.467 

1.000 

0 


DENSITY V T 

LB-SEC2/FT4 FT/SEC DEC F 

0.927E-02 361.2 89.14 

0.893E-02 375.0 88.30 

0.857E-02 390.6 87.30 

0.820E-02 408,7 86.10 

0.779E-02 429.9 84.62 

0.736E-02 455.3 82,75 

0.688E-02 486.9 80.27 

0.635E-02 527.8 76,82 

0.572E-02 585.2 71.50 

0.493E-02 678.8 61.65 


0.316E-02 0, 106E>04 6.667 


FILM THICKNESS - 0.70000E-03 INCHES 


X 

P 

HACH 

4FL/D 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


L8-SEC2/FT4 

FT/SEC 

DEG 1 

0 

59.19 

0.368 

2.975 

0.911E-02 

421,5 

85.22 

0. 500E-02 

57, 10 

0.381 

2.677 

0.880E-02 

436.1 

04. 17 

O.IOOE-Ol 

54.90 

0.396 

2.380 

0.84BE-02 

452,5 

82*96 

0. 150E-01 

52.57 

0.413 

2,082 

0.8I4E-02 

471.3 

81.51 

0.200E-01 

50.08 

0.433 

1.785 

0.778E-02 

493.2 

79.76 

0.250E-01 

47.38 

0.457 

1.487 

0,7396-02 

519.2 

77. 57 

0.300E-01 

44.41 

0.486 

1.190 

0.697E-02 

550.9 

74.74 

0.3506-01 

41.07 

0,524 

0,892 

0.649E-02 

591.4 

70,89 

0.400E-01 

37.15 

0,576 

0.595 

0.593E-02 

646.9 

65.17 

O.450E-01 

32.09 

0.660 

0.297 

0.523E-02 

734.5 

55. 10 

0.500E-01 

20.16 

1.000 

0 

0.362E-02 

0.106E«04 

6.667 


FILM THICKNESS - 0.80000E-03 INCHES 


X 

P 

INCHES 

PSIA 

0 

57.74 

0.500E-OZ 

55,05 

O.IOOE-Ol 

53.85 

0.150E-01 

51.74 

0.200E-01 

49.48 

0.250E-01 

47.04 

0.300E-01 

44, 35 

0.350E-01 

41.32 

0.400E-01 

37.75 

0.450E-01 

33.15 

0. 500E-01 

22.24 


MACH 

4FL/0 

NUMBER 


0,415 

2.053 

0.429 

1.848 

0.444 

1.642 

0,461 

i.437 

0.481 

1.232 

0.505 

1.027 

0.534 

0.821 

0,571 

0.616 

0.622 

0.411 

0.701 

0.205 

l.OOO 

0 


DENSITY 

LB-SEC2/FT4 

0.895E-02 

0.867E-02 

0.839E-02 

0.808E-02 

0.776E-02 

0.741E-02 

0.702E-02 

0.659E-02 

0.609E-02 

0.545E-02 

0.400E-02 


V 

FT/SEC 

473.3 

488.3 

505.0 

524.0 

546.0 
571.8 

603.0 

642.3 

695.1 
776.5 
0.106E^04 


T 

OEG F 
81.36 
80.16 
78.77 
77.14 
75.19 

72.79 
69.74 
65.67 

59.79 
49.83 
6.667 


FILM THICKNESS - 0.90000E-03 INCHES 


X 

p 

MACH 

4FL/0 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


LB-SEC2/FT4 

FT/SEC 

DEG 1 

0 

56.37 

0.456 

1.500 

0.880E-02 

517.9 

77 * 67 

0.500E-02 

54.65 

0.469 

1.350 

0.855E-02 

533.0 

76. 36 

0. lOOE-Ol 

52.84 

0.465 

1.200 

0.829E-02 

549.7 

74.85 

0.150E-01 

50,91 

0.502 

1.050 

0.801E-02 

568.6 

73.09 

0.200E-01 

48.86 

0.522 

0.900 

0.7T2E-02 

590.2 

71.01 

0.250E-01 

46.63 

0.546 

0.750 

0.740E-02 

615.5 

68.47 

0. 300E-01 

44.17 

0.575 

0.600 

0.705E-02 

645.8 

65.29 

0.350E-01 

41.41 

0.611 

0,450 

0.667E-02 

683.5 

61. 12 

0.400E-01 

38.15 

0.659 

0.300 

0.621E-02 

733.5 

55.23 

0.450E-01 

33.93 

0,734 

0.150 

0.563E-02 

809.0 

45.53 

0.500E-01 

23.93 

1.000 

0 

0.430E-02 

0.106E»04 

6. 667 



FILM THICKNESS « 0.100006-02 INCHES 


X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


L6-SEC2/FT4 

FT/SEC 

OEG F 

0 

55,12 

0.491 

1.143 

0.865E-02 

556.5 

74.22 

0. 5006-02 

53,54 

0.505 

1.029 

0,0436-02 

571.4 

72.82 

O.IOOE-Ol 

51.68 

0.520 

0.915 

0.619E-02 

58B.0 

71.23 

0. 150E-01 

50.12 

0.538 

O.BOO 

0.794E-02 

606,5 

69.38 

0.200E-01 

48.23 

0.550 

0.686 

0.7676-02 

62 7,7 

67.21 

0.250E-01 

46.18 

0.581 

0.572 

0.7396-02 

652,2 

64.60 

0.300E-01 

43.93 

0,609 

0.457 

0,7076-02 

6B1.4 

61.36 

0.350E-01 

41,39 

0,643 

0.343 

0,6716-02 

717.4 

57.17 

O.^OOE-Ol 

36.40 

0.690 

0.229 

0,6306-02 

764.6 

51.34 

0.450E-01 

34.53 

0.760 

0.114 

0.577E-02 

834,9 

41.99 

0.500E-01 

25,30 

l.OCO 

0 

0.4556-02 

0,106Et04 

6.667 

FILM THICKNESS » 0, 

,110006-02 INCHES 




X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


LB-SEC2/FT4 

FT/SEC 

□ EG F 

0 

54,23 

0.515 

0.950 

0.855E-02 

582,7 

71.75 

0.500E-02 

52.75 

0,529 

0.855 

0.834E-02 

597.4 

70.30 

O.IOOE-Ol 

51.19 

0.544 

0.760 

0.812E-02 

613.7 

68,66 

0.150E-01 

49,53 

0.562 

0.665 

0. 789E-02 

631.9 

66.77 

0,2006-01 

47.76 

0.581 

0.570 

0.764E-02 

652.6 

64.56 

0.250E-01 

45.84 

0.604 

0.475 

0.737E-02 

676.5 

61.91 

0,300E-0l 

43.73 

0.631 

0.380 

0.707E-02 

704.8 

58,66 

0,350E-W 

41,34 

0.665 

0.285 

0.674E-02 

739.5 

54.49 

0,400E-01 

38.52 

0.710 

0.190 

0.635E-02 

784,7 

48.75 

0.450E-01 

34.88 

0,777 

0.950E-01 

0.586E-02 

B51.3 

39.69 

0.500E-01 

26.16 

1.000 

0 

0.471E-02 

0. 1066+04 

6.667 

FILM THICKNESS - 0. 

12000E-02 INCHES 




X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


LB-SEC2/FT4 

FT/SEC 

OEG F 

0 

53,86 

0.525 

0.881 

0.851E-02 

593,3 

70.71 

0.500E-02 

52,42 

0.539 

0.793 

0.831E-02 

607,9 

69.25 

O.IOOE-Ol 

50,90 

0.554 

0.705 

0.809E-02 

624.1 

67.59 

O.I50E-01 

49. 29 

0.571 

0.616 

0,787E-02 

642.1 

65,69 

0.200E-01 

47,56 

0.591 

0.528 

0.762E-02 

662.6 

63.46 

0,2506-01 

45.69 

0.613 

0.440 

0.736E-02 

666.2 

60,81 

0.3006-01 

43,63 

0.640 

0.352 

0.707E-02 

714.1 

57.56 

0.350E-01 

41.31 

0.674 

0.264 

0.675E-02 

748.3 

53.40 

0.400E-01 

38.56 

0.718 

0.176 

0.637E-02 

792.7 

47.71 

0.450E-01 

35.01 

0.783 

0.B81E-01 

0.589E-02 

857.8 

38.77 

0,500E-01 

26.53 

l.OCO 

0 

0,477E-02 

0.106E+04 

6.667 

FILM THICKNESS > 0, 

13000E-02 INCHES 





X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


LB-SEC2/FT4 

FT/SEC 

OEG F 

0 

53.51 

0.535 

0.820 

0.647E-02 

603.1 

69,73 

0. 500E-02 

52,11 

0.548 

0.738 

0.BZ7E-O2 

617.6 

68.25 

0,1006-01 

50,63 

0,563 

0.656 

0.807E-02 

633.7 

66.58 

0. 1506-01 

49.05 

0.580 

0.574 

0.784E-02 

651.6 

64.66 

0.2006-01 

47.37 

0.600 

0.492 

0.761E-02 

671,9 

62.43 

0.250E-01 

45.55 

0,622 

0.410 

0,735E-02 

695,2 

59,78 

0.300E-01 

43,54 

0,649 

0.328 

0.707E-02 

722.7 

56.53 

0,3506-01 

41.27 

0,682 

0.246 

0.676E-02 

756.4 

52,39 

0.400E-01 

38.59 

0.725 

0.164 

0.639E-02 

800,0 

46.74 

0.450E-01 

35,12 

0,790 

0.820E-01 

0.592E-02 

863.6 

37.93 

0.500E-01 

26,85 

1.000 

0 

0.483E-02 

O.I06E+04 

6,667 

FILM THICKNESS » 0. ! 

14000E-02 INCHES 





X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


LB-SEC2/FT4 

pT/SEC 

DEG F 

0 

53,07 

0,546 

0.751 

0.8426-02 

615,3 

68,49 

0,500E-02 

51.71 

0.560 

0.676 

0.8236-02 

629.7 

67,00 

O.lOOE-01 

50.28 

0,575 

0.601 

0.8036-02 

645.6 

65,31 

0. 150E-01 

48.76 

0,591 

0.526 

0.7826-02 

663.3 

63.38 

0.200E-01 

47,13 

0,611 

0.451 

0,7596-02 

683,3 

61,14 

0.250E-01 

45,36 

0,633 

0.375 

0.7346-02 

706.3 

58,48 

0.300E-01 

43.42 

0,659 

0.300 

0. 7076-02 

733.3 

55,24 

0,3506-01 

41.22 

0.691 

0.225 

0.676E-02 

766.3 

51.13 

0,4006-01 

38.62 

0,734 

0.150 

0,6416-02 

808.9 

45,55 

0.450E-01 

35.26 

0,797 

0,7516-01 

0.595E-02 

870,8 

36.89 

0,5006-01 

27.23 

l.OOO 

0 

0.490E-02 

0,1066404 

6,667 

FILM THICKNESS - 0.15000E-02 INCHES 





X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


LB-SEC2/FT4 

FT/SEC 

OEG F 

0 

52.63 

0,558 

0.686 

0.837E-02 

627.7 

67,21 

0.500E-02 

51.31 

0,571 

0.616 

0.819E-02 

641,9 

65,71 

O.IOOE-Ol 

49,93 

0,586 

0.549 

0.799E-02 

657.6 

64,01 

0.150E-01 

46,45 

0,603 

0.480 

0,7796-02 

675.1 

62.07 

0.2006-01 

46,88 

0,622 

0.412 

0.757E-02 

694.8 

59,83 

0.2506-01 

45,17 

0,643 

0.343 

0.733E-02 

717.4 

57,17 

0.300E-01 

43.29 

0.669 

0.274 

0.707E-02 

743.9 

53,94 

0.350E-01 

41.16 

0.701 

0.206 

0.677E-02 

776.2 

49,86 

0.406e-01 

38,64 

0.743 

0.137 

0.643E-02 

817.7 

44,35 

0,450E-01 

35,39 

0.804 

0.686E-01 

0.599E-02 

877.9 

35,86 

0.500E-01 

27.61 

1.000 

0 

0.496E-02 

0.1066404 

6.667 



FILM THICiCNESS « 0,16000E-02 INCHES 


X 

p 

MACH 

4FL/D 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


LB-SEC2/FT4 

FT/SEC 

DEG F 

0 

52*21 

0.568 

0*631 

0.833E-02 

639.1 

66*01 

0.500E-0Z 

50*93 

0.502 

0*568 

0.815E-02 

653.1 

64.50 

0, lOOE-Ol 

49,59 

0.597 

0*505 

0.796E-02 

668 • 7 

62.79 

0.150E-01 

48*17 

0.613 

0.442 

0*7766-02 

685.9 

60*85 

0*200E-01 

46*64 

0.632 

0,378 

0,7546-02 

705.3 

58*60 

0.250E-01 

44,98 

0.653 

0.315 

0. 7316-02 

727.5 

55.95 

0.300E-01 

43, 16 

0.679 

0.252 

0.706E-02 

753.6 

52* 74 

0* 350E-01 

41,09 

0.710 

0.189 

0.6786-02 

765.2 

48*69 

0.400E-01 

38.66 

0,751 

0.126 

0.644E-02 

825.7 

43*26 

0.450E-01 

35.50 

0.811 

0.631E-01 

0.602E-02 

884. 3 

34*92 

0. 500E-01 

27.95 

l.OCO 

0 

0.5026-02 

0.106E+04 

6.667 

FILM THICKNESS » 0, 

17000E-02 INCHES 





X 

P 

MACH 

4FL/D 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


LB-S6C2/FT4 

FT/SEC 

DEG F 

0 

51.82 

0.578 

0.583 

0.828E-02 

649.6 

64. 88 

0.500E-02 

50,50 

0.592 

0.525 

0,8116-02 

663.6 

63.35 

O.IOOE-Ol 

49,28 

0.6C6 

0.466 

0.792E-02 

678*9 

6 1 * 64 

0. 150E-01 

47,89 

0.623 

0*408 

0.773E-02 

695.9 

59.70 

0. 200E-01 

46.41 

0.641 

6.350 

0.752E-02 

715.0 

57*45 

0.250E-01 

44.80 

0,663 

0*291 

0,730E-02 

736.9 

54* 81 

0.300E-01 

43.03 

0.6B8 

0.233 

0*7066-02 

762. 5 

51*62 

0.350E-01 

41.03 

0,718 

0.175 

0.678E-02 

793.4 

47*61 

0*400E-01 

38.66 

0,758 

0.117 

0.646E-02 

833*1 

42*24 

0*4506-01 

35.60 

0.817 

0.583E-01 

0.604E-02 

890.1 

34*06 

0.500E-01 

28,26 

1,000 

0 

0.508E-02 

0. 106E+04 

6* 667 


FILM THICKNESS « 0. 

18000E-OZ INCHES 





X 

p 

MACH 

4FL/0 

DENSITY 

V 

T 

INCHES 

PSIA 

NUMBER 


LB-SEC2/FT4 

FT/SEC 

DEG 1 

Q 

51*45 

0.588 

0.541 

0.824E-02 

659,5 

63*80 

0* 5C0E-02 

50*24 

0.601 

0,487 

0.8076-02 

673.3 

62*27 

0. lOOE-Ol 

48*96 

0.615 

0.433 

0.7896-02 

688.4 

60* 56 

0.150E-01 

47.63 

0.632 

0.379 

0.7716-02 

705.2 

50*61 

0*2006-01 

46*19 

0.6 50 

0.325 

0.7506-02 

724.0 

56*37 

0*2506-01 

44.63 

0.671 

0.271 

0. 7296-02 

745.6 

53*74 

0* 3006-01 

42.91 

0.696 

0.216 

0,7056-02 

770.7 

50, 57 

0* 3506-01 

40.96 

0.726 

0.162 

0. 6706-02 

801.1 

46* 60 

0.4006-01 

38.66 

0.765 

0.108 

0.6476-02 

839.8 

41.30 

0.450E-01 

35.68 

0.822 

0. 5416-01 

0.6076-02 

895.5 

33*27 

0.500E-01 

28.54 

1,000 

0 

0.5136-02 

0, 106E<»'04 

6 • 667 
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TABLE 1. - PARAMETERS CALCULATED BY QUASC 


Parameter 

Formula^ 

Units 



SI 

U.S. 

Gas constant 

^ Universal gas constant 
Molecular weight 

J/kg-K 

Ibf-in. 

lbm-“R 

Seal surface area 

A = !T (r| - Rj) 

2 

m 

in. 2 

Flow length 

AR = Rg - Rj 

m 

in. 

Flow width (mean) 


m 

in. 

Mass flow rate 

M = Whpu 

kg/ sec 

Ibm/min 

Volume flow rate 

Q = 

Re(p) = pu2h/p 

scms 

scfm 

Reynolds number due 
to leakage flow 

Reynolds number due 



^Re(r) = pVh/ p 



to rotational flow 




Knudsen number 

2.96 M^^^ 
max 



Kn - 

1 Re(p) 



Mean free path 

X = Kn X h 

m 

in. 

Viscosity of air 
(Sutherland’s law) 

C2T^* ^ 
'"a'’’ " T +Cg 

2 

N- sec/m 

2 

Ibf- sec/in. 

Power 

(lAV^/h 

W 

bp 


c^ X Power 


On 

Apparent temperature 

T = -^-r 

MCp 

K 

R 

rise due to power 
dissipation 




Shear heat 

H = Cg X Power 

W 

Btu/min 

Torque 

Power/ Speed 

N-m 

ft-lb 

Sealing dam force 

/'AR 

F = W / (P - P,)dx 

N 

Ibf 

Center of pressure 

X =W/'"''(P-Pg)xdX 

m 

in. 

Axial film stiffness 

S = -dF/dh 

N/m 

Ibf/in. 


^Constants in equations are as follows (for SI and U. S. units, respectively). 
Cj = 5. 051554x10-3 (13, 083); Cg = 1- 4591x10'® (1. 57639x10-^°); 

Cg = no. 33333 (198.6); = 1. 0 (42. 42); and Cg = 1. 0 (42. 42). 

'’where p is density at midseal, p viscosity at midseal, and V is mean 
rotational velocity. 




TABLE II. - VARIABLES IN NAME LIST/ SDATA/ 


Variable 

Description 

Units 

name 


SI 

U. S. 

RlIN 

Inner radius of seal 

m 

in. 

R2IN 

Outer radius of seal 

m 

in. 

RDIFIN 

Flow length 

m 

in. 

WIDTH 

MOLWT 

Mean flow width 
Molecular weight 

m 

in. 

CP 

Specific heat 

J/(kg)(K) 

Btu/(lbm)(°R) 

MU 

GAMMA 

Reservoir viscosity. The program will calculate MU 
for air but not for other gases. 

Ratio of specific heats 

N-sec,/m^ 

Ibf- sec/in. ^ 

SPEED 

Rotational velocity 

rps 

rpm 

CAPV 

Seal face speed 

m/sec 

ft/sec 

XLAM 

Exponent in friction factor * Reynolds number relation 
for laminar flow (eq. (10)) 

— 

— 

CONLAM 

Constant in friction factor - Reynolds number relation 
for laminar flow 

— 

— 

XTURB 

Exponent in friction factor - Reynolds number relation 
for turbulent flow 

— 

— 

CONTRB 

Constant in friction factor - Reynolds number relation 
for turbulent flow 

— 

--- 

RE LAM 

Maximum Reynolds number for laminar flow 

— 

— 

RETURB 

Minimum Reynolds number for turbulent flow 

— 

— 

PWRSKP 

Logical variable. It it is set to . TRUE. , calculations 
involving power are omitted. 

— 

— 

NRMSKP 

Logical variable. If it is set to . TRUE. , normalized 
values of F and will be omitted. 

— 

— 

PRSSKP 

Logical variable. If it is set to . TRUE. , printout 
of distributions across face of seal of P, T, p, u, 
M, and 4f(L - x)/D will be omitted. 



PLTSKP 

Array of logical variables. If any element is set to 
. TRUE . , the corresponding plot will be omitted: 
PLTSKP (1) applies to plot of power against h. 
PLTSKP (2) applies to plot of x^ against h. 
PLTSKP (3) applies to plot of F against h. 

PLTSKP (4) applies to plot of P against x, 

PLTSKP (5) applies to plot of T against x. 

PLTSKP (6) applies to plot of p against x. 

PLTSKP (7) applies to plot of M against x. 
PLTSKl' (8) applies to plot of 4f(L - x)/D 
against x. 



SKPH 

Logical variable. If it is set to . TRUE . , no film thick- 
ness data will be read. 

— 

— 

NOST 

Logical variable. If it is set to .TRUE input, inter- 
nal calculations, and output will be in U. S units. If 
it is set to , FALSE . , input, internal calculations, 
and output will be in SI units. 






TABLE m. - VARIABLES IN NAME LIST/PD AT A/ 


Variable 

Description 

Units 

name 


SI 

U. S. 

POIN 

Sealed gas pressure (upstream reservoir pres- 
sure), Pq 

Ambient pressure (downstream reservoir pres- 
sure), Pg 

N/ni^ 

CM 

c 

P3IN 

N/m^ 

Ibf/in. ^ 

PRIN 

Sealed gas/ambient pressure ratio, Pq/P^ 


--- 

TO IN 

Sealed gas temperature (upstream reservoir 
temperature), Tq 

K 


LOSS 

Entrance velocity loss coefficient 



INCODE 

Input code. For running many cases with one 
loading of the program, the input code tells 
the program what new data are expected for 
the next case: 

INCODE == 1 means a new title card and 
new SDATA are expected. 

INCODE ^ 2 means new film thickness 
data are expected. 

INCODE - 3 means new PDATA are ex- 
pected. 




TABLE IV. - DATA SHEET FOR SAMPLE PROBLEM 


Card column 


1-6 


7-12 


13-18 


19-24 : 


25-30 


31-36 


37-42 


43-48 


49-54 


55-60 


61-66 


67-72 


73-80 


SAMPLE PROBLEM - IMTERMATIONAL UNITS 
iSDATA R1IN= ,082931, R2IN= .084201, ROIFIN=*0., WIDTH=0., M0LWT*28. 966 , 
CP=I.3047788E3, GAMMA=1.4, MU=0. , SPEE0=0., CAPV^60.96, XLAH=1« , XTURB=0.25, 
C0NLAM=24., C0NTRB-.079, RELAM=2300., RETURB=3000. , PWRSKP=F, PRSSKP=F, 
NRMSKP=F, PLTSKP=8*F, NOSI=F, SKPH=F $ 

18 

.00000254 .00000508 .00000762 .00001016 .00001270 .00001524 

.00001778 .00002032 .00002286 .00002540 .00002794 .00003048 

.00003302 .00003556 .00003810 .00004064 .00004318 .00004572 

iPDATA POIN=448158.57, P3IN= 103421.21, PRIN*0. , TO IN»310. 92778 , LOSS=l., 
INCODE*! $ 

SAMPLE PROBLEM - U. S. CUSTOMARY UNITS 
iSDATA RlIN=3.265t R2IN=3.315. RDIFIN*0., WTDTH=0., M0LWT*28. 966 , CP*. 24, 
GAMMA*1.4, MU*0., SPEED=0., CAPV=200., XLAM*1. , XTURB*0.25, C0NLAM=24., 
C0NTR8=.079, RELAM*2300., RETURB=3000., PWRSKP*F, PRSSKP*F, NRMSKP*F, 


PLTSKP*8*F 

18 

, N0SI=T, 

SKPH=F % 




.COOi 

• 0002 

.0003 

.0004 

.0005 

.0006 

.0007 

• 0008 

• 0009 

.0010 

• 0011 

.0012 

.0013 

• 0014 

.0015 

.0016 

.0017 

• 0018 


iPDATA POIN*65., P3IN*15., PRIN=0., T0IN*100., L0SS=1., INCODE*! $ 





TABLE V. - STEPS IN SOLUTION OF 
EQUATION (C3) BY METHOD 


OF ITERATION 


Guessed value 
of M 

G{M) 

Are M and G(M) 
equal? 

0. 1 

0. 3841 

No 

.3841 

.4926 

No 

. 4926 

. 5078 

No 

. 5078 

.5087 

No 

. 5087 

. 5087 

Yes 


TABLE VI. - PARAMETERS IN COMMON BLOCK /CONSTS/ 


Word 

number 

Symbol 

Description 

1 

y 

Ratio of specific heats 

2 

AR 

Flow length 

3 


Exponent in friction factor - Reynolds number 
relation for laminar flow 

4 


Constant in friction factor - Reynolds number 
relation for laminar flow 

5 

"t 

Exponent in friction factor - Reynolds number 
relation for turbulent flow 

6 


Constant in friction factor - Reynolds number 
relation for turbulent flow 

7 

'I'o 

Sealed gas reservoir temperature 

8 

Pq 

Sealed gas reservoir pressure 

9 

P 3 

Ambient pressure 

10 

Ptol 

Pressure tolerance to stop iteration on exit pres- 
sure for subcritical flow 

11 

S 

Gas constant 

12 

^L 

Entrance velocity loss coefficient 

13 

(Re); 

Upper limit on Re for laminar flow 

14 

(Re)t 

Lower limit on Re for turbulent flow 

15 

^0 

Sealed gas reservoir viscosity 

16 

T* 

Temperature at point of choking 





TABLE Vn. - PARAMETERS IN COMMON BLOCK /ARRAYS/ 


Array 

number 

Symbol 

Array di- 
mension 

Description 

1 

X 

(101) 

Distance across face of seal 

2 

P 

(20, 101) 

Pressure 

3 

M 

(20, 11) 

Mach number 

4 

u 

(20, 11) 

Leakage flow velocity (x-direction) 

5 

T 

(20, 11) 

Temperature 

6 

p 

(20, 11) 

Density 

7 

B 

(20, 11) 

Friction parameter, B = 4f(L - x)/D 

8 

h 

(20) 

Film thickness (first element is h*) 

9 


(20) 

Dimensionless center of pressure 

10 

L 

(20) 

Distance from entrance to point of choking 

11 

p* 

(20) 

Choking pressure 

12 

Re 

(20) 

Reynolds number 

13 

f 

(20) 

Mean friction factor 

14 

— 

(20) 

Power 

15 

F 

(20) 

Force 


TABLE Vm. - PARAMETERS IN COMMON BLOCK /CONVRT/ 


Word 
number 
in array 
CONV 

Variable to which 
constant applies 

Word 
number 
in array 
CONV 

Variable to which 
constant applies 

1 

Temperature 

7 

Reynolds number 

2 

Sutherland’s law 

8 

Power 

3 

Sutherland’s law 

9 

Temperature difference 

4 

Speed 

10 

Torque 

5 

Mass flow rate 

11 

Velocity 

6 

Standard volume flow 

12 

Density 





SPRING 


DETAIL OF SEAL RING FACE 



Figure 1, - Pressure-balanced face seal with a gas bearing added for axial film stiffness. 



Figure 2. - Pressure-balanced face seal with no axial film stiffness. 










Figure 5. - Block diagram of solution scheme for finding choking film thickness. 
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Figure 8. - Diagram of subroutine calls. 
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